



Introduction 


SCOLE stands for "Spacecraft Control Laboratory 
Experiment". The objective of the SCOLE Program is to 
provide an example configuration and control objectives 
which enables direct comparison of different techniques in 
modeling, systems identification and control. The SCOLE 
Design Challenge" was formulated in 1983 by L. W. Taylor and 
A. V. Balakrishnan. The details of this challenge are reprinted 
at the end of this document. 

Annual SCOLE Workshops have been held for specialists 
to share and compare their research results. This proceedings 
is a compilation of the material presented at the 4th 
Workshop held at the USAF Academy on November 16, 1987. 
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PARAMETER IDENTIFICATION USING EXPERIMENTAL MODAL DATA 
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The motion of a distributed structure with zero damping is governed 
by a PDE of the form 
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Discretization (cont'd) 

Rayleigh-Ritz Type Parameter Expansion: m(P) 




OR****' 


r««ui> '» 




7 . ; 

•I 


-J* 



OJ 


i 


"Q 




O 












<D 

f,* • 




V 


f— 

<d 



CD 

»— 

, ’ 


XJ 

3 

✓O 


o 

4-> 

<T „ 


2EZ 

co 




O 



r— 

Cl. 



fd 

''•W* 

L* 


3 


* 

» 

>- 

4-> 


1 r* 

J 3 

U 


CO 


<£ 

o 

r d 

CD 

s./ 

X i 


!U 


o 

/i 

x e 


X 

<- I! 


S- 


<'\ 


X < 

II 

O ! 

i c c— o ;i 

:x 

- i 

s- 


u 

.'- 


II 

o 


M 


X i 

c 1 

J- 

O 



X J 

2 : 

O) 

+~> 

2 : 

4. 

> 

S *— ^ 

C\J u 

CM o 

*«— 

CTI 

3 

3 

-M 




rd 




v- 




CD 




4~> 




•— < 





CD 

X 





jz 

4 






3 






4-> 




f— 

+-> 

n3 


i— 



rC 

c: 


»— i 


4- 

Jd 



-C 


o 


X 


o 


3 


CD 


QQ 



XT 

4 




* 

u 

3 


• 


• 

3 

GO 


• 


• 

o') 

fd 


• 




CD 



H— 

CM 

E 

E 


t-H 


O 

aj 



o 

4— 

3 

-C 

CD 


CQ 

3 


4-J 

-C 








cr 

• 

o 

CD 



o 

* 

3 

+-> 

o 


a 

♦ 

- 

fO 

4-> 





X7 



« 

CM 

!l 

Cl 

CD 


« 

3 


3 

cn 


* 




4 



f-H 

O 

CD 

CD 


CM 

3 

Ol? 

£ 

> 

• 

O 




C 

CD 

0 


o 

«*• 

o 

4 


li 

3 i 

4. 

u 

3 

r— < 



O 


4-> 

o 

3 ? 


QCJ 

GO 

<J 

o 


• « 


CD 

3 



to 

«« 

•r— 

4 


# « 

CD 

4 

u 

.o 

it 

GO 

•r— ■ 

O 

c 

GO 


CD 

U 

3 

CD 


o 

«f— 

c 


— > 

i — 

OU 

4) 

CD 

GO 

cr 

<d 


SZ 

3 

4 

CD 

3 


QJ 

cr 

CD 

4 

4-> 

9 • 

3 

a 

4-> 

Gi- 

o 

to 

cr 

s- 

CD 


fd 

4. 

CD 

u. 

E 

f'— 


CD 

s- 


fd 

<d 

QJ 

-M 

Lu 

r— 

4 

4 

x: 

CD 


fd 

<d 

3 

*4 

cr 

r— 

4> 

D. 

-4 


r 3 

fd 

3 


<d 

4- 


s- 

+-> 

QJ 

c 

o 

<d 

3 

fd 

-C 



ry 


21 

-4> 

r— * 

CO 


fd 



rd 

CD 

-a 


T3 

o> 

CJ 

m r~ 

CD 


QJ 

c 

*r~ 

u 

4-> 

T3 

4-> 

•f 

X> 

c 

fU 

CD 

fd 


QJ 

<D 

r— 

4. 

r— 

CD 

4 

13 

3 

3 

3 

u 

O 

Or 

-M 

GO 

4-> 


CD 

CD 

CO 

<d 

GO 

o 

sz 

4 

O 

CD 

o 

4- 

•W 

4- 

a. 

2E 

Q- 


6 


__ pA0& ® 


I 



Iterative Approach (cont'd) 
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Iterative approach (c on t 1 d) 
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SCOLE MODEL: 
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We used an n = 4 degree of freedom model: 
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Indirect approach Is used to obtain physical parameters 
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INFINITE-DIMENSIONAL APPROACH TO SYSTEM IDENTIFICATION 
OF SPACE CONTROL LABORATORY EXPERIMENT (SCOLE) 


S. A. Hossain and K. Y. Lee 

Department of Electrical Engineering 
The Pennsylvania State University 
University Park, Pennsylvania 16802 


ABSTRACT 

The identification of a unique set of system parameters in large space struc- 
tures poses a significant new problem in control technology. This paper presents 
an infinite-dimensional identification scheme to determine system parameters in 
large flexible structures in space. The method retains the distributed nature of 
the structure throughout the development of the algorithm and a finite-element 
approximation is used only to implement the algorithm. This approach elimi- 
nates many problems associated with model truncation used in other methods of 
identification. The identification problem is formulated in Hilbert space and an 
optimal control technique is used to minimize weighted least squares of error be- 
tween the actual and the model data. A variational approach is used to solve the 
problem. A costate equation, gradients of parameter variations and conditions for 
optimal estimates are obtained. Computer simulation studies are conducted us- 
ing a shuttle-attached antenna configuration, more popularly known as the Space 
Control Laboratory Experiment (SCOLE) as an example. Numerical results show 
a close match between the estimated and true values of the parameters. 
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DISTRIBUTED PARAMETER IDENTIFICATION 


TWO APPROACHES « 

F1NITE-DIMESIONAL METHOD 
INFINITE-DIMENSIONAL METHOD 



Table 1 


Literature Surveyed on the Parameter Estimation of Large Space Structures. 


Reference 

Approach 

Wells and Spalding (1977) [2] 

A finite-dimensional design approach 

Tung (1981) [3] 

where the structural model is trunc- 

Balas and Lilly (1981) [4] 

ated and the estimator is designed 

Balas (1981) [5] 

based on the reduced-order model. 

Lee and Bitter (1981) |6| 


■ 

Banks (1982) (7] 


Hendricks et al (1982) [8] 


Hendricks et al (1984) [9] 


Banks and Rosen (1984) [10] 


Rajaram and Junkins (1985) [ 1 1 ] 


Lee, Walker and Hossain (1985) [12] 


Lee (1986) [l3j 


Spalding (1976) [14] 

An infinite-dimensional design appr- 

Burns and Cliff ( 1977) ( 15] 

oach where the PDE model is retained 

Sun and Juang (1982) [16] 

as long as possible and truncation is 

Lee (1986) [13] 

carried out only after the estimation 


algorithm is developed. 
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Real World I Conceptual World 
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SYSTEM IDENTIFICATION SCHEME 




A Distributed System Model 


3^ \l dxi 

m(z)-^y(x,f) + D 0 — {x,t) + A 0 u(x,t) 


x € fi, t € (0, T], 


F{x,t) = Fg{x,t) + F c (i,0 + F D (x,t), 


M 

Fc(x,t) = B 0 f = Y,b l (x)f l {t), 

t= 1 


jV 

F B {x,t) = B b g = 

i=i 


y = Cqu + Eout, 


F{x,t), 
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Basic Problem Formulation 


,(•) - i i (y - z) T R[t){y - z)dt, 


where z is the measurement of output vector y given as 

z{x,t) = y[x,t) + e{x,t) 

with a measurement error e{x,t). Also, it is defined that 

(y - z) T R{t)(y - z) = J \y{x,t) - z{x,t)\ T R{x,t)[y{x,t) - z{x,t)}dx, 


22 



Infinite-Dimensional Formulation 


|U(o + + = mm 

dxi 

ua 2 (0,T;V), *eLi(0,T;H), 

ot 

where f(t) is given in 1/2(0, and the initial conditions are 

u( 0 ) = u 0 , u 0 given in V, 


and 



= ui, 


Ui given in H. 


The output function is 

y(0 = Cuft), 

The identification problem can now be formulated as an abstract problem of deter- 
mining the parameter vector q‘ (x) 6 Q that minimizes 

T 

JM = ^ f l»(<) - *(i)l T *(0|y(0 - »(0I*. 

where z{t) is the observed data belonging to Y 
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Development of Infinite-Dimensional Identfication Algorithm 


THEOREM : Given a state equation (18) with initial conditions given by Eq. (19) and 
the cost function by Eq. (22) with y(t) satisfying Eq. (21), then the optimal parameter 
vector q * satisfies the state equations (18)-(19) and the foilwing system of equations : 

£p(0 - D‘ 4p(0 + A-p[t) = -ic T R(Cu - z), (23) 

at * at l 

with the final conditions 

p(T) = j t P{T) = 0, (24) 

and the first variation of an augmented cost functional is 

6J a — [ p T ^-\D^ + Au - Bf] 6qdt = 0, (25) 

Jo d( 1 dt 

where p(t) is a costate variable also belonging to the Hilbert space V . 

PROOF : By combining Eqs. (18) and (22) an augmented cost functional can be defined 
as 

T* 

= ^J mo - z(t)\ r jnt)m - z(t)\dt 
+ j f p(t) T l^u(() + Dj t + A«(0 - B/(t)] <4!. (26) 
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Parameter Identification of Vibrating Beams 


Case I : A Simply-Supported Beam 


d-u d A u 

pA dti +E1 d^ = xe(o,i],t >o, 


d 2 

“( x ’0 = °. xed[o,L\, t > o, 
u(x,°) = _ u( x 1 0U = O, . x€[0,Lj. 


»(0 = «(§*<)• 
2 (0 = y(0 + e(t). 


J = j 0 \y ~ Z \ T R \y - z \ dt - 

25 



DISPLACEMENT 


m s 6 7 kg/m 

El * 23000 N/m 2 


STEP LOAD 




rr ig. 2 (a) Simply-supprted beam with step load, 
(b) resultant displacements. 
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d 2 U 
dt 2 


Eld 4 u 1 r/ L w , . 
ma? + m <(l "2 )/(,) ' 


or 


a 2 u 

at 2 




where 


El 

9i = — , 

m 


1. 


92 = 


m 


and the parameter vector is defined by 9 — [9 i,92) T 


a 2 P a 4 P R. , . L, 

W ~ ~ qi dx* + r^ u ^ z 2 ’ 


x e [o,L], t € [0,7), 


p(^^) = ^p(x,f)| t=r = °, ie[o,L]. 


p(*.0 


£ t p(x,o = °, 


x 6 a[0,L],t 6 [0,7). 


5J a _ f T d 2 pd 2 u 

s qi ‘ y 0 ax 2 ax 2 


<592 


= - J n pH x ~ ^)/(0*- 


Thus, parameters can be updated by the steepest descent algorithm 


k+1 it ( 

(j- 


* = 1,2- 
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Table 2 


Performance Data for Case I 


Iteration 

9i 

92 

m 

El 

/ 0 T error 2 dt 

1 

301.58 

0.0158 

63.00 

19000.00 

0.20178E-01 

2 

311.68 

0.0154 

64.84 

20207.72 

0.16640E-01 

3 

325.65 

0.0151 

66.27 

21582.87 

0.71528E-02 

4 

338.48 

0.0149 

66.91 

22647.68 “ 

0. 57208 E-03 

5 

342.11 

0.0149 

66.98 

22913.64 

0.34507E-04 

6 

342.96 

' 0.0149 

66.99 

22973.13 

0.26524E-04 

7 

343.19 

0.0149 

66.99 

22989.35 

0.24364E-04 

True 






values 

343.28 

0.0149 

67.00 

23000.00 
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Case II : A Cantileverd Beam 


d 2 u / — — d 3 u 

m-r~x - 2 ZvmEI- 


dt 2 


dx 2 dt 


■Ei a ^ = b(z)m, 


x € [0, L\, t > 0, 


where £ is a damping coefficient. 


y(0 = u{L,t). 


The boundary conditions for a cantilevered beam are : 


d 2 

«(< M ) = ^ u ( z ’0 


= o, 


i=0 


d 2 

dx 2 




\z — L 


dx 3 


u(x,t) 


= 0, 


t > 0, 
t > 0. 


\z~L 


The beam is initially at rest and hence the initial conditions are 


U ( I,G ) = at U ( I,f )l«=o = °’ xe (°’ L l- 


du(x, t) 

dx 


\t-o 


dx 2 


u(x,t) 


= 0 , 


€ [0, L\. 


i t=o 


d 2 u d z u d 4 u 

d?- q3 d^3, 


where 


9i = 


El 

1 

m 

i 


92 - — , 

m 

93 = 2 

V m 


and the parameter vector is defined by q = [gi,< 7 2 , 93 ] T 
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DISPLACEMENT 


STEP LOAD 



o 



TIME (ms) 
(b) 


Fig. 3 (a) Cantilevered beam with a step load, 
(b) resultant displacements. 
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Table 3 


Performance Data for Case II 


Iteration 

9i 

92 

93 

m 

El £ 

Jq error 2 dt 

1 

301.59 

0.0159 

2.08 

63.00 

19000.00 0.060 

0.7843E-01 

2 

318.01 

0.0155 

2.76 

64.50 

20512.92 0.077 

0.1664 E-01 

3 

326.20 

0.0153 

3.06 

65.26 

21286.34 0.085 

0.9432E-02 

4 

331.30 

0.0152 

3.23 

65.73 

21775.47 0.089 

0.4498E-02 

5 

334.72 

0.0151 

3.34 

66.05 

22107.25 0.091 

0.2284 E-02 

6 

337.10 

0.0151 

3.43 

66.27 

22339.85 0.093 

0.1238E-02 

7 

338.83 

0.01.51 

3.49 

66.43 

22510.09 0.095 

0.6649E-03 

8 

340.07 

0.0150 

3.53 

66.55 

22632.61 0.096 

0.3796E-03 

9 

340.96 

0.0150 

3.55 

66.65 

22721.78 0.096 

0.2199E-03 

10 

341.56 

0.0150 

3.56 

66.71 

22783.93 0.096 

0.1519E-03 

11 

342.06- 

0.0150 

3.59 

66.76 

22835.96 0.067 

0.8667E-04 

12 

342.44 

0.01497 

3.62 

66.80 

22874.19 0.098 

0.4961 E-04 

13 

342.73 

0.01496 

3.64 

66.83 

22903.95 0.098 

0.1165E-04 

True 







values 

343.28 

0.0149 

3.70 

67.00 

23000.00 0.1 
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Case III : A Simply-Supported Beam with Spatially Variable Parameter 


d 2 u 


dt 2 




6(x)/(«), 


x € [0,L], t > 0, 


where q = El{x). 


d 2 P 
dt 2 


dx 2 




xG [0,1], t 6 [o, r), 


g. w _ r gpg± A . 

s<; ax 1 ax 2 


32 



SECTION 1 


STEP LOAD 


SECTION 2 


E 




|- L/2 

SECTION 1 
m = 67 kg/m 
El = 23000 N/m 2 


(a ) 


— L/2 

SECTION 2 

m = 67 Ic g/m 

El = 40000 N/m 2 


>— 

o 

z 

i 

LU 

UJ 


iO 

uj £ 

<0 

o • 

to 

< *7 


-j ! 

• 

a. £ 

UJ o 

CO 

-J R 

•— 

< c 

a 

o - 

CO — 




TIME (■•) 



Fig- 4 (a) Simply-supported beam with spatially variable flexible rigidity, 

(b) resultant displacements at L/4, 

(c) resultant displcements at 3L/4. 




Table 4 


Performance Data for Case III 


Iteration 

q[se c. 1) 

q( sec. 2) 

/ o r error 2 dt 

1 

21000.00 

39000.00 

0.76863E-02 

2 

21686.15 

39286.47 

0.40601E-02 

3 

22261.30 

39521.45 

0.14636E-02 

4 

22630.66 

39688.41 

0.41442E-03 

5 

22814.49 

39759.51 

0.13629E-03 

6 

2291 1 .04 

39793.95 

0.51 151E-04 

7 

22968.10 

39815.07 

0.1 1209E-04 

True 

values 

23000.00 

40000.00 
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^'8- 1 The shuttle j antenna configuration of the spacecraft 
control laboratory experiment ( SCOLE ). 




YAW 

Z 


t 



PITCH 


Fig. 2 Drawing showing the direction of “roll" bending, 
“pitch" bending and “yaw” twisting. 
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The SCOLE Mathematical Model 


A. Dynamic Equations 


Roll Beam Bending Equation in v-z Plane 




n= 1 ' 


(t)5(s 3 n ) + 9t,n[t) $ { 3 


t > 0 , -CO < 3 < 00, 0 <3 n < L. 



Pitch Beam Bending Equation in x-z Plane 


pA-j^ - 2 id\fpAEld 


d 3 ug 

da 2 dt 


+ 


EI 9 


da 4 


— 3 n) + — -*n)^ > 

t > 0, — oo < a < oo, 0<a n <L. 


Yaw Beam Torsion Equation for z-Axis 


4 

= — 3 *)> 
n=l 

t > 0, —oo < a < oo, 0 < s n < L. 
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B. Forcing Functions 


The forcing functions on the right side of each equation ire dependent on boundary 
conditions and proof-mass actuators. 


Forces at s = 31 = 0 ( shuttle body forces ) 


The forces at si =0 involves the shears at that point which are equal to the shuttle 
mass mi times the corresponding component of acceleration. 


ft »i (0 Qt% 0 » 

d 2 

= -mi— u e (0,t). 


Forces at s = sa = L [ reflector body forces ) 

d 2 d 2 

ft 1^(0 ~TT%4 0 _ t) — F y , 

d 2 d 2 

/ m (0 = - m 4ft2 u o{L,t) + m ^y—V. i> {L,t) + F z , 


where is the reflector mass, (r s ,ry) is center of reflector mass from the beam tip at 
s = L, and F x and F y are the applied forces at the center of the reflector mass. 


Forces at s = 37 f proof-mass actuator forces ) 

d 2 d 2 

ft, 2(0 = ~ m 2 ^u*(«S 2,0 -)- m 2 — &, t>t2 , 
. . , d 2 . d 2 

fe, 2 [t) = -m 2 —u 9 {3 2 ,t) + m 2 — A* >2 , 


where A and m denote displacement and mass of the proof-mass actuator. 

Forces at s = sa f proof-mass actuator force ) 

d 2 d 2 

ft, a(0 = “*(*3,0 + m 3 — A* 3 , 

d 2 d 2 

f«, 3(0 = ->713^2 U 9 (s 3 ,f) + m 3 — A*, 3 . 
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C. Momenta 


Moments at s = 0 f shuttle body moments ) 

1 9 +' 1 \ 

I 99,i I = — [hwi + wi ® 4- Mi(t) + Mo{t), 

\H, i) 


where I\ is the moment of inertia of the shuttle body, Mi(t) and Mx?(t) are control and 
disturbance moments, respectively, applied to the shuttle body, and ® denotes the vector 
product. 

Moments at s = L ( reflector body moments ) 


99 , A j = “ f-^A^A + w 4 ® I4W4 - M 4 (t ) + r ® F 4 [t)\ “ m * r ® -fftT' 

Ha / ' 7 

where M4 and F4 are the control moment and force applied at the reflector center of the 
ma as and £4 is the coordinates of the beam tip. 

Also, J 4 is the moment of inertia of the reflector, and I4 is that with respect to the beam 
tip given by 

( r l ~ r * r v 0 \ 

I4 = I4 + m 4 -r s r v r\ 0 

V 0 0 r l + r lJ 
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Abstract Formulation of the SCOLE Problem 


M 0 f (t) + A 0 r{t ) + B 0 F{t) + K 0 {r{t) 2 ) = 0 , 


where Mq is the 17 x 17 matrix specified by 


*3 

*3 



*4 

*5 

*T 

*• 

*0 
*10 
*11 
*13 
*13 
*14 
* 1 * 
*10 
*1T V 


*3 

0 

0 


*3 

0 

0 

pA 


*4 


mi 

0 


*6 X<j X 7 


0 

mi 

nu 0 

0 m* 

0 0 

0 0 

0 0 

0 0 

0 0 

r * — 


*• *o *10 *11 


0 0 0 0 

0 0 0 0 

0 

h 0 

0 

0 0 0 

0 0 0 

0 0 0 


*13 *13 *14 *is *ia *1 


0 m«r, 

0 -m^r,. 

0 0 

0 0 

0 0 

A 

mj 0 

0 rr%2 

m 3 0 
0 m 3 
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Table 2 


Performance Data for Case I : Nonlinear SCOLE Model 


Iteration 

7 1 

72 

m 

El 

/ Q r error 2 dt 

1 

4.2353E+08 

11.765 

0.0850 

3.6000E+07 

0.1341E-04 

2 

4.2531E+08 

11.067 

0.0904 

3.8427E+07 

0.1896E-05 

3 

4.2612E+08 

10.809 

0.0925 

3.9422E+07 

0.3038E-06 

4 

4.2743E+08 

10.705 

0.0934 

3.9929E+07 

0.9131E-07 

5 

4.2704E+08 

10.665 

0.0938 

4.0041E+07 

0.8031E-07 

6 

4.2617E+08 

10.637 

0.0940 

4.0064E+07 

0.6646E-07 

7 

4.2508E+08 

10.611 

0.0942 

4.0059E+07 

0.4938E-07 

8 

4.2390E+08 

10.589 

0.0944 

4.0029E+07 

0.3237E-07 

9 

4.2277E+08 

10.574 

0.0945 

3.9982E+07 

0.2032E-07 

10 

4.2174E+08 

10.561 

0.0947 

3.9934E+07 

0.1465E-07 

11 

4.2072E+08 

10.546 

0.0948 

3.9894E+07 

0.1332E-07 

12 

4.1959E+08 

10.525 

0.0950 

3.9867E+07 

0.1365E-07 

13 

4.1819E+08 

10.495 

0.0953 

3.9845E+07 

0.1676E-07 

14 

4.1783E+08 

10.471 

0.0955 

3.9903E+07 

0.7113E-08 


True 

values 4.1858E+08 10.465 0.0956 4.0000E+07.00 
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Table 3 


Performance Data for Case I : Linearized SCOLE Model 


Iteration 

Si 

?2 

m 

El 

Jq error 2 dt 

1 

4.2353E+08 

11.765 

0.0850 

3.6000E+07 

0.1341E-04 

2 

4.2531E+08 

11.068 

0.0904 

3.8427E+07 

0.1894E-05 

3 

4.2613E+08 

10.809 

0.0925 

3.9422E+07 

0.305 IE-06 

4 

4.2744E+08 

10.704 

0.0934 

3.9929E+07 

0.9142E-07 

5 

4.2705E+08 

10.665 

0.0938 

4.0041E+07 

0.8040E-07 

6 

4.2618E+08 

10.637 

0.0940 

4.0064E+07 

0.6660E-07 

7 

4.2509E+08 

10.611 

0.0942 

4.0060E+07 

0.4955E-07 

8 

4.2392E+08 

10.590 

0.0944 

4.0030E+07 

0.3255E-07 

9 

4.2277E+08 

10.574 

0.0945 

3.9982E+07 

0.2032E-07 

10 

4.2175E+08 

10.561 

0.0947 

3.9933E+07 

0.1478E-07 

11 

4.2066E+08 

10.546 

0.0948 

3.9888E+07 

0.1378E^-07 

12 

4.1957E+08 

10.524 

0.0950 

3.9567E+07 

0.1361E-07 

13 

4.1816E+08 

10.495 

0.0953 

3.9846E+07 

0.1667E-07 

14 

4.1773E+08 

10.471 

0.0955 

3.9893E+07 

0.8371E-08 

True 






values 

4.1858E+08 

10.465 

0.0956 

4.0000E+07.00 
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FLEXURAL RIGIDITY 



0 2 4 6 810 121 4 

ITERATION 
( a ) 



0 2 4 6 8 10 12 14 

ITERATION 

(b) 


Fig. 5 Convergence of parameters for the SCOLE problem in 
Case I using nonlinear model. 

(a) For pA / unit length. 

(b) For flexural rigidity, El. 
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I 



1 



TIME (»•) 



1000 


Fig- 4 Roll, pitch and yaw displacements when a damping 
lactor of 0.003 is added to the system. 
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Table 4 


Performance Data for Case II 


Iteration 

El 

pA 

pI'P 

e 

1 

32000000.00 

0.0750 

1.1000 

0.00250 

2 

37582208.00 

0.0881 

0.9690 

0.00250 

3 

39607652.00 

0.0928 

0.9428 

0.00250 

4 

40157244.00 

0.0941 

0.9362 

0.00250 

5 

40372424.00 

0.0946 

0.9247 

0.00280 

6 

40338068.00 

0.0945 

0.9183 

0.00288 

7 

40309128.00 

0.0945 

0.9127 

0.00293 

8 

40282296.00 

0.0943 

0.9088 

0.00289 

9 

40260868.00 

0.0943 

0.9047 

0.00287 

10 

40243016.00 

0.0942 

0.9008 

0.00285 

11 

40226244.00 

0.0942 

0.8983 

0.00281 

12 

40202376.00 

0.0942 

0.8979 

0.00286 

13 

40181576.00 

0.0942 

0.8977 

0.00292 

• 

14 

40162324.00 

0.0941 

0.8974 

0.00299 

15 

40149520.00 

0.0941 

0.8973 

0.00298 

16 

40140040.00 

0.0941 

0.8973 

0.00295 

True 

values 

40000000.00 

0.0956 

0.9089 

0.003 
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UNIT LENGTH MASS / UNIT LENGTH 


0.095 




0.085 


0.075 


r * 


j i i i 


J L 


0 2 4 6 8 10 12 14 

ITERATION 

(a) 


\ 

* 


1.10 
1 .00 
0.90 


h* 


* * 

*_*_*_* *-*_*+-*--* # 


I I I I I L 1 1 m 

0 24 6 8 10 12 14 

ITERATION 

(b) 

Fig. 7 Convergence of error criterion for Case II. 
* 

(a) For pA / unit length. 

(b) For pl^ /unit length. 
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CONCLUSION 


Infinite-dimensional identification method presented in this paper shows a significant 
promise in the parameter estimation of flexible structures with great potentials for appli- 
cations to LSS’s. The basic approach is the abstract formultion of the system dynamics 
in function spaces and then applying optimal control theory to adjust system parame- 
ters so that the error between actual and model data is minimized. The use of partial 
differential equation for the purpose of estimation elimainates many problems associated 
with model trunction in the finite dimensional approach. Based on partial differential 
equation models and a quadratic performance index an algorithm to estimate the opti- 
mal parameters has been developed. The numerical results show the effectiveness of the 
algorithm in estimating parameters of the flexible beam in the SCOLE problem. The 
results show fairly good match between the model and the estimated parameters. How- 
ever, as the number of parameters to be identified increases it becomes increasingly time 
consuming and difficult to solve. Also, due to model mismatch, slightly less accuracies 
are expected if experimental measurement data from physical beam were used. 
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Abstract 

We introduce a class of nonlinear damping models with application to flexible flight 
structures characterized by low damping. We are able to obtain approximate solutions of 
engineering interest for our model using the classical “averaging” technique of Krylov and 
Bogoliubov. The results should be considered preliminary pending further investigation. 
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1. Introduction 


The problem of characterizing the damping mechanism in flexible structures has 
received renewed attention in recent years in connection with the need to stabilize flexible 
flight structures such as antennas deployed in space. The damping models even when 
simplified to be linear appear to lead to rather complex mathematics if the structure is 
described by partial differential equations and much progress has been made (the analyticity 
of the generated semigroup has been shown to be essential). But experimental evidence as 
in SCOLE [ 1 ] seems to support the need for nonlinear models — the decrement is much 
smaller than predicted by linear models. Some of the difficulty inherent in handling 
nonlinear models is offset by the fact that damping, whatever its nature, is still small. This 
opens up in particular the feasibility of obtaining approximate solutions using the classical 
averaging method of Krylov-Bogoliubov [2], 

In this paper we study a class of nonlinear models and approximate the response by 
the Krylov-Bogoliubov technique. We use a modal expansion and neglect off-diagonal 
terms. The emphasis is on useful engineering solutions rather than abstract mathematics. 

We begin in Section 2 with the primary nonlinear damping model for the simplest 
system — the one-dimenisonal or single-mode case. We emphasize in particular one 
feature that emerges, viz., the potential lack of identifiability from response data. In Section 
3 we generalize to the multi- (non-finite-) dimensional case. In Section 4 we show the 
relevance of the Krylov-Bogoliubov technique for approximating solutions to nonlinear 
boundary feedback. We may mention that there is much work — even classical in nature 
— on nonlinear oscillations such as the nonlinear pendulum where the spring constant is no 
longer linear; however, relatively little attention appears to have been paid to the small 
nonlinear damping term case. 
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2. Single-mode Example 

To illustrate ideas, let us begin with a one-dimensional (single-mode) example: 

x(t) + zD(x, x) + to 2 x(f) = 0 (2.1) 

where the dots indicate time-derivatives, as usual. We assume that: 

D(x, x)x > 0 (2.2) 

so that for E{t), the energy 

E(t) = f (i(r) 2 + co 2 x(f) 2 ) (2.3) 

we have 

J^-£(f) = -eD(x, x)x < 0 (2.4) 

satisfying the energy nonincrease requirement. The particular choice for D(x, x) we 
shall make is: 

D(x, x) = 2co£i + yx 2 " 1 |x| a + (2.5) 

where m, n are nonnegative integers, 

0 < a, (3 and 0 < a + P < 1; 0 < ^ < 1, 0 < y < 1 . 

For small enough e we may apply the averaging method of Krylov-Bogoliubov [2, 5]. 
Thus, we write for the approximate solution: 

x(f) = a{t) sin ((Of + <J>(f)) (2.6) 

where the amplitude function a(t) and the phase function <|>(r) are slowly varying over 
the period T = 2kJ (0. According to the K-B approximation [2]: 

f, - -I'M-) <- 2 - T > 

where 

l r 2n 

Kq ( a) = ■=— } D(a sin <J>, a( o cos <|)) cos <() d§ (2.8) 

K o 

S£ 



-3- 


and 


d§ _ £ 

dt (0 a 


Po(a) 


2 n 


1 f 

Po(<* ) = 2 ^ J P*( a sin < t ) > aa) cos 40 sin <(> d<j> 


Now we can readily calculate that for our choice, because of (2.2), 


Po(a) = 0 


K 0 (a) = 2(0 


1 2 n 

2 ~ J co cos 2 <|> d<$> 


1+P 2m+2n+ l+a+| 

+ (0 Y a 


2 n 


2 ^ J sin 2m <f> cos 2n + 2 <(> (sin <t>| a |cos tyf 


. _ 2 r 2m + 2n+l + a + P 2 « + 1 + B 

- or £ a + a yco 




where 


2n 


j AH 

(J. = 2n J sin 2m <J> cos 2n <() cos 2 <(> |sin <)>| a |cos <}>|^ d§ 

7 o 

and is a constant less than 1/2. Hence letting p = 2/n + 2n, we have 

da * p+l+a+B 2n + B s 

Ji - -e(coCa + <r p to p y^) 

We may set e = 1 without loss of generality since we may absorb it into £ 
Then 

*(') 


r = - 


da 


a<0) 


77T p + 1 + a + B 2n + B 

+ a CO 


YM' 


yielding 


a(r) = a(0)e~ 7 “ C [l + a(0) p+a + P 0) 2 "" 1 + P ^ (1 - e~“° ;(p+ a + P) )j 
We can readily verify that for £ = 0, we have 

a(r) = a(0)[ 1 + a(0) p+ “ + ^co 2n + Py|i(/> + a + |3)r] F+a+T 


" p + a + 3 


F 7 


(2.9) 

( 2 . 10 ) 


d< t>j 

( 2 . 11 ) 

( 2 . 12 ) 

and y. 

(2.13) 

(2.14) 


(2.15) 



The case y = 0 is even more obvious. One salient fact that emerges immediately from 

(2.14) and (2.15) is that it would be difficult to resolve p + a + p into its components 

from response data, unless we can change (0. 

Note also from (2.15) that the rate of decay is not exponential in t and further the 

decrement over any integral multiple of the period depends on the initial amplitude as well 

as the frequency of oscillation. Finally for inegral k and 

_ Ink 

co 

we have, taking logarithms and setting c = p + a + (3 


log 


4—1 

l CO j 
0(0) 


-2nkC, 


c ,0g 


1 + 


0(O) c CO 2 " -1 + Pym. 

C 


(1 -e 


-271 kt,c 


) 


(2.16) 


For small £ this is well approximated by 

-Inkt. - ^log (1 + 2nkXa(0) c ) , 


where 


X = 


2 n + B - 1 
(0 K 


Y H • 


The slope (as a function of k) 


= -2kC> - 


2nXa(0) c 
c( 1 + 2nkXa(0) c ) 


(2.17) 


and hence the linear damping term is yielded by the asymptotic slope as while 

for small k there is a marked curvature which depends also on the initial amplitude 
a(0). The inital (at k = 0) slope 


_ 2k c > _ 2 ^MQ) C 


is larger (in absolute value). The second derivative being positive, the curve is convex 
— CUP. This is in excellent qualitative agreement with SCOLE damping data: see 
Figure 1 where amplitude is plotted on logarithmic scale (period =* 5 seconds). 
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To get another version of (2.16) we may replace (2.7) by the more exact formula 



gjt+T} - a(r) 

T 

-§K 0 (a(t)) 

(2.19) 

and hence using 

a k - 

a(kT) 

(2.20) 

we would have 

a k + 1 = a k 

- («*) 

(2.21) 

so that 


r er 

l 1 (0 a k ) 


which under our 

“small damping” assumption, may be replaced by 



a t + l 

108 o, • - 

eT k o (a*) 
to a k ’ 

(2.22) 


= — 

e2jc(C + a c k yp<o 2 "" , + P ) . 

(2.23) 


3. Multidimensional Generalization 

Analogous entirely to the one-dimensional case, we may write the general nonlinear 
dynamic equation for flexible structures [2] as 

Mx{t) + 3)(x(t), x(t)) + Ax(t) = 0 (3.1) 

where the state x(t) ranges in a separable (real) Hilbert space X; M is a self- 
adjoint positive definite (with bounded inverse) operator on X onto X; A is a self- 
adjoint nonnegative definite closed linear operator with domain dense in X and with 
compact resolvent; we shall (for simplicity) assume that zero is in the resolvent set of 
A. In the linear case 

$(x(t), x(t)) = Dx(t) (3.2) 

where D is also a self-adjoint nonnegative definite closed linear operator whose 
domain includes that of Va . In the most important case we further specify that 

60 


I 
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D$ k = 2^ k a k M$ k 


where {<(>*} are the AZ-orthonormalized eigenfunctions of A with eigenvalues ( 0 \ such 


A$ k = i£r k M§ k . (3.4) 

Here ^ is the damping ratio. If C* = C an d we have strict proportional damping — 

see [3] for more — D is then essentially the positive square root of A (except for M). 

More generally we require that lim ^ > £ > 0. In the nonlinear analogue of (2.5) we 

k 


<t>jk) = 0 j* k 


and more generally for x, y such that 


where 


y 1 2 .2 y 1 4m + 2 a .2/1+2 + 23 2 
I (0* b k + I a k b k y k < 


a k = l*> 4>*1 ; b k = Ly. <t>*] 


we define: 


ISC*. >). <M = y k *\ m K\ a b l n+ ' \h? + 


where, as before, m and n are nonnegative integers and that 


0 < a , p < 1 ; a + P < 1 ; 0<y 


Note that 


[ 3 )(*, y), y] ^ 0 


for every x and y. Hence if 


E(t) = j{\Ax(t), *(/)] + [Mx(t), ±(r)]} 


we have that 


f t m = -[SDW0,i(0). i(0)I ^ o 


( 3 . 10 ) 


Or, the energy is nonincreasing. Using the modal expansion 
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x(t) = X a k (t) 

we see that for each k 

a k (t) + CO \a k (t) + a k {t) 2 \a k {t)f a k (t) \a k {t)f + 2^ k (O k a k (t) = 0. 

We can therefore invoke the K-B averaging procedure obtaining the approximate 

a k (t) = A k (t) sin (CO k t + <J>*) 




fi . c 2 n + P-1 YtP’/ 

[i + At(0) — ^( 


i - «f 


AtM — A k ( 0) e 
And for ^ = 0, 

i_ 

A k (t) = A k (0 ) l + A*(0) c coJ" +fJ y*^cr] 

where, as before, 

1 2 n 

(X = ^ J sin 2m <{> cos 2n <i> l s in §(* |cos <t>f cos 2 (j) d§ . 
o 

For a + (3 = 0, we can give a kernel representation. Thus 


r 


Z = T(x, y) = I y, <|>, M»,-, x] 2m [<|>-, y] 2n+l 


where 

X y? < °° ; y, > 0 ; 

and for the concrete realization X = Z> 2 (0, L ), the corresponding “kernel” would be 
W(s, o,, .... o 2m ,s,, .... 5 2n+1 ) = X y, ♦,•(«) <t>,(Oj) - ♦ / (o 2m ) ♦ i (si) - 

i 

and 

L L 

z(s) = f ■■■ f W(s, o,, o 2m , s,, s 2n + l ) x(o,) x(a 2 ) ■■■ x(a 2m ) 

0 0 

X y(s t ) ■■■ y(s 2n+i ) do, do 2 - da 2m - ds 2n+ , . 

A plausible model in this case would be to rewrite (3.1) as 

Mx{t) + h(x(t), Dx(t)) + 2 ^Dx(t) + Ax(t) = 0. 



I 


(3.11) 

(3.12) 
solution 

(3.13) 

(3.14) 

(3.15) 


(3.16) 

(3.17) 

(3.18) 
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which will satisfy (3.10), since the y, in (3.16) are nonnegative. In the notation of 
(3), the “roll” equations for example will have the form: 


pAu^t, s) + El^”(t,s) 


2£Vp AEI^ 


d 3 u^(t, j) 

dTWs 


L. U 

- J -J W(s, 


>1- 


>2 m > *1 > 


*2» + l) X “*('* <*l) - °2m) 


d\('. Jl) 

x ~jtw~s 


3 \<*‘ ^2 n + 1 ) 
dtd 2 S 


do l "• d<S 2 m 


ds 2n+\ 


= 0. 


It is clear that we may generalize (3.17) without recourse to modes. The “nonlocal” 
nature of the operator should hardly be surprising, since this is already so in the 
linear case if we want strict proportionality (£* = Q for example. 


4. Application to Nonlinear Boundary Feedback 

In this section we shall apply the K-B averaging technique to obtain approximate 
solution to the response of a flexible structure to nonlinear boundary feedback control. The 
control effort is small so that the K-B approximation is reasonable. We follow [4] for the 
model where the “boundary” is finite-dimensional. Thus we have in the same setting as 
Section 3, but omitting the natural damping term: 

Mx(t) + Bf(B*x(t)) + Ax(t) = 0 (4.1) 

where B means R m onto H and /(•) maps R m into R m and is such that 

[/(m), u] > 0 for u * 0 . (4.2) 

Using the modal expansion as in Section 3: 

x(t) = I a k (t) 4)* 


we obtain 


d k (t) + CO \a k {t) + [/(I dj(t)bj), b k ] = 0 

€3 


(4.3) 
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where 


Taking the approximation 


we see that setting 


that 


B** k = K . 

[/(I aj(t)bj), b k ] = [f{a k {t)b k ), b k \ 
a k (t) = A k (t) sin (G)*r + $*(0) 

f <t >*(') - o 

*o(A*(0) 

5r^ (,) * 


where 


1 2n 

K 0 (a) = 27 c J [/(aco* cos b k ] cos <J> 


To simplify matters further let us take 


m - 1 . 


Then 


We shall take: 


[/(a<n*f> t cos 4»), = b k f(a b k (ti k cos 4>) 


f(u) = X. tan - 1 u 


which is consistent with (4.2). Then (4.3) becomes 


2 2 

a*(0 + co*a*(r) + Xb k tan -1 a k (t) = 0 


Hence 


j 

/Co (a) = J "kb k tan -1 (ab k i o k cos $) cos 4> 
71 o 

- 2^W> + ^ - >1 

A k (t) dA k (t) 

[Vl + a 2 b 2 k (iil - 1 



To solve this, let 

z(r) = 1 + A k (tfb 2 k (£>l 

64 


(4.4) 


(4.5) 


I 
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so that 


<m 

2(V*0 - 1) 


-b\ \ dt . 


(4.6) 


Let 


F(z) = e 1 (VF - 1) , z > 1 . 


Then 


F’(z) >0 for z > 1 


and hence we may define the inverse function 

F{z) = y ; z = F _ 1 (y) . 

Thus (4.6) has the solution: 


(4.7) 


z(0 = F" 1 [F(z(0))e“ 6 * ,t '] , 


(4.8) 




(4.9) 


where 

z(0) = 1 + >4* (0) 2 b\ (61 > 1 


unless A k (0) = 0. Note that 

F~\y) -* 1 


as y -4 0 


and hence z(t) decreases monotonically to 1 and hence the amplitude A k (t) decays to 
zero asymptotically. 

Note that the decay rate depends on the control effort \b\ as well as the initial 
amplitude. Of course we have in (4.4) yet another nonlinear damping model. Following 
(2.22) we have: 


log 



2nX 

2to? 


f Vl + dfojfoj - 


n 


where 


a j = A k(J T ) ; 
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Shuttle (and Reflector) Body 

®J= -Ii “ M i “ M l,BeaiJ 

▼l = ( F 1 + F t,BeamVm, 


»i»T jji ipT 

T, = -idiT, 

Roll (and Fitch) Beam Bending 

paA-C^S^ + EI A =D r «.n 6(s ' s n) 


♦ dt 2 *ds 2 dY « ds 4 fcr l 


Yaw Beam Torsion 


+ e 0 , n ^ ( s-s n ) 


plA + CI^V-GI A = ^8 V n 6(s-s n ) 

*dt 2 *ds 2 dt as 2 n^l 


Beam Elongation 


p aA+c z a^-ea^ Z 


dt 4 


ds 4 



*z,n^ s ~ s n^ 
n= 1 


no 




Error in Static Deflection 






.05 .17 % 

.10 .67 % 

.20 2.7 % 

.30 6.0 % 

.40 10.6 X 

.50 16.4 1 
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Lumpec^Mass^odeJ 

• Exact Static Deflection 

• Approximates Low-Frequency Modes 

• Nonlinear Kinematics 

• Linearized State Space, Modal Model 

• Classical Damping(Working Proportional) 

• Extended to n-Body Network 
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•4 El 2WL 
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0 


Mu = 


0 


4 El 2WL 
L 15 


M u = 


0 

_ 6EI.W* 

0 l 2 10 

6E_I,W* o 
L 2 10 


0 


0 


0 


0 


0 


0 


# Gravity Effect 
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1 2EI 6W* 
L 3 5L 


0 


0 


1 2EI 6W* 
L 3 5L 



0 

0 


EA 

L 



0 

0 
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Asymptotic Approximation 





. Motion Approaches Clamp-Clamped System 
as Mode Number Increases 

• Accuracy Increases with Increasing 

Mode Number 

• Explicit Expressions for Modal Frequencies 

and Mode Shapes 

> First Variation Approximation for Motion of 
End Bodies 


> Singular Perturbation Technique can be used 
to Improve Approximate Solutions 


77 



Comparison of Modal Frequencies 



EXACT 


FINITE ELEMENT 


MOPF no 

IFF 5*6 

IFF 9 

REF. 10 

REF. II 

1 

.278025624 

.278 

.277 

.2740 

2 

.313776751 

.317 

.314 

3??9 

3 

.812326353 

./2b 


.7494 

4 

1.18366347 

1.226 

1.175 

1.244 

5 

2.05047101 

2.069 

2.028 

2.052 

6 

4.75561758 

4.77 

4617 


7 

5.51248431 

552 

5388 


8 

122598619 

124 

11.782 


9 

12.8877037 

130 

12.513 


10 

235359367 

24.2 

14.670 


11 

24.2568205 

24.7 

22.968 


12 

26.4794890 

26.2 

23490 


13 

38.9199260 

454 

37.568 


14 

39.4643489 

459 

38.146 


15 

45 1313668 

56.3 

\ 44.653 


16 

57.90 


\ 45161 



17 57.92 

18 80.72 

19 80.72 
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Comparison of Modal Frequencies 


EI ACT 

MODE NO REF.5&6 

1 .278025624 

2 .313776751 

3 .812326353 

4 1.18366347 

5 2.05047101 

6 4 7556 1758 

7 5 51248431 

8 12.2598619 

9 12.8877037 

10 235359367 

1 1 24.2568205 

12 26.4794890 

13 38.9199260 


Lumped Mm Aryololic 



11.88805 
11.88805 
23 313674 
23313674 

38.534998 


« / •yrs 4/ 

15 45.1313668 

16 57.90* 

17 57.92* 

18 80.72* 

19 80.72* 

* - Uncoupled (Reference 3) 
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Modal Frequency - Hz 


Comparison of Modal Frequencies 



Mode Number 
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3-Dimensional Beam 


Equation | 

■ii. ZSSSSSBSZSSSESal 


9 Includes Nonlinear Kinematics 


• Makes No Small Angle Approximation 


i <\ a ( 9 Hirort Jnnc) Tore i no if- Flnna^t i no 







It foliowsjhat the deflection vector of the be»m is: 


R(s) 



du x /ds’ 
dUy/ds' ds' 

1 + d / d s* 


and the direction cosine of the cross section axes is 
given by: 


0 -dty/ds dG/ds 
dT(s) _ T ( c ) dip/ds 0 -d0/ds 
ds -d0/ds d0/ds 0 


The forces and moments applied to the beam are 
related to the beam deformations by: 


dFx 

dFy 

dFz 

dMx 

dMy 

dMz 



du x /ds 

d Uy / d s 

du z /ds 

d(|)/ds 

d0/ds 

dlp/ds 


ds 


dF = F u du «■ F c de 
dM = M u du + M e de 

Where 

du x /ds 

du= dUy/ds ds 

du z /ds 


d<|>/ds" 

de = d6/ds ds 
d!|J/ds 
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The incremental force can be related to delormation of 
the beam. 


dR 

ds 


R(s+ds)-R(s) 

ds 



0" 


" du x / ds ” 

= T 

0 

+ T 

dUy/ds 


J. 


du z /ds 


r du v /ds 

dUy/ds 

riu z / as 


^ <n 
ry i U IV 

ds 


O' 

0 

1 


dF_ F(s + ds) - F(s) _ F T TdR 
ds " ds u ds 


u 


o 

0 

1 


+ F — 

Tr e ds 


Similarly, for the incremental moment.. 


dM T TdR M 
_ = Mu T __ Mu 


o 

0 

1 


. w de dF 
^ ^eds ^ ds 


The equations of motion for the beam element are 

. 2 , 


m 


d_R - T d.rSf 1 + TF 
dt2 ds Ids J 


and 


i d 2 e _ 
lo d"t2 “ 


_ + <-\m + M 

dt °dt ds Lds J 
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The equations of motion for the beam element become 


„,d 2 R _ d 

m— -n = T — 

dt^ ds 


F J T 


T dR 
ds 


oi|+, 4t 


+ T F 


2 

I d e _ _ de , de d 
0 dt2 dt °dt + ds 


+ F u( 


T dR 

ds 


Mu(T 

0 


T dR 

ds 


iOi v 

- Oj +M 
1 


d e 
e ds 


) +F < 


de 
e ds 


+ M 


Where 

e = the arc direction cosine[T] = 


de 

dt 


d<t>/dt 

dG/dt 

[dip/dtj 


<t> 

e 

V 


d T _ y de 
dt dt 

dT _ jde 
ds ds 


m = the mass per unit iength 

I 0 = the moment of inertia per unit length 


XX 

“Ixy 

-^xz 

xy 

1 y y 

-lyz 

xz 

~lyz 

*zz 


8r 



For the case in which deflections are small and the 
end of the beam is aligned with the inertial axes: 


T ~ I + 6 - 


1 -dU de 
dlj) I — d (J) 

-de d<t> 1 


dR _ T 
ds 


du x /ds 

dUy/ds 

du z /ds 


+ 


n n 

w 

0 

1 


rl ii /flc 
A 

du v /ds + 
du z /ds 



R ~ Ro + u + 

s 

r* 

e" 

-0 


0 

1 

L. -J 


The linearized equations become: 


m 


d^R _ d_ 
dt 2 ds 




+ F, 


de 

ds 


+ F4- e F 


1 

°dt 2 _ ds 


e 


m 


u'-ds : , r e ds 


+ Fu 


/dR 

Us 



ae 



For bending only, in a single plane the equations of 
motion become those for the Timoshenko beam. 


m 


dRi 


— x — d_r kG A to 

dt2 “ ds[ Ids 


-e) 


+ F, 


I 


d 2 e 

XX 


ds L 


EI XX — x 
xx ds 


]-ICA(^X-b)+ m : 


If rotary inertia effects are neglected the result is the 
Bernoulli-Euler beam. 


___ d z R x _ d z 

”772 - “ di2 


PI d R y 
E, **dP 


+ Fv 
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Kelvin- Voight Damping 


Bernoulli-Euler Beam Equation with Kelvin-Voight 
Damping 

EIu + Cu"" + mu = 0 


Allows Separation of Variables 


Theoretical Basis for Damping 


Locus of Modal Characteristics 
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Proportional Damping 


Bernoulli-Euler Beam Equation with Proportional 
Damping 

EIu' '' + Cu" + mil = 0 

Allows Separation of Variables for Pinned and 
Infinite End Conditions 


Lacks Theoretical Basis for Damping 
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Piano-Wire Damping | 



Viscous Damping Ratio 


Smaller Mass 

s = 

.0015 

Larger Mass 

§ = 

.0013 

General Mass 

£ _ 
O — 

.uO i 3 


Nonlinear Damping 


mx = -CjX - c 2 lx|x - kx 
A n+ , = A n - A n (.00 1 38)2tr - ( 0012)An 


Determined to be 
Air Damping 
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SCOLE Experiment 

Modal Characteristics 




Large 

Amplitude 


Small 

Amplitude 


Imaginary 
Fart, jui 

30 



20 


10 


.03 -.02 

Real Part, O' 
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SCOLE DAMPING 


j 


V iscous 

Damping Ratio 


Mode 

Configuration * 1 

Configuration *2 

1 

.0016 

.0013 

2 

.00 1 1 

.0009 

3 

.00058 


4 

.00 1 1 


5 

.00084 



! Nonlinear Damping is Evident for 
! Large Amplitude Motion. Analysis 
is Underway. 

i i 
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Nonlinear Damping 


Mass, Spring. Nonlinear Damper 

mx = -c|x| a |x| b x - kx 
Considering Only Light Damping ... 


ui-fk/m 
For Free Decay 

x(t)= A(t)sin(iutJ 


a - -^ixr d iii° 


c. u, u A tt+0 
m 


Solving 


ci t 


m d A 


r m 
V/ iXJ 


b^a+b 


L T l 


m 

0 c(a+h- 1 )w b A a+b 1 


A ( t) - 


m 


|_c(a+b- 1 )ujk + *o) 


1 


a+b- 1 




A(n) 


m 


c(a+b- 1 )ui b n J 


_J 

a+b- 1 


Where t + t 0 = n^IT , hi = [k/m] 


1/2 




1 


A (n) = 


m 


_c(a+b-l)k^ u i// ^2TTn. 


a+b- 1 


For Example 



c • 

a=b=0 

dA 

2iri = 

cir 

• 

X 

m 

A 

“^km 


-\i\i 

m 

a=0 

dA 

4cAui 

4cAk 1/2 

• 

b= 1 

A 

3m 

3m^/2 


-Ix|x 

m 

a= 1 

dA 

2cA 


• 

b=0 

A 

3m 
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Lop, Amplitude 


OmOMAL PACE IS 

•C-BQQK guMcr 


Log Amplitude Response 





Number of Cycles 



Lor Amplitude 




v . t s 

VHUiV ! 


Log Amplitude Response 


r 


lu^hU^’ HTvrtTi r^> 


r 

i 


Linear Only 


Linear plus 
Quadratic 

6 200 dOO 

Number of Cycles 



i 


SCOLE DAMPING 





C3 

o 

vC> 


O 

o 

ro 


o 

o 

CM 
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Number of Cycles 




Number of Parameters 


0WQ *W- page is 
6E fiPOfi QUALITY 


: <n*r?i T ; .jKjywhJs aaaaaar ^^us^mww'tmhhbwmhmm 

j^The Curse of Dimensionalit y | 


! «<f.n 


independent Modal / 
Characteristics 


800 


600 h 


40 Oh 


2 0 Oh 



1 Frequency 1 
1 Damping 
32 Deflections 

Total = 34/Mode 


Interdependent Modal Characteristics 
Distributed Parameter Model 


1 5 10 15 20 25 30 35 40 45 50 

Number of Modes 
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Concluding Remarks 


• The Accuracy of the Modal Characteristics of the 
SCOLE Configuration were Examined Using Exact 
and Approximate Solutions. 

• Sixty-Seven Modes are Required for a Static 
Deflection Error of less than 1 %. SCOLE Model 
Requires Hundreds of Modes. 

• Exact Solutions Encounter Numerical Difficulties. 

• Asymptotic Solutions in Combination with Limited 
Exact Solutions Enable Generation of a Proof 
Model with the Required Accuracy. 

• Damping Must be Incorporated into the Model 
from the Start. Proportional Damping is Not 
Adequate. 
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’ll will never be possible to have the absolute 
conviction before flight that a valid mathematical 
model has been devised for a space vehicle. 

we surely must make every effort to ensure that 

failures do not come from inadequate analysis of 
the best models available.’ 


Peter Likins 
1971 
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Direct Fixed Order Compensator Design 


KiGH-ORDbH 




SUBOPTIMAL 



IOG 




Given ... 


x = Ax + Du + 
y = Cx + w 2 , 

... design an n l c h order robust , zero set-point compensator 

x c — A c x -j- F y 
u = —Kx c 

to minimize 

J = T lim ~ Jq {x T (t)Rix(t) + u T (t)R 2 u(t))dt , 
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LQG Solution 


K = R 2 1 B T P , 

F = QC r V 2 "> , 

A c = A-BK-FC , 

P and Q positive definite solutions to 

PA + A t P + - PDR 2 l B T P = 0 

QA t + AQ +Vi~ QC t V 2 -'CQ = 0 


I 08 



0 = PA, + AfP + '£ A ?PA i - PjR- 2 ]P. + fli 

1=1 

4 - YJiAi - Q.VfSCifPiAi - Q.VjCi) + t!PBR 2 'B t Pt x 
1 = 1 

0 = A,Q + QA] + £ ArQAj - Q.VfSQl + Vi 

t = l 

+ £>,- - B t R^P,)Q(A t - B t R 2 ]P,) T + t x QC t V 2 ~'CQtI 

t = l 

0 = PA Q s + A t qs P + PJR^P.-tIPBR 2 1 B t Pt 1 
0 = ApsQ + QAls + Q^-'Q* -t±QC t V 2 -'CQtI 
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Application to SCOLE 


• Reflector Panel 


• Objective: Vibration Damping 

• 3 inputs (reaction wheels at hub) 

• 5 outputs (gyros at hub, accelerometers at reflector 
center) 

• 10 modes 


no 



MEOP procedure 


• Robustification with respect to modal frequencies 


• Robustness measure: e 


• Results: 


Order 

M 

— 2*7%)— 

Cost 

20 

0.0138 

~20 to +4 

0.229 

12 

0.0141 

-30 to +20 

0.231 

10 

0.0153 

-45 to -f 30 

0.231 

8 

0.0140 

-9 to +30 

0.235 


Optimal Projection Design 


Order 

M 

t(%) 

Cost 

20 

0.0148 

-25 to +40 

0.407 

12 

0.0156 

-50 to +50 

0.311 

10 

0.0154 

-50 to +50 

0.319 

8 

0.0154 

-40 to +40 0.322 


MEOP Designs 


in 






Open loop outputs. 



Closed loop outputs. 





TIME IN SEC 

Vibrational energy profile 
(full order designs) 


1.5 
1.0 
0.5 
0 . 

0. 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.8 2.8 3.0 

TIME IN SEC 



> I 1 t l la — sk m - ==»• — ua m -Mrf m 


Vibrational energy profile 
(reduced-order designs) 
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Torque input on x-axis at the hub. 



Torque input on y-axis at the hub. 
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Model Reference Adaptive Control (MRAC) 


Procedure: 

Find largest possible family of adaptation laws assuring 
stability, select specific adaptation law for particular ap- 
plication. 


Methods: 

Hyperstability and Positivity Concepts 


IIS’ 



Control Approach 



PI Adaptive Model Following Feedback 


lie 
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Control Objectives 


Control designed for first five modes 


2-10% damping required 



mode 

number 

Frequency 

Hertz 

Desired 

Damping 

1 

.964 

10% 

0 

A* 

.964 

10% 

3 

7. IT 

2% 

4 

7.51 

2% 

5 

9.6 

2% 


Effects of actuator dynamics not included 


1)8 
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Application to MiniMast 


ORIGINAL 

nF poor Q ° AU1 


Mini Mast Modal Frequencies 

! mode 

Frequency 

Mode 

j number 

Hertz 

Type 

i 1 

.964 

x - y 

2 

.964 

y - z 

3 

7.17 

plate 

4 

7.51 

torsion 

5 

9.6 

y - z 

6 

9.8 

X - z 

7 

10.2 

y - z 

8 

12.1 

mix 

9 

16.08 

mix 

10 

16.8 

mix 


Mini Mast Model Actuators 

number 

FEM 

Actuator 

Coordinate Frame 

Force 


Point 

Type 

X 

Y 

z 

Limitations 

1 

334 

Linear Actuator 

1 

0 

0 

30 newtons 

2 

336 

Linear Actuator 

1 

0 

0 

30 newtons 

3 

335| 

Linear Actuator 

161 


6 

30 newtons 

4 

337 

Linear Actuator 



0 

30 newtons 

5 

338 

React ion Wheel 

0 

0 

1 Rot 

50 Ft. - lb. 
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e.s 

0. 

-.6 

0. e.s l.e 1.S 8.0 8.5 3.0 3.5 4.0 4.5 5.0 


Mode 2 Reference Error 



Seconds 


II 9 













Seconds 


Input 3 Magnitude 



.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5. 

Seconds 


Input 4 Magnitude 



- 20 . 


-40. I 1 1 1 1 1 1 1 » l I 

0. 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Seconds 
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Mode 6 Position 




.5 2.0 2.5 3 . 

Seconds 

Mode 6 Velocity 


Seoonds 


Mode 9 Position 


Seconds 

Mode y Velocity 


2.0 2.5 3 . 

Seoonds 
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Outlook 


• Effects of Actuator Dynamics 

• Refinement, of STAC 

• System Identification 

• Experimentation 
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Abstract 


large space structural systems, due to their inherent flexibility 
and lew mass to area ratio, are represented by large dimensional mathe- 
matical models. For implementation of the control laws for such systems 
a finite amount of time is required to evaluate the control signals; and 
this time delay may cause instability in the closed loop control system 
that was previously designed without talcing the input delay into con- 
sideration. The stability analysis of a simple harmonic oscillator 
representing the equation of a single mode as a function of delay time 
is analyzed analytically and verified numerically. The effect of in- 
herent damping cn the delay is also analyzed. The control problem with 
delayed input is also formulated in the discrete time domain. 

I. Introduction 


Large flexible space structures have been proposed for possible use 
in camnunicaticns , electronic orbital based mail systems, and solar 
energy collection. The size and the low mass to area ratio of such 
systems warrant the consideration of the flexibility as the main contri- 
bution to the dynamics and control problem as compared to the inherently 
rigid nature of earlier spacecraft systems. For such large flexible sys- 
tems, both orientation and surface shape control may often be required. 


The equations of motion describing the shape of any large space 
structure are either represented by a few partial differential equations 
or a large number of ordinary differential equations. As the partial 
differential equations are difficult to solve for control system design 
purposes, the structural dynamics are commonly described using Finite 
Element Methods (FEM) . Two typical large space structures narrely the 
Hoop/ Column antenna^ and the Space Station initial operational configura- 
tion (IOC) 4 are both described using 672 degrees of freedom. Thus the 
dynamics of a large space structure can be written as^: 


M Z + K Z = U c 


( 1 ) 
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where 


M = NXN mass/inertia synmetric matrix 
K = NXN stiffness symmetric matrix 

Z = NXl generalize coordinates representing the degrees 
of freedom 


Uc= 


influence of the external forces in each degree of 
freedom = B’U. 


With the modal transformation 
Z = <t> q 

and the properties of the modal transformation such as 


= I 

iTvi _ 


2 2 


2-i 

... « n ] 


()> K<{> — diag [a>-| co ^ 
and neglecting the higher modes, equation (1) can be written in standard 
state space form as 

X = AX + BU 

where 

X = 2nxl state vector representing modal coordinates 
and their velocities [q,q] T 


( 2 ) 


U = mxl 


A = 


0 


control vector 


nxn 


-to. 


i 0 


21 


— 0 ) 


“n| 


system matrix 


B = 


0 


nxm 


(J) T B' 


nxm 


control influence matrix 


II. Control with Delayed Input 


The proposed control systems for large space structures are based 
an state variable feedback of the form: 


U = -FX 


(3) 


ns 



STABILITY ANALYSIS WITH DELAYED INPUT 


and the control gain matrix, F, is designed using techniques such as the 
linear quadratic regulator (LQR) theory^, pole placement', and/or linear 
quadratic Gaussian/loop transfer recovery (LQG/LTR).® 

For the case when the^complete state is not available for feedback, 
an estimate of the state, X, is obtained using an appropriate estimator 
from the measurements of the form 

Y = CX (4) 


where 

Y = £xl measurement vector 

C = £xn sensor influence matrix 

In general, it is assumed that the estimated* state, X, is instanta- 
neously available. As the state estimator is implemented using a 
digital computer and the number of the status (2n) is of the order of 
hundreds for a large space structure, the computational time becomes 
appreciable. Thus, in the present paper, the stability of the closed 
loop control system, with the control as given in equation (2), is 
analyzed as a function of the delay time (h) using the modified control 
law of the form: 

U ( t ) = -FX(t-h) (5) 

The characteristic equation of the closed loop system 


X = AX ( t) -BFX(t-h) (6) 

is given by 

— cih 

G(s,h) = det (sI-A+BFe ) = 0 (7) 

which, in turn, can be written as 

G(s,h) = P.(s)e“ shi = 0. (8) 

i=0 


The roots of the characteristic equation, (8), as a function of the delay, 
h, are obtained from the corresponding auxiliary equation^ 


G'(s,h) = I P. (s) (1-Ts) 21 (1+Ts) 4n 21 = 0 
i=0 


where 


-sh 


= [ 


1-sT -I 2 
1+sT J . 


(9) 

( 10 ) 
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The value of T for which the roots of the equation (9) cross the 
imaginary axis in the s-plane is obtained and the corresponding h is 
evaluated using the relation, (10) . 

III. Example of a Harmonic Oscillator 


The equation of motion representing the ith structural mode is 
the familiar harmonic oscillator and is given by 


x. + oj. x. = f. 
1 li l 


( 11 ) 


Considering the delayed velocity feedback of the form 


= -2^0)^ (t-h) (12) 

with 

“i = 6' = °* 5, 

the characteristic equation is given by 

G(s,h) = s 2 +36+6se -sh = 0 (13) 

= Jq P i (s)e“ shi = 0 

where 

P (s) = s 2 + 36 
P 1 (s) = 6s 

The corresponding auxilary equation is given by 


I P. (s) (1-Ts) 2l (1+TS) 2—21 =o (14) 

i=0 

(s 2 +36) (1+Ts) 2 +6s (1-Ts) 2 = 0 

T 2 s 4 + (2T + 6T 2 ) s 3 + {I+IQV 1 -!??) s 2 

+ (72T+6) s+36 =0 (15) 

Using the Routh-Hurwitz criterion, it can be found that the roots of 
equation (15) cross the imaginary axis at oj = 9.7 for T = 0.0426. 

The corresponding delay (h) can be calculated from the relation (10) 
with s = jo) and is 0.16. This result can also be verified directly 
for this sinple system with the substitution s=j(u into equation (13) , 

resulting in the value of oj and delay h for which the roots of the 
characteristic equation cross the imaginary axis. 
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Thus, equation (13) can be written as (keeping ^ and u^) : 
(uk 2 -^ 2 ) + j ( 2 ^ i (i) i o))e = 0 


or 


2 2 

(u) . -o) +25 .oj.o> sinwh) 4* j25.o).o) cosojh = 0 . 

i li J 1 i 

For equation (17) to be satisfied 

cos u)h = 0 or oih = ^/o 


and 


or 


2 2 

0). ~o) + 25 .oj.o) = 0 

i i l 


o) = 5.0). + o) . / 1+5 . 2 
ii—i l 


Taking the positive value for o), the delay h, is given by 


it/. 


h = 


oj. [5 .+ /ThT" ] 

11 1 J 


(16) 

(17) 


(18) 

(19) 


( 20 ) 


The value of h for 5 . =0.5 and o), = 6 is 0.16 and thus the earlier 

i i 

result is verified. It is observed that an increase in damping reduces 
the tolerable delay (h) in the input. 

The equation of motion of a single mode with inherent (natural) 
damping and velocity feedback can be written as: 

X+ 2 c[cu i X-Ho 1 2 X = f = - 2 ? i u i X(t-h) ( 21 ) 

where 5 ^ is the inherent damping ratio. 


The corresponding characteristic equation is given by 



s 2 +2c !uj . s-Hjj . 2 +2£ .(jj.se = 0 . 

11 1 l i 

( 22 ) 

After substituting s = jo>, equation (22) can be written as: 



2 2 t 

(o) . -o) +25 .03 . 0 ) simoh) + j (25 .oj .o )+25 .oj .ojcosoih) = 0 
1 ii J i i l l 

(23) 

For 

equation (23) to be satisfied for all o) and h, we have 



25 . 0 ) .+ 25 . 0 ) . cosojh = 0 
s i 1 11 

(24) 

or 




cos o)h = -5 . / 
15 

(25) 


i 
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Thus, for coswh = < 1, the inherent damping must be less than damping due 
to control for instability. For the system will always be stable. 

With the value of u)h from equation (25) the frequency w can be 
calculated as : 

co — ( jj , [ / £ £ 1 2 + / 1+^2 -£ 1 2 - ] ( 26 ) 

1 i i i i 

and selecting the positive value of ui, h is given by: 

cos 

h = : 


a ),[✓ ± / (27) 

1 i i i i 

For £ . = C * it can be seen that the delay, h, is half the undamped 
natural period of vibration. As the damping due to control increases, 
the tolerable delay (h) decreases and is in accordance with the observa- 
tion made in the case without the inherent damping. The effect of inher- 
ent damping in the system is to increase the amount of delay that the 
system can tolerate without become unstable as compared to the case 
without damping 


IV . Discrete Time Domain 

As the controller is implemented on a digital computer, it may be 
more natural to consider the delayed input problem in the discrete time 
domain . 

The equations of motion as given by equation (2) can be written in 
the discrete time domain as 


X( i+1) - A d X(i)+B d U(i) 


where 


A = e AA , B = / e A(t A) B dt 

a a 

o 

A = discretization time. 


(28) 


The delayed input . problem can be considered in discrete time in one of 
the two following ways: 

i) Designing the controller of the form U = -FX(i) without taking 
into consideration the delay and then examining the effect 
of delay on the stability of the closed-loop control system. 

The control gain matrix F is designed such that the matrix (A^-B^F) 
has the eigenvalues within the unit circle. Then the delay is introduced 
into the control law as: 
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U(i) = - FX(i-£) 


(29) 


and 


X(i+1) = A d X(i)-B d FX(i-£) . 


(30) 


The stability of equation (30) can be studied using the augmented system 
given by 


X(i+1) 


A d 

0 

0 

0 

- B d F 



X(i) 

X(i) 

= 

I 

0 

0 

0 

0 



X(i-l) 

X(i-£+l) 


0 

1 

0 

0 

I 

0 



X(i-£)j 











or 


Z(i+1) 


Z (i+1) = A d Z(i) 


'x, 

A, 


Z(i) 


(31) 


(25) 


(ii) Designing the control by taking into account the delay in 
the input. 6 ’ 11 

Equation (28) can be modified as : 

X(i+1) = A d X(i) + B d U(i-£) (32) 

The control law of the form U(i) = -FZ(i) can be designed from the 
augmented system: 


X(i+1) 

U(i) 

U(i-l) 

U(i-£+l) 


A, 0 0 

d 


0 0 

I 0 

0 0 


0 B 

c 

0 0 

0 0 

1 0 


X(i) 

U(i-l) 


U(i-£) 

'Xj 

Z(i) 


0 

I 

0 


U(i) 


(33) 


or 


Z(i+1) = A d Z(i) + B d U(i) 
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Thus the input U(i-£) is a function of the previous inputs, U(i-£-l), 
U(i-£-2) , . . . , and the previous states X(i-£). Though this design can 
take delay into consideration, the sequence of the control signals: 
U(i-£), U(i-£+l),... must be generated at an interval of one step and, 
thus, the original delay problem is not completely solved. 


Conclusions 


The effect of delay in the input on the stability of the continuous 
time controller that is designed without taking this delay into consider- 
ation is presented. The closed-loop control system of a second order 
plant becomes unstable for a delay of 0.16 seconds, which is only 16 per- 
cent of its natural period of motion. It is also observed that even a 
small amount of inherent (natural) damping in the system can increase the 
amount of delay that can be tolerated without the system becoming unsta- 
ble. The delay problem is formulated in the discrete time domain and an 
analysis procedure is suggested. 
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Abstract 

The minimum time attitude slewing motion of a 
rigid spacecraft with its controls provided by 
torques and forces, which have their upper and low- 
er limits prescribed, is considered. The two-point 
boundary-value problem is derived by applying the 
Pontryagin's Maximum Principle to the system and 
solved by using a quas i 1 i nearization algorithm. 

The nominal solutions to the problem as well as 
the starting values of the total slewing time and 
the unknown initial costates for this algorithm 
are generated by using Euler's eigenaxis rotation 
theorem. It is pointed out that one of the four 
initial costates associated with the quaternions 
can be arbitrarily selected without affecting the 
optimal controls and, thus, simplifying the compu- 
tation. The minimum slewing time is determined by 
shortening the total slewing time until at least 
one of the controls becomes a bang-bang type. 
Several numerical tests for the rigidlzed SCOLE 
model are presented to show the applications of 
the methods. 

1. Introduction 

The problems of large-angle attitude maneuvers 
of a spacecraft have gained much consideration in 
recent years . 1-8 In these papers, the con- 

figurations of the spacecraft considered are: 

( 1 ) completely rigid, ( 2 ) a combination of rigid 
and flexible parts, or (3) gyrostat-type systems. 
The performance indices usually include minimum 
torque integration, power criterion, and frequency- 
shaped cost functionals, etc. Also some of these 
papers used feedback control techniques. In this 
paper, we try to concentrate on the minimum time 
slewing problem of a rigid spacecraft. 

In Ref. 2 , the author studied the rapid 
torque-limited slewing of SCOLE about a single 
axis (x-axis) about which the spacecraft has a 
small moment of inertia. The control torque about 
this axis is of a bang-bang type or a bang-pause- 
bang type. The author computed the slewing motion 
on the simplified model of the rigidlzed SCOLE^ , 
then worked on the practical rigidized model (with 
nonzero products of inertia); hence, this leads 
to a large error of the attitude after the slewing. 
Also it seems that no details were given for the 
controls about the other two axes (y, z). 

In the present paper, we apply optimal control 
theory (Maximum Principle) to the slewing motion 
of a general rigid spacecraft (include the rigid- 
ized SCOLE, without simplification). The slewing 
motion need not be restricted to a single-axis 
slewing. The attitude error at the end of the 
slewing can be made as small as required. 

*Research supported by NASA Grant NSG-1414 
^Graduate Research Assistant, Dept, of Mechanical 
Engineering, Student Member AIAA 
**Professor of Aerospace Engineering; Associate 
Fellow AIAA 


All the controls (torques and forces) are computed 
and the minimum slewing time is found by using the 
quasilinearization algorithm for the resulting two- 
point-boundary-value problem. 

2. Attitude Description and State Equations 
2.1 Attitude Description and Euler Rotation 

Let a - [a^ a^ represent a set of unit, 

orthogonal vectors of an Inertial reference system, 
and b = [Fj F~] a set of unit, orthogonal vec- 
tors of a Dody-fi/.ed coordinate system of a space- 
craft. Then, the attitude of the spacecraft rela- 
tive to a can be described by a direction cosine 
matrix C with C satisfying the relation 

F ■ C a ( 1 ) 

and 

q 0 +q l' q 2" q 3 2 ( q iVV3) 

C= 2(q l q 2* q 0 q 3 ) q S +q i- q T q f 2(q 2 q 3 +q 0 q,) 
Z(q l q 3 +q 0 q 2 ) 2(q 2 q 3' q 0 q l) q 0 +q 3' q l' q 2 

( 2 ) 

where q = [q Q q. q^ q.]^ is the attitude quaternion 
vector and subject to a constraint equation 

q T q = 1 (3) 

It can be seen that q can be used not only to 
represent an attitude orientation of a spacecraft, 
but also to describe a rotation of a rigid body 
(spacecraft). For example, when a rigid spacecraft 
rotates about an axis defined by a unit vector 
€ = ^ e l c 2 ^ xe< ^ in both 1 and *5, the quaternion 

describing this rotation is 

q 0 = cos (e/2), q. = t- sin (e/2) i = 1,2,3 (4) 
where e is the rotation angle. 

The Euler rotation theorem tells us that an 
arbitrary orientation of a rigid body can be accom- 
plished by rotating it about a certain eigenaxis, 
c = [e, e£ £ 3 ] » through 9 angle from its initial 
positiAn. By means of this theorem we can find 
the desired rotation quaternion, q, between the 
initial position q( 0 ) and the final orientation 
q(tf) by the relation 


q o ■ 


q 00 q 10 q 20 q 30 " 


q 0f~ 

q l 


“ q 10 q 00 q 30 ~ q 20 


q H 

q 2 


” q 20 ” q 30 q 00 q 10 


q 2f 

q 3 


q 30 q 20 ' q 10 q 00 


q 3f 


'(5) 

where the second subscript " 0 " and "f” represent 
the initial time and final time, respectively. ' 
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The associated c and 6 can be obtained by the 
following relations 

e * 2 arc cos q Q . e, * l-q^, 1 * 1,2,3 

(6) 

2.2 Kinematical and Dynamical Equations 

The attitude quaternion and the angular velo- 
city of a rigid spacecraft satisfy the following 
kinematical and Euler dynamical equations. 


q = (1/2) u q 

(7) 

I in * iii I w + B u 

(8) 


where 

<*> _ angular velocity vector in the body system, 

U> * 

u _ control torque and force vector, 
u • [u lU2 u 3 . . . u n ]T 

and 


0 • 

“(D-j 

-“2 

•“3 


0 

“3 

*“2 

“1 

0 

“3 

'“2 

hi * 

*“3 

0 

“1 

“2 ' 

’“3 

0 

“1 

» 

*2 

'“I 

0 

.“3 

“2 

-“1 

0 





'hi 

- 

h 2 * 

r !3 





_I 12 


\ n - 

! 23 





**13 

- 

! 23 

*33 






and B is a 3xn alignment (control influence) 
matrix. Eq. (8) can be rewritten, by pre- 
multiplying by the inverse of I, as 

i - I- 1 u i u + r 1 b u (9) 

The associated initial and terminal boundary 
conditions of the states, q, w, are prescribed: 

q(t*0) , u>(t*0) ; q(t f ), ai(t f ) (10) 


For this case the minimum time rotation of the 
spacecraft about one of its principal axes will 
yield the following results: the control torque 

about this axis Is of a bang-bang type, while the 
other two torques remain zero. Otherwise, if the 
slewing motion is not about a principal axis, none 
of the 3 controls remain zero, but we can reason 
that at least one of the 3 control inputs reaches 
Its bounds, except some jumps at the switching 
points during the period, tf. As for a general 
case where the control torques are about a body 
axes system which does not coincide with the prin- 
cipal axes (non-diagonal I) and some additional 
control forces, u^ u$, ..., u n , are available, 
the control laws become more complicated. 

To handle the problem in which some controls 
reach their bounds and others do not, we Introduce 
an additional cost function 

J * x / tf u t Ru dt (13) 

2 0 

where u Is the control vector, R Is a proper 
weighting matrix. From Refs. 3 and 8 , we 
can see that, for the case of rest-to-rest (l.e. 
<i>(0) ■ 0, di(tf) * 0) slewing with only 3 control 
inputs involved. If we use only (13) as a criterion 
and tf Is long enough, the control torques are 
approximately linear functions of time, and the con- 
trols will not reach their saturation levels. But 
if we shorten tf in order to find a minimum time, 
some of the controls must reach their bounds and, 
thus, contribute more effort to the slewing. By 
continuing the^shortening of tf, we can get a parti- 
cular value, tf, during which at least one of the 
controls remains as bang-bang with one switching 
point, while others are generally not of the bang- 
bang type. This value, t* is called the minimum 
time which is required. 

The motivat iors for using (13) as our cost 
function are: 

1) Ease of using the quasilinearization 
algorithm 

2) No need to determine the switching points. 


3. Optimal Control - Two Point 
Boundary Value Problem 


3) Easy to guess the unknown initial values 
of the costates. 


In this paper, we try to minimize the slewing 
time tf, 

t f = / tf dt (11) 

0 

under the constraints that the elements of the 
control vector u have their upper and lower limits, 
respectively 


jmin 


u. < u, , 
j - jmax* 


1,2,3, 


n (12) 


Generally speaking, minimization of tf under 
the constraints (12) will result in a so-called 
two-point boundary-val ue problem in which several 
controls (at least one) will reach their bounds 
during the slewing time, tf. To explain this 
point, let us first consider a well-known special 
case where there are only 3 control torques, u^ , 
U£» u 3 , about the 3 principal axes of the space- 
craft, respectively (i.e. diagonal matrix I). 


3.1 Necessary Conditions 

The Hamiltonian, H, for the system (7), (9) 
and (13) is 

H = (l/2)u^Ru + p^q + r^iL 

* (l/2)u"^Ru + {l/2)p^u>q + r^( 1 ulw+l B u) 

(14) 

where p and r are costate vectors associated with 
q and w, 

P = [p 0 p 2 P 3 ] T . r = [r, r 2 r 3 ] T . 

By means of the maximum principle, the nec- 
essary conditions for minimizing J, are 

p = - (3H/3q) , ==> p = (1/2) w p (15) 

r = - (aH/ao) , ==> r = g(u,r ) + (l/2)[q]p 

(16) 
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where g( w ,r) is a 3x1 vector function of and r, 

and the detailed form of gfui.r) can be found In 
Appendix I; [q] Is a 3x4 matrix 

' q l * q 0 - q 3 q 2 ' 

[q] - q 2 d 3 ■% 

q 3 -q 2 q, *q 0 

The initial values of p, r are unknown, p(t=0), 
r(t-O). 


± ? _ JjX! •• rli cs of t hn T i a 1 JVI ue: of 

The key to settle this problem is to find the 
unknown initial values of the costates 

p( 0) ■ [ p 00 p 10 p 20 P 30 ] T and 

r(0) * [r 1Q r ?0 r 3Q ] T 

Notice that the coefficient matrix of the right side 
of Eq. (15) is anti -symmetric , so. 


If u is a 3x1 control torque vector and B is 
a 3x3 nonsingular matrix, R can be a positive- 
definite matrix defined by 

R = B T B O 7 ) 

From 

=0, *■*'- R u + B t I" 1 r = 0 

du 

or u * - R ’ 1 B t I " 1 r 

= - B* 1 T 1 r (18) 

we have 9 


u jmin’ 

1f u j < u jmin ; 

-(B* 1 I' 1 r) j 

■ 1f u jmin - u j - u Jmax ; 

u jmax • 

1f u j > u jmax. 


j = 1,2,3. (19) 


If u is an nxl (n>3) vector, B is a 3xn matrix, 
the R formed by (17) is a semi-positive-definite 
matrix. To circumvent the singularity of R, we 
introduce a 3>1 vector, v, 

v * B u (20) 

Thef1 1 /olX 
( 3H/3v ) * 0, ==> v * - 1" 1 r (21) 


p T p * 0 i.e. p T p = constant 

The extra constant is usually treated as an unknown 
and is determined by iteration. This results in 
more computational effort. However, as we shall 
prove, this unknown constant can be arbitrarily g 
selected without changes in the optimal controls. 

Compare Eqs. (7) and (15), they have the same 
coefficient matrix on the right sides. Therefore, 
they have the same state transition matrix. Let 
Q represent this 4x4 matrix, then the q and p at 
any instant can be obtained by 

q 1 Q q(0), p = Qp(0) (24) 

We know that Q satisfies the following matrix dif- 
ferential equation 

Q * (1/2) «Q (25) 

Ref. 10 shows that Q, the solution of (25), has 
the form 

q ll - q 12 - q 13 - q 14 

q 1 2 q ll q 14 - q 13 

q 1 3 ' q 14 q ll q 12 

q 1 4 q 13 * q 12 q 11 




By means of pseudo-inverse of matrix B, B + , we can 
get u , _ T i 

u = B + v = B t (8B t ) _i v 

= - B T (BB T )- 1 (I- 1 r) (22) 


The control laws are 


I 

[ u jmin* 

if 

Uj < Uj min; 

U j " | 

' -(B + I' 1 r) j , 

if 

u , J < u • < u • 
jmin - j — jmax 

1 

u jmax* 

if 

u j ” u jmax. 



j * 

1.2 n. (23) 


Note that Eq. (23) is reduced to Eq. (19) if B" 
exists. 

In summary, we seek the function q(t), u>(t), 
u(t), p(t), and r ( t ) which satisfy the equations 
(7), (9), (15-16), (23) subject to the boundary 
conditions (10). 


On substituting Eq. (26) into Eq. (25) we can 
verify that only 4 of the 16 qy are independent. 

We rewrite and rearrange the first equation of (24) 
as 



where q Q0 = q Q (0), Q 10 s qi(0), etc. It is clear 
that the matrix, Qq, in Eq. (27) is orthogonal, so 

q = QT q (28) 

From Eq. (28) we get 

q ll f q U + q 13 * **14 = ' 
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This means that Q Is also orthogonal. On the other 4< initial Values of Costates 

hand, we have a similar equation for p, and tKe~SVewing Time 


p 0 a 


(29) 


where Pq has the same form as Qq In Eq. (27). 
After substituting Eq (28) into Eq. (29) and eli- 
minating q, one arrives at 


(30) 


A 1 


d 0 

- d i 

* d 2 

- d 3 ' 


Po 1 

p 1 


d l 

d o 

” d 3 

d 2 


q l 

p 2 


d 2 

d 3 

d 0 

- d 1 


q 2 

A 


_ d 3 

” d 2 

d 1 

d 0 . 


q 3. 

where the constants 

d 0* d 1 

» ^2* 

d 3 

, are 

M 

[ 

O 

O 

or 

q 10 

q 20 

q 30 

1 

[ p 00 


' q 10 

~ q 20 

" q 30 


q 00 

q 30 

" q 20 


q 30 
q 00 
q l 0 


q 20 


-q 


'10 

q oo 


p 10 

p 20 

p 30 


(31) 


Eq. (30) represents the relationship between the 
quaternion and the associated costates. Eq. (30) 
can be rewritten as 


o 

CL 


p l 

3 

p 2 


P 3 

L 


^0 


-q 


’ p 2 

q l 

<ln 



CL 

o 

J 


d 1 


d 2 

. 

1 

CL 


(32) 


By means of Euler's eigenaxls rotation theorem, 
from the known attitudes at the initial and final 
time, q( 0 ) and q(tr), we can find a unit vector 
(eigenaxls), e, which is fixed in both the body 
axes and Inertial coordinate system, and a rotation 
angle, 0 *. Then the attitude changes from q(0) to 
q(tf) can be realized by rotating the spacecraft 
about the axis, e, through the angle, e . 

Theoretically, there are many ways through 
which we can change the attitude from its initial 
value, q(0), to its final value, q(tf). For exam- 
ple, this change of attitude can be achieved by 
successively rotating the spacecraft about the x,y, 
z axes (l.e. 1 - 2-3 rotations) through certain dis- 
placements In the angles, 0 ,, 02 * 63 , respectively. 
To do this way, we need to Speed upland slow down) 
the spacecraft 3 times, and the total rotation 
angle Is, e,+ 0 ~+ 0 3 . On the other hand, for the 
Euler rotation /we only need to rotate the space- 
craft about e once through the angle 0* which Is 
less than the total angle required by any other way. 
Since the Euler rotation is simple and requires a 
smaller angle, it may take less time and consume 
less energy (torques and forces). Therefore, in 
view of our cost functions, (11) and (13), it is 
reasonable to think that the optimal slewing is 
near the Euler rotation. We shall call this rota- 
tion the ‘'expected rotation," which is determined 
only from the initial and final attitude of the 
spacecraft and will be used in obtaining a set of 
approximate unknown initial values of the costates 
and the starting solution of the quasi 1 inearization 
algorithm. 


Substituting of Eq. (32) into Eq. (16) results in 


4.1 Initial Values of Costates 


r * gfu.r) -(1/2) C d (33) 


where d * [d, d. d 3 3 T , C is just the attitude matrix 
given by Eq. (2). J It can be seen that r is inde- 
pendent of d Q , from Eq. (33), and u depends only 
on r, from Eq. (23). Therefore, u is also inde- 
pendent of d Q . This means the arbitrary selection 
of the valueof dg yields the same extremum con- 
trol, u . Now we can explain the results in Ref. 

11 . In view of Eq. (31), we have 

d 0 +d l +d 2 +<J 3 S P 00 +P 10 +P 2Q +P 30 

If we set d n = 0 the norm of the initial costates 
in Eq. (34)reaches its minimum, the solution of 
which is considered in Ref. [3]. From Eq. (31) we 
can also know that dg r 0 means 

p( 0 ) T q ( 0 ) = 0 

which is the orthogonality condition obtained in 
Ref. 11 


Before starting the quasi 1 i near ization algo- 
rithm, we need to guess the unknown initial values 
of the costates, p and r. Considering the analyt- 
ical solution about a single principal axis maneuver 
in Ref. 3 we define a rotation angle e(t), about 
an arbitrary axis e . 

e(t) = e(0)+e(0)t + ] j e(0)t 2 + ^ e(o)t 3 (35) 

where e(0), e(0), e(0), and e(0) are constants to be 
determi ned . 

For simplicity, here we only consider the solu- 
tion of 0 with the following boundary conditions 

0 ( 0 ) =■ 0, 0 ( 0 ) i 0, e(t f ) = e* e(t f ) = 0 (36) 

These conditions correspond to the boundary condi- 
tions of the states 

q(0), 40 ) f 0, q(t f ), 4t f ) = 0 
Substituting Eq. (36) into Eq. (35) yields 

0(0) = (6e*/t^) - (4e(0)/t f ) (37a) 

0(0) = -( 1 2e*/t|) ♦ (6e(0)/tf) (37b) 


(37 



Itv ,nH r ^c%o U ] er I? rotatfon - the angular veloc- 
ty and its derivatives are expressed as follows 

W 3 C0, (I) *e6, <I) * ee (33J 

To approximately determine the initial values 

° f /n^ r, „ We need t0 use the rfynamical Eqs. ( 9 ) 
and ( 3 3) Upon using Eqs. (20) and (21), substi- 
tuting v (u) into Eq . (9) and solving for r, we 

r - I « I» - l2 i {39) 

and the derivatives 

r = Jf (I U I U ) - I* u (40) 

that h C TS fc-f ,n,e * fr ° m £q ‘ (33) ’ and not ’ n 9 

d = 2 C T [g(u, r) - r ] ( 41 ) 

Foc^nq’a?! 1 " 0, by puttin 9 E( 1S- (37-38) into 

r?oi inHH 1 ‘ac e r Can the a PP roxim ate values of 
r(0) and d. As for p(0), we can set p n (0) » 0 

d ^nH d °fn? n ^i a o b ^ trar il y chosen > a ^ d for 

do and Pl (0), 1<=1,2,3, by using Eq. (32). 

Now we determine e(0) from known initial value 
^ nera,1 >* “(° Is not equal to e9(0) since 

e = E e(0) -u(0) 

To find a minimum value of e T e, we differentiate 
e E e with respect to e{0) and note that Jr . ? 
we get L L 1 » 


a™ 9 the initial values p(0) and r(0) 
obtained above, and integrating the differential 
quations (7) (9), (15-16), with Bu in Eq. ( 9 ) 

replaced by v in Eqs. (20) and (21), we can at it a 

uhi ^iJ U k S ’ q P(t), and r(t), 0< t <tf 

which will be used as the starting values of~the f ' 
quasilinearization method. * T lne 

4,2 Initial Value of t f 

closed t S ^ r I in? Val » e t f ( ° ) *«««<•* to be made as 
the minimum time, t f , as possible. This 
can be done by using the techniques similar to 

an r e i S » r ’ be t ! b ° Ve ‘ Suppose the slewing motion 
Lu Lr? tat TK n about a vector, E , through 
1 angle, eft). Then, by putting the first two 
nations of Eq. (38) and Eq. (20) into Eq. (8), 

1 £ 9 = E ' 1 E + v (43) 


For simplicity, we only consider the case 
♦ i£ 1n kT ’ v imax* Then, let c^ * v im , v and v=c-t - 
the above vector equation can be eip^ssed as thj” 
following 3 similar equations for e(t): 


9 = b< 


+ c,- Ti 


i = 1,2,3 


?!! ere i a, ' J b1, and c i are constants, Tj is the 
normalized control about the 1 th bod^axis and 

hill 1 1* 1.2.3 (45) 

Uni Eo^nfi? ° f ^ i 44) Wlth the boundary condi- 

:? fclffjsa 

see that the control for this problem is of a 
bang-bang type and the the problem can be solved 
2 l ° 9Ct th $ minimum time t* ff (i«l 

Since the only minimum time, t/, that everv 
»' '«• («> c. accept it it,’ 

miJui ““ *’* «" •* 

el’t™ l:“f, "> 

equations and it converges quadratical ly. 

In the quasilinearization algorithm the lfn 
? an « d sta ^ a nd costate equations are solved us - 

C^n.^ r meth ? d 0 , f P artic ular solutions. 1 ? The 
computational values of u, which satisfy Eq ( 23 ) 
re ieter.,„, d by , t , ch , que g 

thp fnii 9 '- TbC n " nlmum slewing time is obtained by 

ewi f J t°ime n9 t Pr ° Cedure - F ? r a " a ^*ed given * 

the I V ’ t f’ as 8 result of the iterations 
the arr1ves at the (converged) values for 

L whe h Ml costates - Th en. a check is made as 
to whether one of the control inputs is of a bana- 

the* 9 ■ y ? € ' yes » thls slewing time is designated 
the minimum time. If not. the assumed t f should 

oH VTe^ thC Ueration c > cle restarted in 
° r tl t L ° deter "" ne new values for the initial co- 
states and new time histories for the control effort 

S. Numerical Results 

Finally, we apply these methods described in 

E a P \7',r TSTS t0 the SC0L£ slewing motion. 1 
Fig la shows the SCOLE configuration It is com 

SSJ" » S r* , Sh " U ! t *" d * reflecting^*"" 

by a Jlexihlf UH' tenna c ! s attached to the Shuttle 
y a flexible beam. Since we only consider the 

bilUv of the ^ 9id . SC0LE in E1 >is paper, the flexi- 
Dinty of the beam is ignored. The X, y Z axes 

Ditch L SPaCCCraft axeS corre sponding to roll, 

P tch and yaw axes, respectively. The controls 

abou^theVV^ Paper include three aments 
about the X, Y, Z axes of the system and two forces 

applied at the center of the reflector in the^ y 

OLE and n \he n,y \ ^ fnertia «ters of the* 

SCOLE and the saturation levels of the controls 

! 11 = 1 .132,508 , I^^ 1 7,007,447, l 3J - 7,113,962 ^ 

I 12 = -7,555, I 23 = 115,202, i^. 52,293 ' 

T ° rqueS: l u i Lax = 10 .°00 ft-lb. i = 1,2,3. ^ 

Reflector actuators: | fj - 800 lb. . j. i 
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We have done some numerical simulations for 
the following cases: (a) A diagonal Inertia 

matrix I Is used. The control Is assumed to be 
provided only by torquers on the Shuttle. No con- 
trol forces on the reflector are assumed. The ex- 
pected rotation Is a rotation about one of the 
three principal axes, through 20 deg., from rest 
to rest. The result Is exactly the same as that 
of the theoretical analysis discussed earlier In 
this paper, i.e,, the control torque about the 
slewing axis is of a bang-bang type while the 
others remain zero. 

(b) Extend the Inertia matrix In case (a) to 
a non-diagonal form. The expected rotation Is a 
rotation about one of the three spacecraft axes 
and the rotation angle Is 20 deg. 

Figs. 2-3 give the control torques and atti- 
tude angles (1-2-3 Euler angles) for the expected 
rotations H X-axis slewing" and "Z-axIs slewing,* 1 
respectively. Fig. 2a shows that u Is nearly ofa 
bang-bang type, while u y » u 2 are not. The non-zero 
contributions of the u/and u 2 are due to the off- 
set of the inertia distribution of the SCOLE con- 
figuration (non-diagonal matrix I). Similar situa- 
tions are shown in Fig. 3a, where u 2 is nearly of a 
bang-bang type and the others are not. 

The starting value of ti°) for these slewings 
(X- and 2- axes) are U u ^= 12.5749 sec. and 
tl°)= 31.5166 sec., respectively, by using the 
method in section 4.2. ^The minimum time, t*, we 
actually obtained are t f = 12.57 sec. and tt = 

31.33 sec., respectively. These results Indicate 
that the estimated values for tf are very accurate. 

In Fig. 2b, e x changes from zero to 20 deg., 
but 6y and e» change very little during the slew- 
ing and finally approach zero. The non-zero changes 
in 6y and e 2 are also due to the offset of the 
structural distribution of the SCOLE. 

In Fig. 3b, unlike the case in Fig. 2b, the 
e x changes greatly. This change is due to the 
differences in the moments of inertia about the 
X-axis and Y-axis. 

Fig. 4 shows the control torques for the 
"X-axis slewing" with a slewing time t f = 15.37 sec., 
which is 2.8 sec. more than the minimum time ti = 
12.57 sec. (Fig. 2). The controls are almost lin- 
ear functions of time (rest-to-rest slewing), u 
is less than the saturation level, and u y , u z are 
ne :ero. From Fig. 4 and Fig. 2a, we see that 
mu * more control effort (approximate 50%) is 
ed If we increase the slewing time a little, 
ther feature of using a longer slewing time in 
c computation is that it needs less number (4 
^es) of iterations for convergence than by using 
minimum slewing time tf (12 times). These pro- 
erties suggest that, in practical applications of 
his problem, it is not necessary to seek the mini- 
jm time, tf, and the associated extremum controls. 

It is enough to know approximate values of the tj 
and the controls. The results of Fig. 5 for the 
”Z-axi$ slewing" are similar to those of Fig. 4. 

(c) Following the case (b), we now add two 
control force actuators on the reflector, f y and 

v 


The associated alignment matrix.B, in Eq. (8) is 



‘ 1 

0 

0 

0 

130 

B - 

0 

1 

0 

-130 

0 


. 0 

0 

1 

32.5 

18.75 


Figs. 6 show the control torques, forces and atti- 
tude angles for the "X-axis slewing" motion. The 
slewing time tf Is greatly shortened, tf * 3.988 
sec. (about one third of the slewing time without 
the forces, f x and fy ). 

Figs. 7 and 8 show the controls and attitude 
angle changes for the "Z-axIs slewing." For the 
sake of comparison, we use 2 different t f in the 
computation, t f « 27.5 sec., and tf* 20.0 sec. 
(minimum time; recall that tf *31.33 sec. without 
f x and fy). By comparing Figs. 7 with Figs. 8 
we can see that the control torques approach the 
bang-bang type when the slewing time Is shortened, 
and the maximum amplitude of the control forces 
Increases gradually. From Fig. 7c and Fig. 8c, we 
can also see the obvious Increases In e x and e y . 

This Is due to the increases In u x , u y , f x ,and J f y . 

(d) Now we consider a general case. Suppose 
the SCOLE is in an Earth orbit and we need the line 
of sight to be directed toward the center of the 
Earth. The orbital coordinate system (x,y,z) is 
shown In Fig. lb. Suppose, before the slewing, the 
Y axis of the spacecraft coincides with the orbital 
y axis, and the angular difference between X and x 
(or Z and z) axes is a* 7.897224212 deg. Thus, the 
Initial attitude quaternion of the spacecraft is 
q(0)= [cos ( a/2 ) 0 sin(a/2) 0]^. According to Ref. 

1, the unit vector along the line of sight in the 
rigid SCOLE coordinate system is R LQS 

* [0.1112447155 -0.2410302170 0.9641208678] 1 

The direction cosines of the orbital z axis in the 
body system at the initial time are z/ B = [sina 0 
cosa]. The angle between R LQS and z/g at the 
initial time is e LOS (0)= R L0S - z/ B = 20 deg. The 
eigen axis of the expected rotation in the body 
system is determined by 

e = ^L0S X Zy V^L0S x */ B l 
Thus, the quaternion for this rotation is 

q c = cos(20°/2 ) , q 1 = c i sin(20°/2), i = 1,2,3 

From 15) we can get the final attitude quatern- 
ion, q(tf). Fig. 9 shows the control torques, 
reflector forces, and attitude angles for this 
slewing motion. The 0 lo$ in Fig. 9c is the angle 
between the line of sight and the line of the 
target direction (from the spacecraft to the center 
of the Earth) . 

6. Concluding Remarks 

(1) There is a good agreement between the guessed 
value of tf and the value of tf to which the 
algorithm converges in the case (b). 

(2) The guessed initial values of the costates 

here: p(0), r(0) are adequate for the algo- 

rithm to converge. If the slewing time, tf, 
is sufficiently larger than the minimum time. 
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tJ, then, the converged values of p( 0) and 
r f 0 ) are very close to the guessed values and 
less number (4 times) of iterations Is needed 
(Fig. 4). The same situation was observed 
In Ref. 8 

(3) The control profiles obtained in this paper 
give us a good reference for future use. 

For example, an extension to the minimum 
time slewing motion of the SCOLE model con- 
taining both rigid and flexible components 
Is planned. 
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Appendix I The Term g (u> , r) in Eg. (16) 


The term I~1 q>Iu> in dynamical equation (9) can 
be replaced by 


[aj.j a! 


^ 2^3 




= [ A:B] u 


where the a^j and b-jj (the elements of 3x3 matrices 
A and B) areconstants associated with the Inertia 
parameters of the spacecraft. Then the term 
r^l''ulw of the Hamiltonian, H, in Eq. (14) has 
the form 

h*r^I ^ojIuj = [R-j R2 R3 R^ R^ Rg] u 


where R^ are 

[R^Rj] 1 = fl T r, [R^Rg] 1 * B T r 
The term g(o>, r) in Eq. (16) is obtained by 


g(w,r) = -(ah/3o>)* - 

2R i 

R 6 

R 6 R 5 ' 
2R 2 R 4 


'“1 ' 
“2 


R 5 

R 4 2R 3 


.“3 . 


Appendix II Solution of Eg. (44) 


Eq. (44) can be rewritten as 

a i e = b.e^ + c.t i (44) 

For simplicity, we only consider the solutions for 
the following boundary conditions 

e(0) = 0, 0(0) * 0; e(t f ) = 0 *, e(tf) * 0 

(IH) 

Suppose a* f 0, b.^0 and let b=bj/aj, c=c^/a.j 
(suppose c>0), we can rewrite Eq. (44) as 

8 = be 2 + ct ( 11-2) 

Since the control for this problem is of a bang- 
bang type with only one switching point, then, 
by integrating Eq. (11-2) and using (II-l), we get 


e = /c(e 2u6 -l )/b, for t=1 ; ( 1 1 - 3 ) 

e = /c(l-e^(e-®*))"/b # for t* -1 (1 1-4) 

By equating Eqs. (II-3) and (II-4), we get e=e s 
and 6=e s at the switching point, t=ts, 

e s - { 1 /2b)1og[2/(l +e _2be * 0 s =^c(e 2be s -1 )/b 

Finally, by integrating (11-3,4) and using (II-l), 
we get 

t cosh~^ (e’ be s) / ^-bc, b<0; 

[(ir/2)-sin _1 (e _be s)]//bc7 b>0 
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l f= 


f t s + C(Tt/2)-sin- 1 (e b (°* e s5 )//Tb^, b<0 ; 

[ t.+cosh* 1 (e b (°^ 0 s))//bc’ i b>0 


For the case 0(0) f 0, more complicated solutions 
can be obtained, but are not given here. 





Fig. la Spacecraft Control Laboratory Experiment Configuration (SCOLE) 


TOWARDS 



Fig. lb. Attitude of the SCOLE Showing Antenna Line of Sight 


(b) TIME (SEC) 


(b) TIME (SEC) 


2. X-Axis Slewing, t f = 12.57 Sec. 


Fig. 3. Z-Axis Slewing, = 31.33 Sec 
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ABSTRACT 


Tne study of the dynamics of the Spacecraft Control Laboratory Experi- 
ment (SCOPE) 1 is extended to emphasize the synthesis of control laws for both 
the linearized system as well as the large amplitude slewing maneuvers required 
to rapidly reorient the antenna line of sight. For control of the system 
through small amplitude displacements from the nominal equilibrium position 
If)R teclmiques are used to develop the control laws. Pontryagin's maximum 
principle is applied to minimize the time required for the slewing of a 
general rigid spacecraft system. Tie minimum slewing time is calculated 
based on a quasi- linearization algorithm for the resulting two point boundary 
value problem. 2 Tie effect of delay in the control input on the stability of 
a continuously acting controller (designed without considering the delay) 
is studied analytically for a second order plant. System instability can 
result even for delays which are only a small fraction of the natural period 
of motion. 3 
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Reflector 



A. 

l 0 Geocenter 


figure i i i - i ; 


the 3-d geometry of the scole 
CONFIGURATION IN ITS DEFORMED STATE 


. Equations of Notion 

‘Derived using a Newton-Euler approach 
. Assumptions 

- Reflector and Shuttle rigid 

- Mast has constant cross-section 

- It is assumed to undergo small elastic deformations only 

- Its modal shapes in orbit are assumed to be the same as 
those of an identical' non-rotating beam. 

Stability Analysis (Rigidized SCOLE) 

A stability analysis of the r igidzed SCOLE was conducted for 
the following configurations: 

a) Rigid - no offset . Pitch motion decouples from roll 
and yaw. in the linear ranges. System not stable 

b) Rigid - with offset parallel to roll axis. Pitch motion 
still decouples from roll and yaw in the linear range. 
System unstable. 

c) Rigid - With both offsets (parallel to roll and pitch 
axes). The motions in all 3 degrees of freedom are 
coupled. System found to be unstable. 

. Control Laws • 

Assumption : All the states of the system are available. 

It was suggested by J.G. Lin that an intuitively appealing 
practical approach to achieve the LOS pointing objective 
is a two-stage procedure, (a) Slew as if rioid then, 

(b) damp-out flexible dynamics. 

.The linear regulator theorv used here to control 
-the linear model of the rigidized SCOLE, 

-The linear model of the actual SCOLE configuration including 
the first foi.r flexible modes of the mast. 

Next 

Preliminary slew maneuvers o f rigidized SCOLE. 
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RIGID SCOLE - CONTROL EFFORTS 


a 



iso 


i 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 




3M 




SCOLE: TRANSIENT RESPONSES 

6.0 DEG. INITIAL PERTURBATION IN PITCH B&ltiJ&Uu h 



TIME (SECONDS) 




SCOLE: TRANSIENT RESPONSES 

B.O DEG. INITIAL PERTURBATION IN TAN fttOM £^Ut uMij* 



TIHE (SECONDS) 



LINEAR MODEL OF SCOLE WITH FLEXIBILITY 

HUNS. IBP. TO A Meg. PERMB.DI PITCH EWUiYi&fUu/^ 





LINEAR MODEL OF SCOLE WITH FLEXIBILITY 



IF7 







LINEAR MODEL OF SCOLE WITH FLEXIBILITY 

TUANS. HEBP. TO A l.Bft.OIST.IN MODE #1 





LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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LINEAR MODEL 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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MINIMUM TIME ATTITUDE SLEWING 
MANEUVERS OF A RIGID SPACECRAFT 

OBJECTIVE 

DEVELOP COMPUTATIONAL TECHNIQUES TO SLEW 
A GENERAL RIGID SPACECRAFT ( INCLUDING 
RIGIDIZED SCOLE ) FROM AN ARBITRARY INITIAL 
ATTITUDE TO A FINAL REQUIRED ATTITUDE 
PRECISELY, AND SATISFYING THE FOLLOWING 
CONDITIONS: 

. IN MINIMUM TIME 

. THE CONTROLS HAVE SATURATION LEVELS 



MET H ODOLOGV 


THE MAXIMUM PRINCIPLE FROM OPTIMAL. CONTROL THEORY IS 
APPLIED TO THE EULER'S DYNAMICAL EQUATIONS AND THE 
QUATERNION KINEMATICAL EQUATIONS OF THE SYSTEM TO 
DERIVE THE NECESSARY CONDITIONS FOR THE CONTROLS. 

THIS LEADS TO THE TWO- POINT BOUNDARY- VALUE PROBLEM. 

AN INTEGRAL OF A QUADRATIC FUNCTION OF THE CONTROLS 
IS USED AS A COST FUNCTION, BUT THE INTEGRATION 
PERIOD OF THIS INTEGRAL, CALLED THE SLEWING TIME, 

IS TO BE CHANGED UNTIL IT REACHES ITS MINIMUM VALUE. 

THE RESULTING TPBVP IS SOLVED BY A QUAS I LI NEARI ZAT I ON 
ALGORITHM < METHOD OF PARTICULAR SOLUTIONS >. 

EULER'S EIGENAXIS ROTATION THEOREM IS USED TO 
APPROXIMATELY DETERMINE THE INITIAL VALUES OF THE 
COSTATES AND THE SLEWING TIME AS WELL AS THE 
NOMINAL SOLUTIONS WHICH ARE UESD TO START THE 
QUASILINEARIZATION ALGORITHM. 

THE MINIMUM SLEWING TIME IS DETERMINED BY SHORTENING 
THE TOTAL SLEWING TIME UNTIL AT LK’*ST ONE OF THE 
CONTROLS BECOMES A BANG- BANG TYPE. 
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CONCLU D I NG REMA RKS 


. THE SLEUING MOTION NEED NOT BE RESTRICTED TO A 

SINGLE AXIS MANEUVER. 

. THE GUESSED STARTING VALUE OF THE SLEWING TIME IS 

VERV CLOSE TO THE CONVERGED VALUE FOR THE SCOLE 
EXAMPLES AND SUBROUTINE USED HERE. 

. THE GUESSED INITIAL VALUES OF THE COSTATES ARE 

ADEQUATE FOR THE ALGORITHM TO CONVERGE. 

. THE METHODS USED HERE MAV BE IMPLEMENTED FOR 
PRACTICAL CONTROL SOURCES WHICH MAT HAVE MORE 
CONSTRA I NTS . 

. AN EXTENSION TO THE MINIMUM TIME SLEWING MOTION 
OF THE SCOLE MODEL CONTAINING BOTH RIGID AND 
FLEXIBLE COMPONENTS IS PLANNED. 
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Appendix - Chapter II 

Stability Analysis of Second Order System 
with Delayed State Feedback 

As a second order differential equation describes the dynamics of 
a single mode of any large space structure, the stability analysis of such 
a system with delayed state feedback is analyzed and the amount of delay 

that can be tolerated by the system without becoming unstable is arrived 
at analytically. 

The differential equation of second order with state feedback can be 
written as : 

*i + 2c i“A + “i^i * -Vi(t-h) - kpX^t-h) (i) 

where 

= i^ 1 modal coordinate 
* i^ 1 natural frequency 
Cj/ - i^ 1 mode inherent damping ratio 
kj. = rate feedback gain 
kp 3 position feedback gain 
h = time delay 

1he feedback gains k,., kp are designed for the required stability and 
transient response specifications without taking the delay into consideration. 

The inherent dasping ratio, q' and the feedback gains, k r and 1^; will 
give rise to five possible combinations as shown in Table 1 and are thus 
analyzed separately for mathematical convenience and easy understanding. 
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Case I: 


‘i- 


0, kp = 0 and k r > 0 


The differential equation of the system can be written as: 


x i + “i x i = - k A(t-h) 


( 3 ) 


Case 


Ci 


I 

- 0 

> 0 

» 0 

II 

> 0 

> 0 

35 0 

III 

= 0 

> 0 

> 0 

IV 

> 0 

= 0 

t 0 

V 

> 0 

> 0 

^ 0 


Note: The remaining three combinations are 

neither feasible nor of interest. 

Table 1: Feasible Combinations of cJ, k_ k 
for Stability Analysis 1 ** P 

and the corresponding characterstic equation is given by 

2*2 cK 

s + oj + 2? i (D i se" sn = o 


k r - 25 iV 


(4a) 


where 

tte value of h for which the roots of equation (3) cross the imaginary 
axis can be evaluated by substituting s = joj. 

Thus 

2 _ 2 

u i “ + j2s^a)sin£uh + 2 ? -uj .oicosuh = 0 , 

1 i (4b) 

For equation (4b) to be satisfied 

sincuh = 0 
and 


2 2 . 

“i cosuh = 0 


( 5 ) 


ivy 


Thus 


“ h = V 2 


h = 


u i C? i + 7 1+ ?i 2 ] 

Case II: 5 » >0, k r = 2^ and kp = 0 

The characteristic equation of the system described by equation ( 1 ) 
is given by 


(o^-w +2t i (u i 0)Sinuh) + j (Zq^u+Zs^wcosuh) = q 


Thus 


cosmh * -qWc. 

l * i 


and h » cos~lfC-i/Cj 

»i C ^ C|-C- 2 + /"l+c|-^2] 3 

For the case where 5^1 thesystem will always be stable since no value 
of h exists for which the roots of (7) cross the imaginary axis. A plot 
°f “i h veT5us 'i for various values of tj is shown in Figure 2.1. 
ease III: - 0, kp « k r > 0 

The characteristic equation is given by 

sZ+uZ+^g-sh^g-sh , g 

or (M^-ul+wkpSinwh+kpCoswh) * j ( wk^osuh-kpsinwh) = 0 (1( 

^ 1US tan o»h = 

^ “ 2 ' I [C2w?*k2) * /[k^^Mk^]] (1] 

Plots of hwj versus V-l for various values of kp/w? are shown in 
Figure 2.2. It can be seen here that these are many combinations of ^ and 
k r for which the roots of Eq. (10) can cross the imaginary axis - i.e. 
value of hwj which leads to instability. 
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Case IV: 


Thus 

and 


? i >o, k r = o, kp f o 
TTie characteristic equation is given by 

(u i 2 -w 2 +kpC0S( 1 th) * j(2c i ' v - kpSiftnh) . Q 


sinwh » j C , i “i" 

^tT~ 


( 12 ) 

(13) 


ID 2 » (|) 2 ‘ 


*i a- 2 C : j t 2] 

The plots of hwj versus v / 2 r P 1 U4J 

P> 1 V“l fOT various values of c ! are shown in 

Figure 2.3 1 

CaSeV: 'i » °. *f>0, kyj-p 

The characteristic equation is given by 
(“i 2 '" 2 *”^ sinwh+kpcoswh) 

* ; (2t! •J i uwuk T cosa,h-k p siiw,h) > o 

By equating the imaginary part to zero h u 

/ p«rz to zero, u>h can be evaluated 


as 


* u>h 


sin’* ^ 


2?^ U)jU) 

^ k p 2 +u) k r 2 . 


- tan 


-1 


uik_ 

C ^ 3 


(is) 


(16) 


after substituting «h in the real part of equation 

the following equation in the single unknown variable a, can be obtained 

2 2 i 

“i -u -^k^in (sin - 1 y-tan‘ 1 ( ^ 3 ) 


*P 


+ kp cos (sin'Vtan^C ^ = 0 

K P 


(17) 


J 75 ” 



Using equations (17) and (16), the limiting value for given values of 
k r* ^ 311(1 “i can be determined. As the equation (17) is nonlinear, 
numerical procedures may have to be used and thus the generalized plots 
similar to the other cases may be obtained; 
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Plots of hto^ vs correspondence to Case II with c •' 
as a parameter. ^1 
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Figure 2.2 Plot of vs k r /u^ corresponding to Cas e III ^r^i 
with kp/<o? as a parameter 
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Abstract 


Modal state estimation tests are performed on the SCOLE mast for 
the fixed Shuttle platform case. Kalman filter state estimation results 
from a five mode computer model of the SCOLE mast, developed from a 
finite element analysis, are compared with those state estimates 
obtained from laboratory tests. Two comparison runs are presented, one 
an excitation of the first two bending modes, another, an excitation of 
the first torsional mode of the mast. Results from both runs show poor 
agreement in modal estimation between the computer model simulations and 
the laboratory test data. At present, the reason(s) for this poor 
performance is unknown. Both the laboratory hardware and software and 
the computer model are being checked for possible sources of errors. 
Further computer simulations as well as laboratory testing will be 
performed. 
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MODEL OF SCOLE MAST 


O MODAL DATA FORM FINITE ELEMENT ANALYSIS OF MAST 

O FIVE DECOUPLED MODES (FREQUENCIES . 443-4 . 345 HZ) 

O ACTUATORS: FOUR JETS AND THREE REACTION WHEELS 

G SENSORS: SIX LINEAR ACCELEROMETERS AND 3 AXIS 
RATE GYRO 

G MODAL STATE AND OUTPUT EQUATIONS: 

X(k+1) = AX(k) + BU(k) 

Y (k) = CX(k) + DU(k) 

G STATES ESTIMATED BY KALMAN FILTER 
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KALMAN FILTER 


O EQUATIONS IN BASIC FORM: 

X(k+1) = AX(k) + BU(k) 

A #. A 

X(k+1) = X (k+1) + G ( Y (k) - Y (k) ) 

X - PREDICTED STATE 

A 

X - ESTIMATED STATE 
G - KALMAN FILTER GAIN MATRIX 


O KALMAN FILTER GAINS ASSUMED CONSTANT 

- SENSOR NOISE INTENSITIES ESTIMATED FROM 
MANUFACTURERS' DATA 


0 ABOVE FORM USED IN SOFTWARE FOR LABORATORY TESTS 


/e* 



SIMULATION AND TEST PARAMTERES 


RUNS USED FORCING FUNCTION TO EXCIT SCOLE MAST 

F = A sin (w T) 


RUN 

EXCITATION TIME 

FREQUENCY 

AMPLITUDE 

DURATION 

1 

10 

sec 

.443 HZ 

5.0 

30 

sec 

2 

10 

sec 

1.504 HZ 

2.0 

30 

sec 


O COMPUTER SIMULATIONS USED SAME VALUES 


IQS 




RUI 



MODEL ANi 



EAI 

VIBRATIONRL MODE, FREl 




SPEC 
* 1 


SCOLE 


i 







COMMENTS 


O LABORATORY AND COMPUTER SIMULATIONS OF MODE ESTIMATES 
ARE VERY DIFFERENT ! 


G ONLY GENERAL AGREEMENT BETWEEN LABORATORY AND COMPUTER 
SIMULATION IS IN WHICH MODES ARE "DOMINANT" FOR THE 
RESPECTIVE FORCING FUNCTIONS 


O CONFIDENCE IN LINEAR COMPUTER SIMULATION MODEL 
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SUMMARY 


O MORE LABORATORY TESTS ARE REQUIRED 


© RE-CHECK LABORATORY APPARATUS (SOFTWARE AND HARDWARE) 


© PARAMETER IDENTIFICATION AND RE-DEFINING MODEL OF SCOLE 
SIMPLE LINEAR DECOUPLED MODAL MODEL MAY NOT BE 
SUFFICIENT FOR PROPER LABORATORY EXPERIMENTS 
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SLEWING AND VIBRATION CONTROL 
OF 

THE SCOLE 


JIGUAN GENE LIN 

CONTROL RESEARCH CORPORATION 
LEXINGTON, HA U2i7:> 
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I, Introduction 


OUTLIfilfea 


* »<- 
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II. Structural Vibrations in SCOLE Excited by 

Time -Minimized Rapid Slewing 

-- Pang-Pause-Bang (BPR ' 1 Control (800 it) 
-- Bang-Bang (BB) Control (0, 80, 25 lbs) 

How bad? Always that bad? 

No Force means Least Excitation? 


III. Active Damping of BPB-Excited Vibrations 
using High-Perforamance Modal Dashpots 


Direct Velocity Output Feedback really cannot 

CONTROL EXCESSIVE VIBRATIONS EFFECTIVELY/ QUICKLY? 


IV. Comments 

V. Conclusions 
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VI. Recommendations 



I SCOLE Primary control task is: 

Rapidly Slew or change the Line-of-Sight (LOS), and 
Settle or Damp Structural Vibrations to a required degree 

I The Objective is: 

Minimize the time required to Slew and Settle, 

UNTIL LOS REMAINS WITHIN A SPECIFIED ANGLE. 

I 2-Stage approach: 

First: Slew the whole structure like a rigid body, 

-- IN A MINIMUM TIME, 

UNDER THE LIMITED CONTROL MOMENTS AND FORCES 

Then: Damp the excited structural vibrations 


I Some previous results on Stage-I design 



I Objective of Current Study 

Stage-2 design: Active Control of Excited Vibrations 


I9S 


Structural Vibrations Excited by 
Bang-Bang-Type Rapid Slew Maneuvers 



Note 1 , Time of application is 11734 sec/ only 35 . 32 % of the slew time, 


i?g 





















DEFLECTION FT 

20 . OG G .00 20.00 40.00 


-00 2.00 4.00 6.00 8.00 10.00 

TIME SEC 


X-flXIS 



'o'. 00 2'. 00 4 1 . 00 gVoO 8 1 . 00 


TIME SEC 


Fig. 3-3 Vibratory responses to Rapid Time-minimized Bang-Bang Slew: 80 lb; 
a. Line-of-sight error and Mast tip deflection. 










]. 00 


y . od 


8.00 12.00 

TIME SEC 


PITCH, 


TRW, 


IB. 00 20.00 



y . oo 


8.00 12.00 

" I ME SEC 


16.00 20.00 


Fig. 3-4 Vibratory responses to Rapid Time-minimized Bang-Bang Slew: 25 lb; 

b. Attitude deviations at the Shuttle (S) and the Reflector (R) ends. 


I Using no forces on Reflector DOES NOT mean least excitation! 


t If LOS error is the only concern. Stop using BOOlb force; 
Use 25lb (Vernier thruster level) instead. 

-- But time is equally important!!! 

I Additional time of various length is still required 
for damping out the excited vibrations. 


I VIBRATION CONTROL CHALLENGE: 

Can excessive vibrations, 

SUCH AS EXCITED BY THE 80OLB RAPID SLEWING, 
BE EFFECTIVELY SUPPRESSED TO A REASONABLE LEVEL 
QUICKLY, SAY, IN 5 SEC? 
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MODAL-DASHPOT DESIGN MD1 


Part 1 : Linear Velocity Feedback 


U 


1 

71 


r y i 
ribi 


I _ _ r I 

II b LVR y I 

lJ ei | Y 16I 

J L J 


>71 


Uo 

O I 

L J 


l =F 4 = 


APPLIED FORCE ON REFLECTOR IN X~DI RECTI On] 

I 

APPLIED FORCE ON REFLECTOR IN Y-DIRECTIONj 


r Y, .1 r 


l 


I '151 RATE OF XZ-DEFLF.CTION AT REFLECTOR END | 

1 1 = I | 

jYjgj Irate of YZ-deflection at Reflector endI 

L J L J 


L VR 


, 53420630E+01 , 43392044E+00 

, 42058249E+U0 , 697S6355E+01 


Additional damping ratio designed = 0.6737, Mode 2 

=0.6, Mode 1 


2% SETTLING TIME OF 3 SEC IS DESIGNED FOR MODE 2 
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MODAL-DASHPOT DESIGN MD1 


Part 2: Angular Velocity Feedback 


V 

r Y i 
1 Y 10 
i 

% 

i Y n 

i 


= - G 

“6 

U AVR 1 v 

1 Y 12 

J 

L J 


V 


APPLIED MOMENT ON REFLECTOR ABOUT X”AXIS 

% 

■»4 ' 

APPLIED MOMENT ON REFLECTOR ABOUT Y~AXIS 

“6 

J 

APPLIED MOMENT ON REFLECTOR ABOUT Z"AXIS 


“ - 



YlO 


rate of Reflector Roll attitude deviation 

Y n 


rate of Reflector Pitch attitude deviation 

y 12 


rate of Reflector Yaw attitude deviation 

_ 


^ - 


G 


AVK 


, 24172707E+0A 
. 15734103E+03 
, 13433660E+04 


. 16653096E+03 
•21781213E+U4 
- , 22055215E+04 


■45158162E+03] 

I 

-, 72768193E+03 I 

I 

. A2951681E+04 j 


Additional damping ratio designed = 0,03/ Modes 3, A, 5 
Inherent damping ratio assumed = 0,003 all modes 
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LOS ERROR DEG 
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F0010 


O 

m 



Fig. 6-1 Simulation results of vibration control design MD1; 

b. Histories of la ne-of-s ight error and Mast tip deflection. 







COORDINATE NORMAL COORDINATE *10‘ 

3 0, GO 8. DO -80.00 -40.00 0.00 40.00 



•rSaS ** 1 


MODE 1 



MODE 2 




Summary — MD1 



LOS error < 18,46° t - 2,5 sec 

< 17,46° t > 2.5 sec 

= 11,79° t = 3,1 sec 


Deflection < + 5 ft t - 2 sec 

< + 0,5 ft r = 4,2 sec 

2% settling time <, 2,9 sec 


Last peak Reflector attitude deviation: 

Roll = 0,460° 

Pitch - 0,546° 

Yaw - 1,360° 
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Summary -- MD1A 


Additional damping ratio re-designed = 0,6/ Mode 2 
CORRESPONDING 2% SETTLING TIME FOR MODE 2 IS 5.38 SEC 

|" , 58920557E+01 .95784262E+00 

6 l\/r " | . 42061494E+00 , 62209375E+01 

L J 

LOS error <. 16.66° ( if ) t 2. 1-8 sec 

«= 9.57° (+) t = 3.1 sec 

Deflection <, + 7.35 ft (♦ ) t = 1.3 sec 

<, + 0.75 ft (4 ) t = 3.7 SEC 

2% settling time = 3 sec 


Last peak Reflector attitude deviation: 

Roll =0.719° (f) 

Pitch = U.582° 

Yaw = 1.399° 
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COMMENTS 


9 The Modal Dashpot design met the vibration control challenge 

FAIRLY WELL: EFFEC TIVE/ FAST SUPPRESSION OF EXCESSIVE VIBRATIONS 

9 For complete suppression and precision pointing 
after the quick suppression. 

Either: increase the Modal Dashpot feedback gains 

Or: Switch to integrated design of LQG/LTR and Modal Dashpots 

9 Direct Velocity Output Feedback Controllers 
need NOT be of "low authority", low performance, 

-- additional damping ratio can be designed to be as high 
as to the optimal value 0,707, if neacessary; 
instead of restricting to only about 0.1 

9 No more high-gain problems of original Canavin dedign 

9 Spillover is minimal: performance degradation un-noticeable 
Spillover is beneficial: concomitant active damping of unmodeled 
modes 

9 Systematic design method for modal dashpots works! 
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CONCLUSIONS 


• 2-Stage Approach is feasible and promising for 

RAPID SLEWING AND PRECISION POINTING OF SCOLE 

I Not all bang-bang type of time-minimized slew maneuvers 
WILL EXCITE LARGE STRUCTURAL VIBRATIONS IN SCOLE 

I Modal dashpots can be a concentrated high-power vibration control/ 

AS WELL AS THE USUAL 

DIFFUSE ("BROAD-BAND")/ LOW“POWER ( “LOW - AUTHOR I TY ” ) CONTROL 


RECOMMENDATIONS 

i Limit the magnitude of applied forces on Reflector to 
EITHER 25 LB 

-- LEVEL OF VERNIER THRUSTERS ON THE REAL SPACE SHUTTLE 

OR 150 LB 

— LEVEL EQUIVALENT TO THE COLD'GAS JETS OF LABORATORY SCOLE 

• To complete Stage 2, add an integrated design of 
LQG/LTR (Li near -Quadratic-Gauss i an/Loop-Transfer -Recovery) 

and Modal Dashpots 

• Validate the 2-Stage approach using the SCOLE laboratory facility 

WITH A COMPREHENSIVE SEQUENCE OF INTEGRATED DESIGNS AND EXPERIMENTS 
COUPLING NONLINEAR RIGID-BODY MOTIONS WITH FLEXIBLE-BODY DYNAMICS 
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PLACING DYNAMIC SENSORS AND ACTUATORS 
ON FLEXIBLE SPACE STRUCTURES 


Gregory A. Norris and Robert E. Skelton 


Purdue University 


ABSTRACT 

Input/Output Cost Analysis involves decompositions of the quadratic cost function 
into contributions from each stochastic input and each weighted output. In the past, these 
suboptimal cost decomposition methods of sensor and actuator selection (SAS) have 
been used to locate perfect (infinite bandwidth) sensors and actuators on large scale 
systems. This paper extends these ideas to the more practical case of imperfect actuators 
and sensors with dynamics of their own. NASA’s SCOLE examples demonstrate that 
sensor and actuator dynamics affect the optimal selection and placement of sensors and 
actuators. 
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1.0 INTRODUCTION 


The objective of this paper is to develop and evaluate a method for the selection of 
sensors and actuators in the control of finite-dimensional linear systems using imperfect 
sensors and actuators — devices which do not provide instantaneous responses, but have 
nontrivial dynamics of their own. In addition, the plant noise and the measurement noise 
is assumed correlated. This important case allows the use of accelerometers as sensors 
(this always yields correlated plant and measurement noise). Application of the 
generalized method to practical control problems demonstrates that correlatedness of the 
noise and the dynamics of the actuator and sensor devices can significantly affect the 
optimal selection of both the number and location of sensors and actuators. 

Consider as a starting point the following familiar dynamic system model: 


Xp = A p Xp(t) + B p [f(t) + w(t)] 

(1.1a) 

Yp(0 = C p x p (t) , z(t) = M p x p (t) + v(t) 

G.lb) 


E [w(t)w T (x)} = S(t-x)W , E {v(t)v T (x)} = 8(t-t)V , E {w(t)v T (t)} = 5(t-x)U(i.lc) 

where x p <= R n ‘, fe R"*, we R n *, z,v, e R" 1 and (A p ,C p ) observable, (A p ,B p ) 
controllable and (Ap,Mp) detectable. The vectors w(t) and v(t) are respectively zero 
mean white noise characterizations of the actuator and sensor noise. 

In control of large space structures, the locations of sensors and actuators becomes a 
critically significant "degree of freedom" in control design [14, 20]. Among over 60 
more recent contributions to the SAS problem, only [4], [7], [10], [1 1], and [12] consider 
noisy actuators (W, V nonzero). In all cases, the disturbances are modelled as Gaussian, 
white, and uncorrelated (W, V diagonal, U = 0). Most of the SAS literature takes no 
account of actuator or sensor dynamics. Two exceptions are McClamrock [19], and 
Howell and Baxter, [6]. In [1] the authors extend the cost decomposition approach [2] to 
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accommodate noise correlation between sensor and actuator noise sources (W, V not 
diagonal, U * 0). A key conclusion in [1] is that the proper sensor/actuator selection and 
placement can be drastically affected by noise correlation. For example, the deletion of a 
noise source (by making an actuator or sensor noise free) may degrade performance 
contrary to the usual expectations when noise sources are uncorrelated. 

Very fast actuator dynamics may be neglected in stability considerations, [9]. A 
more thorough discussion of the effect of actuator dynamics is given by Goh and 
Caughey [8]. The analysis of [8] and [9] demonstrates that plant frequencies occurring 
above the actuator bandwidth can lead to closed loop instability, even for co-located 
sensors and actuators. Goh and Caughey do not address the problem of selection of 
dynamic actuators. That is the goal of this paper. 

This paper is organized as follows. First the system model is augmented to include 
sensor and actuator dynamics. The closed-loop input and output costs are then developed 
for the fully augmented system, and they are used to define expressions which reflect the 
effectiveness of each dynamic actuator or sensor in minimizing the cost function. 
Finally, the method is illustrated by application both to small scale numerical examples 
and to NASA’s SCOLE flexible space structure model. It is found that in the selection of 
noisy actuators and sensors, finite dynamics can significantly affect selection results. 

2.0 MODELING DYNAMIC ACTUATORS AND SENSORS 

In [2] the results of Closed-Loop Input/Output Cost Analysis (CIOCA) were 
developed and applied to the Sensor and Actuator Selection problem (SAS) for systems 
of the form (1.1) under closed-loop control. In [1] the control f(t) is the vector of optimal 
state estimate feedback controls: 
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(2.1a) 


f(t)=Gx c (t), G = -R- 1 B p T K. 

x c = ApX c (t) + B p f(t) + F[z(t) - MpX c (t)] , F = [PM^ + B p U]V~ l , (2.1b) 

0 = KA p + A p T K - KBpR^BpT’K + C p T QC p (2.2a) 

0 = [A p - B p UV -I Mp]P + P[A p - B p UV _1 Mp] T - PMj’V^P (2.2b) 

+ B p WBp T -BpUV- 1 U T Bp T 
which minimizes the cost function 

V = £»{ Hy p (t)ll 2 Q + llu(t)ll 2 RJ , EJ^lim E[-] (2.3) 

t— 


where x c € R n * is the vector of state estimates. The conclusion from [1] for this problem 

(1.1), (2.1) (2.2) is that when U^O, the sensor/actuator selection results can be 
drastically different. 


2. 1 Adding Actuator Dynamics 

First the system (1.1) is augmented to include stable, observable, controllable 
actuator dynamics of arbitrary order. 

x a = A a x a f B a (u+w u ) , (A a ,B a ) controllable (2.4a) 

f = C a x a , (A a ,CJ observable, x a € R"\ f e R"* (2.4b) 

Figure 2.1 presents schematic representations for actuator models of varying degrees of 
complexity; Figure 2.1a represents the non-dynamic actuator, while Figure 2.1b 
represents the general model for a dynamic actuator with white noise. Note that for the 
non-dynamic actuator the noise w(t) is purely additive with the input u(t). In the case of 
dynamic actuators the analyst may consider the actuator’s output (into the system) to 
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Figure 2. lb 


Figure 2.1: Actuator Models 
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include additive actuator output noise w^t), or actuator command noise w u (t) which is 
filtered by the dynamics of the actuator, or both. Both types of noise are assumed 
possible in our development. 

Augmenting the system states x p of the original system (2.1) with a vector of 
actuator states x a , we obtain: 


x = Ax + Bu + Dw , y = Cx , z = Mx + v = z p 


(2.5a) 





cx 

v __ 


Wf' 

, A = 

A p 

B p C a 

X = 

. X a. 

»y = 


, W = 

. W 1 

0 

A a 


D = 


W = 


Wf u fu 
u f ; w u 


, c = 


Cp 0 

0 c. 


, B = 


0 

B, 


, M r = 


B p 0 
0 B. 


M p 

0 


where f = C a x a , (A a , C a ) is observable, and Re 
is controllable. 




<0, i = 1, 2, • • • n a . (A a ,B a ) 


First note that since x p is observable from y p , (i.e., (A p ,Cp) is observable) and x a is 
observable from f (i.e., (A a ,C.j) is observable) then from the definitions (2.5) x must be 
observable from y, that is: 


(A,C) is observable 


(2.6a) 


Also note that the actuator dynamics are assumed stable, so that the system (2.1) has 
not been augmented with any unstable states. Therefore, detectability of (A p ,M p ) 
together with stable A a yields 

(A,M) detectable. (2.6b) 
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Finally, Theorem 1 states the conditions for controllability of the system (2.1) 
augmented with actuator dynamics (2.5). Proof of the theorem is contained in the 
Appendix. 

Theorem 1 

Consider the controllable system 

x p = A p Xp + B p (f+w f ) , (A p ,Bp) controllable (2.7a) 

x p e R nXp (2.7b) 

augmented with controllable and observable actuator dynamics of arbitrary order 

x a = A a x a + B a (u+w u ) , (A a ,B a ) controllable (2.7c) 

f = C a x a , (A a ,CJ observable, x a eR nx ‘, fe R nu (2.7d) 

to form the composite system 

x = Ax + Dw + Bu (2.7e) 



Ap B p C a 


1 

DO 

XJ 

O 

i 

A = 

0 A a 

, D = 

r 

o 

DO 

£ 




The system states x p are controllable from u(t) if the number of poles minus the 
number of zeros is the same for each individual actuator’s transfer function. 


Remark 1: Note that full controllability of the augmented-system state vector x is 

not guaranteed under the conditions of the theorem. 


Remark 2: The conditions of the theorem are always met for first order dynamic 

actuators, (assuming no direct input/output "feedthrough'' for the 
actuators), since each actuator will have one pole and no zeros. 
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Remark 3: 


As long as the original system states Xp are controllable through some 
minimum set of actuators meeting the criteria of the theorem above, then 
controllability of Xp will be maintained with the addition of actuators of 
any order and any number of transfer zeros. 

Remark 4: Finally, note that the usefulness of the theorem stems from the fact that 

by meeting certain mildly restrictive conditions, the actuator dynamics 
can be guaranteed not to destroy controllability of the original system 
states, regardless of the pole/zero locations of the plant. 

2.2 Adding Sensor Dynamics 

Next the system (2.5) is augmented to include stable, observable, controllable 
sensor dynamics of arbitrary order. 

x s = A s x s + BsCMpXp+vJ, (A S ,B S ) controllable (2.8a) 

z = C s x s + v out , (A S ,C S ) is observable, x s e R"* (2.8b) 

Re[X i(A s )]<0, i = { 1, 2, ••• nj (2.8c) 

Figure 2.2 presents schematic representations for actuator models of varying degrees of 
complexity; Figure 2.2a represents the non-dynamic sensor, while Figure 2.2b represents 
the general model for a dynamic sensor with white noise. Note from both eqn (2.8) and 
Figure 2.2b that (similarly to the case of actuator dynamics), adding sensor dynamics 
leads to the possibility of both sensor input noise and sensor output noise. The i^ 
sensor’s input noise (v in )i is filtered by the dynamics of the i^ sensor, while the output 
noise (v out )j is purely additive with the sensor output. Both types of noise are assumed 
possible in our development. 
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Figure 2.2a 



Figure 2.2b 


Figure 2.2: Sensor Models 
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The fully augmented system equations have the following form: 


x = Ax + Bu + Dw 

y -Cx 

z = Mx + v 

* T = [x p T , xj, x s T ] , y T = [y p T , f T ] , w T = [w f T , w u T , vfl , v = 


(2.9a) 

(2.9b) 

(2.9c) 



1 

Cd 

O 
O t 


A p B p C a 0 


'c p 0 0 

D = 

o B a 0 

, A = 

O 

< 

o 

, c = 

0 C a 0 


1 

o 

O 

DO 




®sMp 0 A s 




or C = [C 0] 


B = 


0 





w 

u" 

B a 

- y=v out , vr = 

u* 

Vin 

0 





, M =[0 0 C s ] 


The response y (s) of the plant to the input f(s) is given by 


y p (s) = Hp(s)f(s) 


(2.10a) 


where 


Hp(s) = Cp(sl-A p r ‘Bp. (2. 10b) 

is the plant transfer function. The response Mx(s) of the actuator/plant system (2.9) to 
the input u(s) is given by 


Mx(s) = H(s)u(s) , 


(2.11a) 


where 


H(s) = M(sI-A) -1 B. 

Finally, the response z(s) of the sensors to an input Mx(s) is given by 
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(2.1Ib) 



z(s) = H^sJMxfs) , 


(2.12a) 


where 

H s (s) = C s (sI-A s ) _1 B s (2. 1 2b) 

is the transfer function for the sensor dynamics. Minimal systems are controllable and 
observable. Thus, given minimality of the plant/actuator system [(A,B) controllable and 
(A,M) observable], then measurability of the full augmented system is guaranteed 
[(A ,M ) observable] if there are no pole/zero cancellations between H(s) and H s (s). 

2.3 Defining the Cost Function 

With the properties of the augmented system established, optimal control design for 
the augmented system is now considered. Recall that the standard LQG cost function 
(2.3) for the unaugmented system (2.1) includes a penalty on the output regulation error 
y(t), as well as a penalty on the control energy u(t). However, in the augmented system 
(2.5), while the actuator command is given by u(t), the actuator response f(t) (contained 
in the augmented output vector y) is distinct from u(t) due to actuator dynamics. A true 
measure of control energy is more appropriately stated in terms of a weighted sum of the 
variances of f(t) rather than of u(t). It can readily be shown, however, that even in the 
presence of a weighting on the actuator outputs, f(t), some nonzero weighting on the 
actuator inputs u(t) is necessary to avoid an infinite gain solution to the optimization 
problem. For this reason, and in view of the relation of f(t) to the design goals as 
discussed above, minimization of cost functions of the form 

V - [ Hy(t)l Iq + IIu(0IIr] (2.13) 

and 
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Q = diag[Q 0 ,QJ , Q > 0 


(2.14) 


provides a stable optimal closed-loop solution. 


3.0 SELECTION OF DYNAMIC SENSORS AND ACTUATORS 
3.1 Closed-Loop Input/Output Cost Analysis 

In order to write the expressions for the closed-loop input and output costs, it is first 
necessary to put the fully augmented system, under closed loop steady-state optimal 
state-estimate feedback control, in the following state space form: 

x(t) = Ax(t) + Dw(t) (3.1a) 


y(0 = Cx(t) 


(3.1b) 


V=E„V Q (t), V 0 (t) = y*(t)Qy(t) , 


where 


x T = [*VL y T = [y ] T w T = [w T ,v T ] 


(3.1c) 


(3. Id) 



A BG 


D 0 


C 0 

A = 

FM A+BG-FM 

, D = 

0 F 

, c = 

0 G 


.Q = 


Q 0 

0 R 


, w = 


w 

u m 


u 

V 


(3. le) 


G = -R~ l B r K, 0 = KA +A J K - KBR~ x B t K + C t QC (3.10 

F = [PM T +DU]V~ l , 0 = [A-DUV~ l M]P +P[A-DUV- x M) t (3.1g) 


-PM t V~ x MP +DWD t -DUV~ x U t D t 


For the system (3.1) the output costs V, y , defined by 
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v* = (i/2){Ejavy8 yi ) yi } 


(3.2a) 


are calculated as follows [2] 

V> = [CXC T QJ ii 

where X is the stead y state covariance satisfying 

0 = AX + XA t + DWD t 

and where the output costs satisfy the cost decomposition propert y 

Tly 

EVJ = v. 
i=l 


The input costs are defined b y 


Y i w = (l/2){£ oo 0K 0 /aw i )w i } 


and are found from [2] 

V ; w = [D T SDW|jj 


where S satisfies 

0 = A T S + SA + C t QC 

and where the input costs also satisfy the cost decomposition property 

z v r = v' • 


(3.2b) 


(3.2c) 


(3. 2d) 


(3.3a) 


(3.3b) 


(3.3c) 


(3.3d) 


The input and output costs represent the in situ contributions that the noise inputs 
and the s y stem outputs make in the cost function. We ma y also wish to know the amount 
b y which the cost function will be reduced if a noise input is eliminated. This amount, 
AV*, is defined as 
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(3.4) 


AV, W = V - V Ri 

where V Ri is the value of the cost function after the i^ noise input is eliminated, (but the 
controller is not redesigned) and AV* is the cost reduction due to eliminating A 
positive value for AV ; W indicates that elimination of the i^ input will reduce the cost, 
while negative AV^ indicates that a cost increase will follow noise elimination. It was 
shown in [1] that the AV; W may be positive or negative in the presence of noise 
correlation. Partitioning the matrices W and D facilitates direct solution for the cost 
reduction [2], yielding 

AVj w = 2V i w - d ‘Sd.Wii . (3.5) 


The closed-loop covariance X may be written 




P+N N 
N N 


where P satisfies eqn (3.1 g) and where N satisfies: 

0 = N (A +BG ) T + (A +BG )N + FVF 1 


(3.6) 


(3.7) 


Also, S has the following form 


S = 


K+L -L 
-L L 


where K satisfies eqn (3. If) and where L satisfies 


0 = L (A -FM ) + (A -FM ) J L + G J RG 


(3.8) 


(3.9) 


For notational convenience the steady state covariance X is partitioned as follows: 
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■ Xp X 12 x 13 

X=[P+N} = xj 2 X. X 23 (3.10) 

xj 3 x5 x s 

Using the notation of (3.10) and the special structure of the closed-loop system matrices 


in eqn (3.13) we write the following expressions for the output costs 

V,» = [CpXpC^Qplu i = 1, • • • ny p (3.1 la) 

V/ = [C a X a C a T QJii i = 1, • • • nu (3.1 lb) 

V? = [GNG r R ] d i = 1, nu (3.11c) 

and for the input costs 

V i w = [D T (/i:+L)DW , ] ii i=l, ••• nw (3.12a) 

V> = [D T (X+L)D^] nw+iinw+i i= 1, • • nz (3.12b) 

V?™ = [F T LFV] ii i=l, ••• nz (3.12c) 

and the input cost reductions 

AVi W = [D T (K+L)DW -D T LFU T ] ll i = 1, • • • nw (3.13a) 

AV^ = [D J (K+L)DW -D J LFU r ] nw + inw+i i = 1, • nz (3.13b) 

AV i v< ”“ = [F T LFV - F r LFV - F r LBU ] ll . i = 1, • • • nz (3. 13c) 
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3.2 Dynamic Actuator Effectiveness Values 


Now that the closed-loop input and output costs have been determined for systems 
with dynamic sensors and actuators, it remains to use the CIOCA results to define 
expressions which reflect the effectiveness of each sensor and actuator in the cost 
function. This section defines the effectiveness values for dynamic actuators. The 
approach taken in [1] and [2] for non-dynamic actuators was to subtract the contribution 
the i^ actuator’s noise in the cost function from the contribution of its control signal, and 
to label this difference the "effectiveness" of the i^ actuator, Vj act . That is, 

V i act = V i u -AV i w (3.14) 

This subtracts the "bad" from the "good" contributions of the actuator to measure its 
effectiveness. The results of applying (3.14) to sensor and actuator selection for a range 
of small and large scale examples in [21, [3], [4], [17] and [181 have demonstrated the 
utility of this approach. 

Extending the definition (3.14) for applicability to systems with dynamic actuators, 
we proceed as follows. In (3.1) there are two noise sources associated with each 
actuator: command noise, w u , which is filtered by the actuator dynamics; and output 
noise, w f , which is additive with the actuator output. Thus, the noise contribution 
associated with the i^ actuator is given by the sum of AV i w * and AVj Wf . 

The beneficial control cost for each actuator is not immediately evident. First, 
recall that it is the actuator output f(t), not its input u(t), which drives the system. Next, 
note that the contribution of the i^ actuator’s output in the cost function, V*, includes the 
effects of noise w ui . That is, even in the open loop (u=0), v/*0 for [W u ljj>0 with 
dynamics. Hence, to define the beneficial (control) portion of V* it is necessary to 
subtract the portion of V/ which is due to noise. This can not be accomplished exactly, 
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since the actuator command u(t) and the command noise w u (x) are correlated for t > x. 
An approximation is obtained, however, by solving for V/ when u = 0 (that is, in the 
open loop). We define the contribution of w ui to V* and the contribution of Uj to V/ as 
follows, using the open loop covariance of the actuator states X„ : 

[V i f ] w = [C a X a C a T QJ ii (3.15a) 

and 

[Vi f J u = V/ - [V/] w = [C a (X a -X a )C a T QJ li (3.1 5b) 

where Xa solves 

0 = AaX, + XaA a T + B a W u B a T . (3. 15c) 

Finally, the input costs and the decomposition of the output cost V/ are combined in 
an effectiveness formula for dynamic actuators which is motivated by the results of [1] 
and [2]: 

Vr t = [V i f J u -AV i wr -AV i w \ (3.16) 

Note that in the absence of command input noise, [V/] w and Vj w * are both zero. Also, in 
the absence of actuator dynamics, f ; (t) is equivalent to uj(t). Thus the expression (3.16) 
reduces to the original effectiveness formula of [1] in the absence of actuator dynamics. 
Note also that (3.16) is applicable whether or not the actuator noise signals are correlated 
with other noise sources, and it is applicable to systems with actuator dynamics of 
arbitrary order. 


3.3 Dynamic Sensor Effectiveness Values 

Unlike the actuator noise, (which has a direct path to the output, independently of 
the conrollers influence) the noise associated with sensors reaches the system only 
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through the controller. Since the gains in the Kalman filter of the LQG controller 
represent an optimal trade-off of each sensor’s (beneficial) measurement information 
versus the (performance degrading) impact of its noise, then a AV* of large magnitude is 
indicative of a highly effective sensor. That is, the fact that a sensor’s noise is being 
allowed to heavily affect the cost means that its measurement information is even more 
critical to performance. For this reason, the following effectiveness formula for non- 
dynamic sensors, generalized to accommodate the possibility of noise correlation, was 
presented in [1]: 

V i sen = I A V I . (3.17) 

For dynamic sensors there are two possible noise inputs associated with each sensor. 
As in the non-dynamic case, both noise inputs reach the system dynamics through the 
Kalman filter. Thus a straightforward extension of (3.17) to dynamic sensors is 

Vj sen = I AV; Vin 1 + I AV i v “' I . (3.18) 

Note that this formula is applicable in the presence of sensor dynamics of arbitrary order, 
and applies whether or not any of the noise sources are correlated with one another. 

This section concludes with the suggestion that (3.16) and (3.18) provide effective 
measures of tire contribution of each actuator and sensor in a closed loop optimal LQG 
control (with sensor and actuator dynamics properly included). 


4.0 SELECTION OF DYNAMIC ACTUATORS FOR SCOLE 

In this section the actuator selection problem is solved for a model of NASA’s 
SCOLE (Spacecraft Control Laboratory Experiment) system. The SCOLE configuration 
consists of a flexible antenna suspended from the Space Shuttle cargo bay by a 130 ft. 
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flexible beam (see Figure 4.6). The effectiveness values for proof mass actuators 
(PMA’s) located along the beam are calculated and plotted versus position for both 
dynamic and non-dynamic actuators in order to evaluate the dynamic actuator selection 
method and to determine the effect of actuator dynamics on our results. 

4.1 PROBLEM DEFINITION 

A certain 2-dimensional SCOLE model includes four flexible modes and no rigid 
body modes [15-18]. Approximate open-loop mode shapes for the four flexible modes 
are presented in Figure 4.7, and the results of an open loop modal cost analysis are 
presented in Table 4.2. A detailed discussion of the model development is given in [16] 
and [18]. The two sensors retained in the model (using the CIOCA method of selection 
for non-dynamic sensors) are angular position and rate measurements located at the 
center of mass of the reflector [18]. Since there are no accelerometers presents, then the 
sensor and actuator noise is uncorrelated. Noise intensity data for the sensors is given in 
Table 4.3. 

The set of admissible actuators includes both a control moment gyro (CMG) located 
at the reflector center of mass and a set of PMAs distributed along the flexible beam. The 
actuator selection problem is to determine the optimal location for two PMA devices 
along the beam. To this end, the admissible set of PMAs was defined as 20 actuators 
spaced at distances of 6.25 feet apart on the 130-foot beam from a point 10.75 feet above 
the shuttle end of the beam to a point 129.5 feet from the shuttle. The PMA locations are 
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Pilch Ik mling mode ft J. Fitqm.ncy = .29 Hi 



Pjlch landing inode ft 2. Frequency^ 1.65 Hi 



Pjich Lending mode ft 3. Frequency^ 4.97 Hi 



Pitch landing mode ft 4. Frequency^ 12.36 Hi 



Pitch bending mode ft 7. Frequency^ 57 92 Hi 



Pilch Lending mode ft 8 . Frequency— 80.73 Hi 
Figure 14: Flexible Pitch -Direction Mode Shapes for 2 Dimension a! SCOLH 
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Tabic 4.2: 2-Dimensional SCOLE Elastic Modal Cost Analysis 


Mode# 

Frequency (Hz) 

Modal Cost 

Percent Total 

1 

.2S9E+00 

.547E+01 

466E+02 

2 

. 164E+01 

626E+01 

.533E+02 

3 

.4 97 E +01 

105E-01 

895 E *01 

4 

.124E+02 

. 1 80E-04 

153E 03 

5 

237E+02 

201E-06 

.171E-05 

6 

389E+02 

.692E 08 

590E-07 

7 

5S0E+02 

.470E-09 1 

,4fK)E-0S 

S 

SK-E+02 

.51* E 10 

.431 E 09 

0 

. in*E *-03 

.74SE-1 1 

fC'SE 10 

10 

.1 3’iE^03 

U2E-11 

I21E-10 

11 

. 175E+03 

.322E-12 

274E-1 1 

i: 

21 SE-r-03 

; h ; >b-i? 

.72 3 E 12 

13 

259E+03 

1 .243E- 1 3 

207E-12 


Tabic 4.3; Noise Specifications for SCOLE Actuators and Sensors 

Actuators 


Dynamic Noise 


Nooe 


Type Range lmen>;;v Type 


L. tensity 


PM A 10 lb .0001 (lb) 2 Accelerometers v, ■= ()02S (deg/sec 2 ) 2 

CMC 10 5 ft— lb 10,000. tft-lb) 2 Angular Pos .non 10 * (deg) 2 

()04 Meg/sec) 2 
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thus selected by evaluating the relative effectiveness of each of the 20 PMA locations. 


4.2 RESULTS FOR NON-DYNAMIC ACTUATORS 

The PMA selection problem for non-dynamic actuators was solved first, for later 
comparison with the dynamic actuator selection results. In all cases (dynamic and non- 
dynamic) the actuator effectiveness values are calculated following controller design 
which achieves a specified output variance and minimizes the amount by which the 
actuators exceed their specified variances. This type of controller is designed by an 
iterative selection of the control and output weights using the Output Variance 
Assignment (OVA) algorithm (DeLorenzo and Skelton, [3]). The variance specification 
for each actuator was equal to 10 times the intensity of its noise (see Table 4.3). 

The actuator effectiveness values based upon standard Closed Loop Input/Output 
Cost Analysis (CIOCA) [2] for non-dynamic actuators are presented in Figure 4.8. The 
figure portrays PMA effectiveness results for four different controllers, each achieving a 
different steady-state line-of-sight (LOS) error variance. The results provide a vivid 
illustration of how the controller objectives can profoundly influence the actuator 
selection results. For lower gain controllers (lower LOS error) the theory determines that 
the upper tip is the most desirable PMA location. However, as the gain increases 
(controller designed for smaller LOS error) the center of the beam becomes the optimal 
location. 

The results of Figure 4.8 are readily explained via modal analysis. The mode shape 
figures for the four flexible modes retained in the 2-dimensional SCOLE model were 
presented in Figure 4.7. Recall that mode #1, which accounts for 46.6 percent of the 
open loop modal cost, has a maximum amplitude at the reflector-end tip of the beam (i.e., 
at 130 ft.). Mode number #2, which accounts for 53.3 percent of the open loop modal 
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sigma 

0.09 
0.095 
«* 0.1 
•a- o 1015 


24 0 




cost, has a maximum amplitude near the 90 ft. point. And mode #3, which accounts for 
only approximately 0. 1 percent of the open loop modal cost, has a peak amplitude near 
the center of the beam. 

Next note from Figure 4.8 that as the gain is increased in order to achieve a smaller 
steady state LOS error variance, the most effective location for PMAs shifts from the tip 
of the beam to the midpoint. This corresponds to a shift from the peak of mode #1 to the 
peak of mode #3. The shift occurs even though with higher gain the noise in the PMAs 
near the beam midpoint becomes the most detrimental to performance (Figure 4.9). In 
fact. Figure 4.9 indicates the reason for the shift in optimal PMA location: with higher 
gain the third mode becomes the least damped by the control of the CMG, and becomes 
therefore a significant mode to be controlled by the PMAs. Figure 4.10 and Table 4.4 
indicate the motion of the closed-loop eigenvalues from their open loop locations under 
varying levels of gain (output performance). 

Since the control cost of each PMA (V ; u = E^qu 2 ) is equal to its effectiveness value 
V; act minus the cost contribution of its noise, V™, then it is clear from Figures 4.8 and 4.9 
that the PMAs are being used primarily to control mode #3 (i.e., near the middle of the 
beam). However mode #3 is the most lightly damped mode in the closed loop. This is 
true in spite of the fact that in all cases the input variances of both the CMG and PMAs, 
when normalized by their variance specifications, are of like order of magnitude (see 
Figure 4.11). 

The results demonstrate the interesting result that while the PMAs are being used at 
a level similar to the CMGs (in relation to their specified variance levels), they 
nontheless make a small contribution to the closed-loop modal damping. This claim is 
verified by deleting all PMAs from the system and again using OVA to achieve a 
specified LOS error of 0.1 (arc sec) 2 , and comparing the resulting closed-loop modal and 
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Tabic 4 4: Modal Characteristics of (A-t BGj as a Function of Output Performance 



Open 

Loop 

.1015 

(deg) 2 

.10 

(deg) 2 

.095 

(deg) 2 

.09 

(no PMAs) 

.10 (deg) 2 

O) |(r/s) 

1.8 

3.14 

3.14 

3.14 

3 14 

3.135 

Ci 

.005 

.028 

.027 

.0255 

.025 

.01 

tj (sec) 

1111 

11.4 

11.8 

12 5 

12.74 

31.9 

o >2 f r/s ) 

10.3 

56. 

62. 

93 

160. 

63 

r 

*a2 

(H>5 

.675 

68 

.6926 

702 

68 

T 2 (sec) 

19 4 

.026 

024 

.015 

009 

023 

CO3 f r/s ) 

31. 

29. 

28 9 

28.8 

28.8 

28 9 

93 

<X>5 

.0217 

.0185 

Oil 

.009 

.016 

tj(sec) 

6.45 

1.6 

1.87 

3.16 

3.95 

2.16 

co^f r/s) 

78 

77 4 

77. 

77. 

77. 

77 3 


.005 

,0065 

007 

.0072 

.0052 

.0068 

x 4 (sec) 

2.56 

2.0 

1.86 

18 

2.5 

2.49 

Steady-State 

Normalized 

CMG 


46.04 

rS0.6 

183 6 

591.1 

65 2 


Vanance 
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performance data with that obtained from a full set of PMAs and an output variance of 

0.1 (arc sec) 2 (see Table 4.4). 

4.3 RESULTS FOR DYNAMIC ACTUATORS 

In this section we add actuator dynamics to the SCOLE model and then re-solve the 
actuator selection problem solved above. The actuator dynamics are given in NASA’s 
original SCOLE document [15] to be first-order with a time constant of 0.1 seconds. 
That is, for each actuator (both PMA and CMG) the response of the actuator f,(t) to its 
input signal u,(t) is governed by 

fi(s)/ Ui (s) = [l/(.ls+l)] (4.7a) 

or 

t) = -10f,+ 10u, . (4.7b) 

There are several possibilities for the characteristics of the white noise associated 
with the actuators; white noise may be an input to (and thus be filtered by) the actuator, 
or it may be additive with the actuator output (thus unfiltered), or both. In this example 
four different actuator noise models are considered. Recalling that the non-dynamic 
actuators had additive white noise with intensity W, the following noise cases were 
studied for dynamic actuators: 

1. white actuator input noise of intensity W u = W; 

2. white actuator output noise of intensity W f = W; 

3. both input and output noise, each white and of intensity W; 
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4. both input and output noise, each white and of intensity W/2. 

The sensors are assumed non-dynamic (without phase lag). 

First we examine the effect of actuator dynamics on the maximal theoretically 
achievable accuracy. From [3], the lower bound y* on the steady-state variance of the i^ 
output is given by 

yi* = [CPC T ] ;i i= 1, .... ny . (4.8) 

The values of the lower bound on the LOS error for the fourth-order 2-D SCOLE model 
under study were calculated for the four different actuator noise cases listed above, as 
well as for the non-dynamic actuator model examined earlier. The results are shown 
below. 


Table 4.5: Maximal Accuracy for Different Actuator Noise Cases 


Noise Case 

No Dynamics 

£ 

II 

£ 

Wf=W 

WpW u =W 

Wf=W u =W/2 

Max. Acc. 

.086921 

.0691 

.086921 

.10072 

.07926 

(arc sec) 2 







From Table 4.5 it is clear that the addition of actuator dynamics along with 
retention of the white noise input to the system states (actuator output noise only, WpW) 
does not change the theoretical maximal accuracy; that is, y* is equal for the non- 
dynamic and the W f = W case. Also from the table, filtering of the actuator noise by 
passing it through finite actuator dynamics clearly improves the maximal accuracy. 
Finally it is noted that for case (3), W f = W u = W, the minimal LOS error is greater than 
that obtained by all but one of the controllers in the non-dynamic case. Thus for 
purposes of comparison only cases (1), (2) and (4) are studied in further detail. 
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For each of the three actuator noise cases a controller was designed (using OVA) 
which assigned the steady state LOS error variance to 0.1 (arc sec) 2 and minimized the 
sum of the normalized actuator variances among those actuators whose variances exceed 
their specifications (normalized variances greater than unity). For each final controller, 
the dynamic actuator effectiveness values for the PMAs are plotted in Figure 4.12 versus 
the actuators’ position along the 130 ft. flexible beam. 

For each of the noise cases the most effective actuator location is toward the 
reflector-end of the beam, with the highest effectiveness values corresponding to 
actuators located at the beam tip. Recalling the mode shape figures for the open loop 
flexible modes, the results in Figure 4.12 indicate that the PMAs are used by the optimal 
controller primarily for control of mode #1, which accounted for 46.6 percent of the open 
loop modal cost. It is interesting to compare Figure 4.12 with the plot of effectiveness 
values for non-dynamic actuators (Figure 4.8); note that the most effective non-dynamic 
actuators for the controller which achieved LOS error = 0.1 (arc sec) 2 were lo:ated near 
the center of the beam (70 ft from the shuttle). Hence, the optimal beam locations for 
PMAs in controllers which are achieving the same output performance are affected by 
the actuator dynamics. 


CONCLUSIONS 

The Closed-Loop Input/Output Cost Analysis (CIOCA) method of sensor and 
actuator selection (SAS) has been extended for application to systems with dynamic 
sensors and actuators — that is, systems in which the response of the sensors and 
actuators to their inputs is not instantaneous but governed by deterministic dynamics. 
The extended SAS method is applicable to systems in which the deterministic sensor and 
actuator dynamics are of arbitrary order. Application to simple numerical examples 
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demonstrates the utility of the SAS method. The examples also demonstrated that even 
uniform sensor dynamics can affect the optimal selection of sensors. Application of the 
actuator selection method in detail to NASA’s SCOLE space structure demonstrated that 
even uniform actuator dynamics can affect the optimal selection of actuators. 
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Proof of Theorem 1 


Let a = (nx p +nx a ). The composite system (2.7e) has a controllability matrix 

W c 6 R ax(ct ' nu) of the following form 

0 B 0 C a B a (ApB 0 C a B a +B 0 C a A a B a) 

B a A a B a A a B a 

• • • (A p a - 2 B 0 C a B a +A p a - 3 B 0 C a A a B a + • • • +B Q C a A a a - 2 B a ) 

(A,) 

Now noting that C a A a B a = is the ith Markov parameter for the system of actuator 
dynamics (2.7c,d), W cl may be rewritten 

W cl = 10 B 0 M 0 (A p B 0 M 0 +B 0 Mj) (A 2 B 0 M 0 +A p B 0 M,+B 0 M 2 ) • • • 

• • • (A“- 2 B 0 M 0 +Af- 3 B 0 M l + • • • +B Q M a _ 2 )] (A.2) 

The columns of W c span the controllable subspace of the composite system. Linear 
independence of all the rows in W c implies full controllability of the composite system. 
However, controllability of the original S’, stem states, x p , requires only that the columns 
of W cl span the state space for x 0 . This in turn will hold if and only if the matrix W cl has 
rank nx p . 

The proof of the Theorem begins with the proof that (A. 3) implies (A. 4): 

{ det M k * 0, Mj = 0, i =0, 1 k— 1 } (A.3) 

rank[W cl ] = n x (or range space of W cl has dimension n Xp ) < A.4i 

Note that the last block of W cl has the form 

w cl(a) = (A“- 2 B 0 M 0 +A«' 3 B 0 M 1 + • • • + B 0 M a _ 2 ) (A. 5) 

Now let k <, a-2 be the index of the first nonzero Markov parameter, M k . (In this case 


W c = 


w, 

w. 


cl 


c2 
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the first k+1 blocks of W cl are zero.) Next, use is made of two results from linear 
algebra (" R [K]" denotes "range space of K"), 

{detK*0} => {/f[JK]=/?[J]) (A. 6a) 

R [J+K] cR[J] + R [K] (A. 7b) 

(where cz ' means "is contained in") to demonstrate the following results which hold 
when M k is nonsingular 

R [B 0 ] = R [B 0 M k ] (A. 8) 

*[B 0 A p B Q ] = R [B 0 ] + R [A p B 0 ] 

= R [B 0 M k ] + R [A p B 0 M k ] 

= R [B u M k ] + R [A p B 0 M k + B 0 M k+l -B 0 M k+1 1 

c R [B 0 M k ] + R [A p B 0 M k + B 0 M k+1 ] + R [B 0 M k+1 ] 

= R [B 0 M k ] + R [A p B 0 M k + B c M k+1 ] 

R [B 0 A p B 0 ] c R [B 0 M k A p B 0 M k + B 0 M k+1 ) (A.9) 

Eqns (A. 8) and (A.9) lead by induction to the main result 

R [B 0 ApB 0 ••• A 0 a - k - 2 B 0 ]c/?tB 0 M k A p B 0 M k + B 0 M k+1 ••• 

A p a - k - 2 B 0 M k + • • • B G M a _ 2 ] 

that is, 

R [B 0 A p B 0 • • • A“ _k_2 B 0 ] c R [W cl l . (A. 10) 
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Condition (A.3) leads to (A. 10). Thus, given (A.3) together with (A p ,B 0 ) 
controllable, the columns of W cl are guaranteed to span the nx p -dimensional state space 
for Xp as long as 

a-k-2 > nx Q -l . 

that is, as long as 

k < nx a - 1 . (A.l 1) 

In fact, the index k of the first nonzero Markov parameter for the system (2.7) will 
always satisfy (A.ll). To show this, simply note that by observability of (A a ,B a ), the 
observability matrix W oa for (2.7) has full column rank: 

rank(W oa ) = nx a (A. 12) 

From (A. 12), 

{W oa B a = 0} => { B a = 0} => {Contradiction of (A a ,B a ) controllable) (A. 13) 

Thus, 

W oa B a = [M 0 T , Mj , .... * 0 (A. 14) 

and so the validity of (A.l 1) is guaranteed for (2.7) completing the proof that 

(f(A p ,B 0 ) controllable] & [M k * 0, = 0, i = 0, 1, ..., k— 1 ] } (A. 15) 

=> (x p controllable u) 

The usefulness of (A. 15) stems from the fact that by meeting certain mildly 
restrictive conditions the actuator dynamics can be guaranteed not to destroy 
controllability of the original system states x p , regardless of the pole/zero location for the 

plant. 
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It remains to prove the equivalence of the condition (A.3) and the requirements on 
the individual actuators’ numbers of poles and zeros. First, note that since each actuator 
is a single input, single output (SISO) system, then the Markov parameters M; for the 
lumped actuator dynamics (2.7) are diagonal matrices of the following form: 

M i = [ m n> m 2i* m 3i> -> ^nui ] (A. 16) 

where m^ is the (scalar) i^ Markov parameter for the actuator. Thus the condition 
(A.3) is met if and only if the index i of the first nonzero markov parameter is equal 
among all the actuators. 

The input/output transfer function for any n^ order SISO system has the form: 

T(s) = (c n _ 1 s n_1 +c n _ 2 s n_2 + • • • + c 0 )/(s n +d n _ 1 s n "4 • • • +d 0 ) (A. 17) 

The scalar Markov parameters m ; for the SISO system with transfer function (A.20) may 
be shown to be given by: 


m i = c n-2 - dn-l n o 


n 2- c n-3 -d n-2 n o~ <1 n-l n l 

n n = c o ~ ^1 ~^2 n l - * • • -dn-in^ 


(A. 18) 


From (A. 18), n { is the first nonzero Markov parameter for a system when the 
number of zeros in its transfer function is 


z = n - i - 1 (A. 19) 

Letting nj and Zj equal the number of poles and zeros for the actuator, respectively, 
(A. 19) yields the conclusion that 
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{det M k * 0, Mi = 0, i = 0, 1 , k-1 } «=* {(nj-Zj) = (n r z;) V ij.e (l,2,...,nu)) (A.20) 


Thus it is concluded that 

{ (nj-Zj) = (n j-zj) V i,j,e(l,2,...,nu)} => {x p is controllable u}. 


(A.21) 


25 7 



2S8 



N89-1347 1 

SCOLE 87/1 


OPTIMIZATION-BASED DESIGN OF CONTROL SYSTEMS 
FOR FLEXIBLE STRUCTURES 


E. Polak, T. E. Baker, T-L. Wuu and Y-P. Ham 


Department of Electrical Engineering and 
Computer Sciences 
University of California 
Berkeley, Ca. 94720 


Presented at 

The 4th Annual SCOLE Workshop 
Air Force Academy 
Colorado Springs, Colorado 
November 16, 1987 


This research was supported by the National Science Foundation grant ECS-85 17362, 


the Air Force Office Scientific Research grant 86-0116; and the Office of Naval 
Research contract N00014-86-K-0295. 


259 

PRECEDING page blank not filmbb 



SCOLE 87/2 


ABSTRACT 

The purpose of this presentation is to show that it is possible to use 
nonsmooth optimization algorithms to design both closed-loop finite 
dimensional compensators and open-loop optimal controls for flexible 
structures modeled by partial differential equations. 

An important feature of our approach is that it does not require modal 
decomposition and hence is immune to instabilities caused by spillover 
effects. Furthermore, it can be used to design control systems for struc- 
tures that are modeled by mixed systems of coupled ordinary and partial 
differential equations. 
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DESIGN OF STABILIZING FEEDBACK-SYSTEM COMPENSATORS 

The optimization-based design of finite dimensional compensators for 
systems modeled by coupled systems of ordinary and partial differential 
equations is made possible by a generalization of the following necessary 
and sufficient stability test for linear systems described by ordinary 
differential equations. 


THE DYNAMICAL SYSTEM 

Consider a parametrized, linear, time-invariant, interconnected, finite 
dimensional dynamical system, X(p), described by a set of state equations: 

x(t) = A(p)Xj(t) + B(p)u(t), 
y(t) = C(p)x(t) + D(p)u(t), 

We shall denote the characteristic polynomial of X(p) by X( S >P) and 
assume that the coefficients of x(s,p) are continuously differentiable in p. 
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S-STABILITY 

When, it is desired to ensure not only exponential stability of a closed 
loop system, but also to exercise some control over the location of its 
poles, it is convenient to make use of the following definition of S- 
stability. 

Definition (S -stability): Consider a linear, time-invariant, finite dimen- 

sional dynamical system Z of the form (1). Let S be an open unbounded 
subset of C which is symmetrical with respect to the real axis, and such 
that S c z> C + , where S c is the complement of S and C + is the closed right 
half of the complex plane. 

We say that the system L is S— stable if all the zeros of its charac- 
teristic polynomial are in S. ■ 
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A MODIFIED NYQUIST STABILITY CRITERION 

Theorem : Let S<= C be as specified in the Definition and let Be C be 

any simply connected set satisfying (0,0) 4 B. Suppose that 
D(s,q) e C[s] is a parametrized polynomial of degree N, whose 
coefficients depend on the parameter vector q e R" D in such a way that 
for every x(s) e Pn satisfying Z[x(s)]<=S, there exists a q x e R n ° such 
that 

(i) Z[D(s,q x )]cS, (2a) 

(ii) x(s)/D(s,q x ) e B, V s e dS. (2b) 

Then, given a polynomial x( s ) e Pn» ^[X( s )] c ^ ^ and only if there exists 
a q x e R° D such that (2a, b) hold. ■ 
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PROOF OF MODIFIED NYQUIST STABILITY CRITERION 

(= = =>) Suppose that Z[x(s)]c:S. Then, by assumption, there exists a 
q x e lR nD such that (2a), (2b) hold. 

(< = = = ) Next, suppose that (2a), (2b) hold. Then, because B is a sim- 
ply connected set which does not contain the origin, the locus traced out 
in the complex plane by %(s)/D(s,q x ), for s e 3S, does not encircle the 
origin. It now follows from (2a) and the Argument Principle that 
Z[x(s)]cS. ■ 

Comment : It is clear from the Theorem that an acceptable parametri- 

zation of the polynomial D(s,q) depends on the shape of the set S and the 
choice of the set B. A further requirement is imposed by semi-infinite 
optimization: the parametrization must be such that it is easy to ensure 
that the zeros of D(s,q) are in S. ■ 
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SCOLE 87/7 


OPTIMIZATION-BASED CONTROL SYSTEM DESIGN 

Id 



DESIGN VECTOR: x = [ x\ x 2 ,..., x 8 ] . 
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SCOLE 87/8 


DESIGN CRITERIA 

1. The feedback system must be exponentially stable. 

2. The system should have a good step input response. 

3. There should be little interaction between channels. 

4. Plant should not be saturated by command input effects. 

5. System should have high output disturbance rejection. 
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SCOLE 87/9 


MODIFIED NYQUIST STABILITY CONSTRAINT 



o i a » 

Initial Design 



STEP RESPONSES 



Channel 1 


Channel 2 



SCOLE 87/10 


CHANNEL INTERACTION CONSTRAINT 



Magnitude of H ! r i(jo),x) 



Magnitude of H y(j(l),x) 


COMMAND INPUT SATURATION CONSTRAINT 



Singular Values of H U:r (jti),x) 
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SCOLE 87/1 1 


OUTPUT DISTURBANCE SUPPRESSION CONSTRAINT 
Must accept some disturbance amplification outside operating 
bandwidth: 

o[H yd (jco,x) < 1.05, V co e [1,1000] 

COST: OUTPUT DISTURBANCE SUPPRESSION 

Suppress disturbance effects inside operating bandwidth: 

f(x) = max o[FL d (jcD,x) 
toe [0.001,1] * 



Frt^weacy (lad/Stc) 


Singular Values of H y( j(jco,x) 
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INTEGRATED STRUCTURE-CONTROL-SYSTEM DESIGN 


BEAM SHAPE AT TIME T 



♦ 

r 

it 


Vibrating Beam 
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DYNAMICS 

• GENERAL MODEL: Euler-Bernoulli Model, Kelvin-Voigt or 

Proportional Damping, Coupled Axial and Flexural Linear PDE’s. 

• Control Forces F*(t), Actuator Positions a 1 , Sensor Positions s'. 

• SIMPLIFIED MODEL: Decoupled Motion Formulation: 

ni 

mu tt (t , x) + CIu txxxx + EIu xxxx (t , x) = X b i( x - a i) F'(t) . 

i=l 

1 1 

y'(t) = J c s (C - s')u(t , Qd£ or y'(t) = J dj(C - s')u(t , Qd£ . 

0 0 

BOUNDARY CONDITIONS 

u(t , 0) = 0 , u x (t , 0) = 0, Ju ttx (t , 1) + CIu txx (t , 1) + EIu xx (t , 1) = 0 , 
Mu„(t , 1) - CIu txxx (t , 1) - EIu xxx (t , 1) = 0. 
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SCOLE 87/14 


DESIGN CRITERIA 

1. The feedback system must be exponentially stable. 

2. Control system compensator should be finite dimensional. 

3. Actuators should not be saturated by command input effects. 

4. System should have high mechanical disturbance rejection. 

5. Average power use should be low. 

6. Structure weight should be low. 

7. Structure should remain in elastic range. 
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SCOLE 87/15 


DESIGN VARIABLES 

• CONTROL SYSTEM COMPENSATOR 

(i) Coefficients of compensator differential equation. 

• STRUCTURE 

(i) Positions of actuators and sensors. 

(ii) Parameters of damping devices. 

(iii) Parameters of composite materials. 

(iv) Parameters determining shape of structure. 
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SCOLE 87/16 


PRELIMINARY RESULTS 

1. The control system can be stabilized using a finite dimensional 
proportional-plus-integral controller which ensures good distur- 
bance rejection. The use of our modified Nyquist stability cri- 
terion in the design of a stabilizing controller requires only 
evaluations of the system frequency response. Since the frequency 
response at a given frequency can be computed in some cases by 
formula and in the more general cases by solving two-point linear 
boundary value problems, there is no need for modal decomposition 
and hence there are no spillover effects. As in the finite dimensional 
case, time and frequency domain constraints can be treated simul- 
taneously and, in an integrated design approach structural parameters 
and constraints can also be introduced into the optimization problem. 

2. If a sequential design approach is used, an infinite dimensional com- 
pensator can be designed using an H°° frequency domain constraint 
formulation which results in a convex optimization problem and 
automatically ensures exponential stability with stability margin. 

3. An infinite dimensional controller designed as above can be 
approximated by a finite dimensional controller without spillover 
effects. 

4. A special semi-infinite optimization algorithm has been developed 
which is highly effective for design with H°° frequency domain 
design constraints. 
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SCOLE 87/18 


THE DYNAMICAL SYSTEM 


Hollow aluminum tube: one meter long, 2.0 cm diameter, 1.6 mm 
thick. Attached mass weighs 1 kg. We assume that motor torque u(t) can 
be directly controlled. 


Standard Euler-Bemoulli tube equations with Kelvin-Voigt visco- 
elastic damping: 


mw tt (t,x) + CIw txxxx (t,x) + EIw xxxx (t,x) — mQ 2 (t)w(t,x) 

= - mu(t)x, x e [0,1] 


(la) 


with boundary conditions: 

w(t,0) = 0, w x (t,0) = 0, CIvv txx (t,l) + EIw xx (t,l) = 0. (lb) 

M(Q 2 (t)w(t, 1) - w tt (t, 1) - u(t)) + CIvv txxx (t, 1) + EI\v xxx (t, 1) = Qic) 

where w(t,x) is displacement of tube from shadow tube (which remains 
undeformed during the motion), u(t) is motor torque, and Q(t) rad/sed is 
angular velocity. Above: m = .2815 kg/m, C = 6.89x1 0 7 pascals/sec., 
E = 6.89x1 0 9 pascals, I = 1.005 x 10 -8 m 4 , The tube is very lightly 
damped (0.1 per cent ). 
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SCOLE 87/19 


THREE OPTIMAL SLEWING PROBLEMS 


Minimize the time required to rotate the tube 45 degrees, from rest to 
rest, subject to the torque not exceeding 5 newton-meters. 


Minimize the total energy required to rotate the tube 45 degrees, from 


rest to rest, subject to the torque not exceeding 5 newton-meters and 
the maneuver time not exceeding a given bound. 


Minimize the time required to rotate the tube 45 degrees, from rest to 
rest, subject to the torque not exceeding 5 newton-meters and an 
upper bound on the potential energy due to deformation of the tube 
throughout the entire maneuver. 
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SCOLE 87/20 


THE DYNAMICAL SYSTEM 

MATHEMATICAL FORMULATION OF THE THREE PROBLEMS 

• To avoid technical problems associated with variable intervals and 
problems due to discretization, augment dynamics by one state variable 
and introduce scale factor T > 0 so that problem becomes defined on nor- 
malized time interval [0, 1], with T also equal to final time. 

• Tube is at rest when the total energy = energy due to rigid body 
motion + energy due to vibration and deformation = 0. 

(i) To ensure a slewing motion of 45°, we define 

g’(u,T) 4 (e-n/4) 2 (2) 

(ii) Rigid body energy at final time is proportional to the square of the 
angular velocity. 

g 2 (u,T) 4 Q(T) 2 . (3) 
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(iii) Kinetic energy due to tube vibration at normalized time 1 is 

1 

g 3 (t,u) 4 w t (l ,x) 2 dx. (4) 

2 0 

(iv) Potential energy due to tube deformation at normalized time 1 is 

g 4 (l,u) 4 w„(l ,x) 2 dx. (5) 

A 0 

• Potential energy due to deformation of the tube at normalized time t: 

P(t,u) 4 -S-f w xx (t,x) 2 dx. (6) 

2 0 

(v) To limit tube deformation for all t e [0, 1] we define 

T 

g 5 (u,T) ^ J [max{ P(t,u) - f(t) , 0 }] 2 (7 

0 

(vi) To ensure slewing time does not exceed T f seconds, we define 

g 6 (u , T) = T - T f . (8 
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SCOLE 87/22 


FINAL PROBLEM FORM 

P,: Te iRr! n u e G {g 0 (U ’ T)lg,(U ’ T) " e - 0 ’j e n. 2,3,4} }, 
where g°(u,T) 4 T, 1R + 4 (ye Rly>0}and 

G 4 { u e L„[0, 1] I lu(t)l < 5, te [0,1] }. 

Pi TelRr; n ue G (g0(U ’ T)l8j(U>T) “ e - 0 ’j e U. 2.3,4,6 1 }, 

1 

where g°(u,T) 4 J ||u(t)|| 2 dt. 

0 


min 


{ g°(u,T) I g j (u , T) — e < 0 , j g { 1,2, 3, 4, 5 } }, 


T e IR+, u e G 


where g°(u , T) 4 T. 


• All gj are continuously differentiable in L„[0, 1], 
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SCOLE 87/23 


THE DYNAMICAL SYSTEM 
COMPUTATIONAL RESULTS 

IMPLEMENTATION. Because we cannot solve the system PDEs 
exactly, we cannot evaluate gi(u,T) or Vg-*(u,T) exactly. Furthermore, 
since u is an infinite dimensional design vector, it can only be entered 
into a computer in discretized form. We use an implementable algorithm 
which adjusts integration precision and control discretization adaptively. 
To discretize the PDE in space, we use the finite element method. Since 
the PDE is fourth order in space, it is necessary to use elements of at least 
second order. We have chosen Hermite splines as basis elements. The 
input u € G is discretized in time and Newmark’s method is applied to 
evaluate the resulting system of ordinary differential equations. 

LINEARIZATION. The results presented are for the case in which 
the Q 2 (t) terms are neglected in equation (1). Similar results have been 
obtained by performing experiments when the ft 2 (t) terms are included. 
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SCOLE 87/24 


OPTIMAL CONTROL FOR MINIMUM-TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY 


Control 
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SCOLE 87/25 


TIP DISPLACEMENT FOR MINIMUM-TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY 


Meters 
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SCOLE 87/26 


DEVIATION FROM SHADOW BEAM FOR MINIMUM-TIME 
PROBLEM WITH TORQUE CONSTRAINTS ONLY 


w(t,x) (meters) 
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SCOLE 87/27 


CONSTRAINT VIOLATION FOR MINIMUM-TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY: 
DISCRETIZATION EFFECTS 


PSI 



2QS 



SCOLE 87/: 


OPTIMAL TORQUE 

FOR MINIMUM-CONTROL-ENERGY PROBLEM 
WITH TORQUE CONSTRAINTS AND FINAL TIME < 0.8 SEC. 


Control 




SCOLE 87/29 


OPTIMAL TORQUE 

FOR MINIMUM-CONTROL-ENERGY PROBLEM 
WITH TORQUE CONSTRAINTS AND FINAL TIME < 1.0 SEC. 
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SCOLE 87/30 


POTENTIAL ENERGY FOR MINIMUM TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY 

Curve A is potential energy 
Parabola B is deformation constraint. 


Potential Energy 
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SCOLE 87/31 


POTENTIAL ENERGY FOR MINIMUM-TIME PROBLEM 
WITH TORQUE AND POTENTIAL ENERGY CONSTRAINTS 


Curve A is potential energy 
Parabola B is deformation constraint. 


Potential Energy 



239 



SCOLE 87/32 


OPTIMAL CONTROL FOR MINIMUM-TIME PROBLEM 
WITH TORQUE AND POTENTIAL ENERGY CONSTRAINTS 


Note: The optimal final time is 0.8177 seconds, an increase of only 3.7 

percent over the solution of Pj. 


Control 
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EFFECT OF ACTUATOR DYNAMICS ON CONTROL 
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Simulation of Two Aspects of Physical Limitations 



I 



DISCOS SIMULATION: 

BODIES CONNECTED BY HINGES 

FINITE ELEMENT MODEL OF BEAM PROVIDED BY NASTRAN 


REFLECTOR (RIGID), BODY 3 



HINGE 1 Z 


HINGE 1 CONNECTS SHUTTLE 
BODY TO INERTIA SPACE 


Fig. 1— Scole Configuration 
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TWO PARTS OF SIMULATION 


OX) 

c 


<L> 

c/5 


<D 

C/5 

I— 

<L> 

> 

0) 

u. 

C 

<D 




*-> 

o 

*o 




o 


r> 

w 

X 


C/5 

*o 

C/5 

2 


c 

o 

T3 

c 



o 

o 

c 

X 


C-> 

C/5 

o 

<L> 

C/5 


in 

C/5 

c 

t— 

vo 

O 

o 

o 

•> in 

m 

« 


£ 

<D 

C/5 

"S c/5 ^ CO 
-n t, ^ - 

c £ C 

CO to £ 


U, t-H 

” O 

U_H 

— ■ c 


£ 2 
£ jc 

5 ~ 

u X) 

9r^ 

O 4-1 

o 


c 

o 

<D 

O 

»- 

o 


o _ 


o 

° c 

OX) ° 

£ C/5 

CO 


. © ^ 
COX) 

<L> o ^ 

CL, --T _ 

n ° ° 

W in 


<L> 
•*— > 

0) 


CL 

E 

o 

x a> o 

g»l * 

si <2 

1-0 4-* C4> 


u- 

0) 

w 

03 

T3 

C 

03 

£ 

« 

Oil 

c 

• — S 

s- 

S3 


C4 

a 


5 2 
S* 
c 3> 

E 2 

cO 

£ o 

4-1 2 

° -S- 

_ CO 

£ P 
O Q - 1 

• P«4 

w 

c P 

w> g 

(l> c 


(N 


291 


vernier thrusters (60 ft-lb torque) on shuttle and reflector 
2 sets of 2 axes proof-mass actuators on beam 
each set of actuators has 10 lb force in both M x" and "y" 
maximum of 1 ft travel distance of proof-mass 



PROCEDURES OF THIS ANALYSIS 
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Thruster forces on reflector 
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J[x T (s)Qx(s) + uT(s)Ru(s)]ds 





a> 


3 

c 



e 

o 

CJ 


s- 

o 

Is 

IS 

3D 

a> 

O' 


03 

C 

.2 

-4-J 

o 

G 

3 


C/2 

o 

CD 

Ctf) 

C 

N 

• *—H 

s 

G 


CD 

CD 

<D fc 

O £ 

c ** 


c3 


u 

O 

4—* 

03 


*h 

€ = 

S W) 

p <D 
^ Vh 

f-M 

2 o 

03 

H3 ^ 

«>.2 
<u w 
Ui 03 
■*— > 

o *5 

N .S 


C/2 

CD 

> 


■S£ 

a 

*3 O 


CD 

<D 

jd* 

O 


P4 






On 

H 

(3 

3 

Cu 


i 


CD 

S3 

u 

X) 

W) 

r— H 

03 


o 

s-< 

C 

o 


CD 

<D 

> 

'o 

c/3 


< 

G< 

+ 

P-i 

<! 

+ 

O' 

ii 

o 


vO 

i 

o 


s- 


jg 

4-* 

. 

£ 



O 

C/5 





297 


x(t) can include a reduced set of modes 



ASSESS EFFECTS OF 1 FT TRAVEL LIMIT 
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Significantly less damping in run #3 than in the other runs. 
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ROLL OF RNTENNR RELRTIVE TO SHUTTLE VERSUS TIME 
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COMPARISON OF CONTROL EFFECTIVENESS 
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ROLL or ANTENNA RELATIVE TO SHUTTLE VERSUS TIME 


O 



S33y03d 

30& 


I 


CONCLUSION 
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NEED MORE ROBUST REGULATOR to do better 
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COMBINED PROBLEM OF SLEW MANEUVER CONTROL 
AND VIBRATION SUPPRESSION 

Y. P. Kakad 

Dept, of Electrical Engineering 
University of North Carolina at Charlotte 
Charlotte. NC 28223 


ABSTRACT 

In this paper, the combined problem of slew maneuver control 
and vibration suppression of NASA Spacecraft Control Laboratory 
Experiment (SCOLE) is considered. The coupling between the rigid 
body modes and the flexible modes together with the effect of the 
control forces on the flexible antenna is discussed. The nonlinearities 
in the equations are studied in terms of slew maneuver angular velo- 
cities. 


INTRODUCTION 

In this paper, the analytics for the combined problem of slew maneuver and 
vibration suppression are developed. It is assumed that the slew maneuver is per- 
formed by applying moments on the rigid shuttle and the vibration suppression is 
achieved by means of forces on the flexible antenna and the reflector. The slew 
maneuver is considered to be an arbitrary maneuver about any given axis [16]. The 
effect of slew maneuver angular velocity on flexible modes is studied by examining 
the spectral norm of the matrix term associated with the coupling between the 
rigid-body modes and the flexible modes. Also, the kinematic nonlinearities are 
further analyzed in terms of the matrix spectral norm variation of the correspond- 
ing term with respect to slew maneuver angular velocity. 
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ANALYTICS 


The slew maneuver is defined as 
X_ — Axis about which the slew maneuver is performed. 

£ - The slew Angle 

to - The angular velocity of the orbiter in the inertial frame. 
The four Euler parameters can be defined as 


<u = X x sin|- 
£ 2 = X 2 sln|- 

£3 = X 3 sin|- 

£ 

£4 = cos-j- 
2 


( 1 ) 


The four Euler parameters can be related to the angular velocity components 
of the rigid assembly as 




*1 *4 — e 3 c 2 


0 

€2 


£ 2 £3 €4 -«1 


0>i 

*3 


£3 -£ 2 €i €4 


0) 2 

*4 


*4 “€i -€ 2 -€3 




The slewing maneuver can be given in terms of the following equations [16] 


/ 0 to + = G(t )+yy 2 (<o) 

A 1 to + A 3£. + Bg_ + Kg_ = Q_{t ) 


(3) 

(4) 


where, 

G (r ) is the net moment applied about the mass center of the orbiter and is 
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given by the following equation (figs. 1 & 2) 

G (r ) = G^Cr ) + (r. + o)xF 2 (5) 

Also, Q{t ) represents the generalized force vector which is given by the following 
equation 


aco = 


ru 

Z ( Qjxib ) + Qjyfi » + Qx t + Q yi + Q* 
Z C Q,*^ ) + Qjyfi )) + Q X2 + Qy 2 + Q i , 2 


Z c ) + Q, y , (f )) + Q* + Q yi + 


where, the generalized force components are given as 


( 6 ) 


Qjx, — f Fjx (z >t )8(z —Zj )<t> xi (z )dz 
o 

L 

Qjy, = f Fjy (Z ,t )8(z — Zj )0 y j ( Z ~)dz 


( 7 ) 

( 8 ) 


and 


Q ,^(0 = 0 ( 9 ) 

Here, Fj x (z ,t ) is the x component of the concentrated force applied at location j 
on the flexible antenna and is the y component of that force. 

Also, 
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Qyi (t ) — F 2y ( t )0 yi (Z, ) (10) 

) = A/ ^,(? (Z. ) 

Here, F 2 is the force applied at the reflector C. G. 

Thus, 

M ,f,(t ) = F 2 x r y + F Zy r x + M (11) 

The location of reflector C. G. is given by coordinates (r x ,r y ) and M 2 ^ 
represents the external moment applied at the reflector C. G. Also, the nonlineari- 
ties N _2 can be expressed in terms of pure rigid body kinematic nonlinearity and 
the nonlinear coupling term between the rigid-body modes and the flexible modes. 

N.2 = A 4 ( cj, 6 _ ) + A 5 ( ) g_ ( 12 ) 

(a) Slew Maneuver 

If only a slew maneuver is to be considered, then Q{t ) = £ and F 2 =0, and 
only moments are applied at the orbiter C. G. However, the angular velocity vector 
to, is nonzero during the maneuver and the flexible modes will be excited. This 
effect of coupling between the rigid-body modes and flexible modes can be obtained 
by evaluating A 5 which depends on the angular velocity vector. In figure 3, using 
the matrix spectral norm as a measure, the coupling effect is studied as a function 
of slew angular velocity. The first ten flexible modes are considered for this 
analysis. The kinematic nonlinearity is also obtained in terms of matrix spectral 
norm as a function of to. This analysis can be utilized in the linearization of the 
slew maneuver dynamical equations. An example of this is shown in figure 4 
which is a single plane slew maneuver. In this case, it is almost a linear relation- 
ship in terms of a single angular velocity component. 
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(b) Slew Maneuver Control and Vibration Suppression 

If it is desired to design control systems for the simultaneous task of slew 
maneuver control and vibration suppression, then equations (3)-(ll) should be 
used. It can be seen that vibration control forces also affect the slew maneuver 

3 . . , - ■ j i 1 * * t i. 

vaj Liuuu^u tuuuui iiiuiij.cj.ji c cuajjxiu^ lca iu&, 

Thus, these equations would suggest that in order to achieve control efficiency 
and to minimize the line of sight error in minimum time, it may be necessary to 
synthesize control systems for the combined problem of slew maneuver and vibra- 
tion suppression. 
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Figure 1- Position Vectors in Inertial Frame 
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Figure 2- Vectors in Body-fixed Frame 
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ABSTRACT 

The design of a robust compensator is considered 
for the SCOLE configuration using a frequency - 
response shaping technique based on the LQG/LTR 
algorithm. Results indicate that a tenth-order 
compensator can be used to meet stability- 
performance- robustness conditions for a 26th-order 
SCOLE model without destabilizing spillover 
effects. Since the SCOLE configuration is 
representative of many proposed spaceflight 
experiments, the results and design techniques 
employed potentially should be applicable tf» a wide 
range of large space structure control problems. 


Introduction 

Large space structures (LSS) have many properties 
that make them difficult to analyze and control 
[1]. They are mathematically modeled by computa- 
tionally difficult partial differential equations 
or high-order, lumped, ordinary differential 
equations obtained through finite element methods. 
LSS have many low and closely spaced resonant 
frequencies, a number of which typically fall 
within the controller bandwidth. In LSS, vibra- 
tional issues must be treated as a first-order 
effect; it is this characteristic of the LSS 
control problem that most distinguishes it from 
spacecraft control problems of the past. Addition- 
ally, inherent damping is low and/or improperly 
modeled. Coupled with stringent operational 
requirements for orientation, shape control, and 
vibration suppression, these properties present an 
unconventional and unresolved control design 
problem to the system analyst. 

A fundamental issue to be dealt with in any LSS 
control problem comes from the large amount of 
modeling error occurring In finite element models 
of such structures. In general, inaccuracy of 
modal data, such as elastic frequencies and mode 
shapes used to form coefficient matrices of the 
dynamic models, increases with increasing modal 
frequency. Hence, a frequency -dependent constraint 
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is inherently imposed on the design process i\ that 
stabilization and performance requirements must Oe 
met without allowing the input control energy to 
"spill-over" and excite and destabilize the lightly 
damped, poorly modeled high-frequency dynamics. 

At NASA's Langley Research Center, a LSS config- 
uration known as the Spacecraft Control Laboratory 
Experiment (SCOLE) was conceived for the purpose of 
evaluating and comparing large space structure 
control and identification concepts [2], The SCOLE 
configuration (shown schematically in Figure 1) 
consists of a 130- foot flexible beam anchored at 
one end to the cargo bay of the space shuttle with 
an antenna reflector connected to the opposite end. 

The center of mass of the reflector is offset from 
the attachment point. The SCOLE configuration is 
representative of many proposed space flight exper- 
iments and space-based antenna systems. Control 
inputs are available from torque actuators located 
on the orblter and force actuators at the reflector 
center. Attitude sensors are located at the 
reflector center. A typical SCOLE control task 
is to slew or change the line -of - sight of the 
antenna rapidly and damp any induced structural 
vibrations to the degree required for the precise 
pointing of the antenna. 

In this paper we consider a SCOLE large -angle 
slewing maneuver to have been completed and attack 
the problem of designing a model -based compensator 
to attenuate residual structural vibrational motion 
and antenna line-of- sight error. The SCOLE mathe- 
matical model is first discussed followed by des- 
criptions of the design objectives and the compen- 
sator design approach. Finally, results from the 
application of the design methodology to the SCOLE 
problem are presented. 

Mathematical Model 

The basic distributed-parameter mathematical model 
of the SCOLE configuration is described in [2], 
while nonlinear and linear ordinary differential 
equation models are found in [3] and [4], respec- 
tively. A linear finite-element model consisting 
of three rigid rotational modes and the first ten 
structural elastic modes is used in this study. A 
state- space realization of the modal model has the 
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form 


where 


and 


for (i-1, ... ,10) . Equation (4) describes the rigid 
body contribution and equations (5) and (6) des- 
cribe the elastic contribution for ten vibrational 
modes of frequencies (i-1 10). A uni- 
form damping ratio of f • 0.003, (i - 1,..., 

10), is assumed. The eigenvalues of A are given 

Et 

mathematically by 

* J (?) 

and are shown in Table I. 

Five control inputs are generated using three 
torque actuators (one per X,Y,Z axis) on the orbi- 
ter and two force actuators (X and Y directions in 
Figure 1) at the reflector center. Three attitude 
sensors (one per axis) are located at the reflector 
center. Sensor and actuator dynamics were not 
included in this study. Rigid-body inertias, 
mode shape and slope data from the finite element 
analysis combine to define the control effec- 
tiveness matrix B_ and output response matrix C_. 

r r 


*F " A F X F + 1 
y F “ C F X F 


A R“ 


0 *3 
0 0 


A e - diag< Ag, 


l F U 

(1) 


(2) 

a e > 

(3) 


(6) 

6x6 


-*?> 

(5) 


a e“ 


■"i - 2f i“i 


2x2 


( 6 ) 


Design Oblectlves 

The design objectives of this study are to produce 
a multivariable, model -based, feedback compensator 
operating on attitude sensor data which will gener- 
ate force and torque Inputs to stabilize the rigid 
body modes; enhance the stability of lightly 
damped, low-frequency modes without destroying the 
stability of higher -frequency modes; meet pre- 
scribed closed- loop performance (bandwidth) speci- 
fications; and possess some degree of stability 
robustness to unmodeled dynamics. Since a low- 
order controller is sought, it was decided to 
employ full -state controller design with a reduced - 
order plant model. The full -order model is 
reserved for evaluation purposes. Order reduction 
for the design plant was performed using modal 
truncation. Past studies ([6], [7]) have indicated 
that a 0.1 rad/sec closed- loop performance band- 
width is sufficient to maintain antenna pointing 
control, and a design model composed of the rigid 
body plus the first three elastic modes in Table I 
is adequate to achieve this bandwidth. Higher 
bandwidths will typically require the addition of 
extra elastic modes to the design model. Denoting 
the 12th-order design model transfer matrix by 
Gp(s), for a unity-gain feedback compensator with 

transfer matrix G (s), multivariable bandwidth will 
c 

be defined as the frequency below which the 
smallest singular value of the closed- loop response 
matrix 

G CL (S) ■ G p< s > G c < s >( 1 + G p<s> G c (s> l' 1 (8) 

remains above unity for s — Jw. In our case, we 
seek a compensator such that 

- [ G Cl(J w ) 1 ^ for 0 5 w 5 0 .1 (9) 

From the block diagonal structure of A in (1), the 

F 

transfer matrix, G^s), of the 26th-order system 
may now be written as 


Analysis of (1) and (2) verifies that the system is 
completely controllable and observable. Attempts 
to reduce the number of control variables to the 
number of outputs retained controllability and 
observability but, In each three- control input 
combination, introduced lightly- damped, low- 
frequency transmission zeros [3] into the model. 
Since the presence of such transmission zeros has 
been demonstrated to reduce system performance in 
large space structure controller designs [6] , the 
compensator was designed with the original five 
inputs and three outputs. However, in order to 
avoid numerical ill -condi tioning brought about by 
the different physical characteristics of forces 
and torques, the inputs were scaled so that the 

frequency response of the largest (o(jw)) and 
smallest (a(ju>)) singular values of the transfer 

matrix of (1) and (2), denoted by G jr (jco) , were 

nearly equal at low frequencies (as shown in Figure 

2 ). 


G F (S) ~ G p (s) + AG(s) (10) 

where AG(s) represents the transfer matrix of the 
remaining 14th-order (residual) modal system. In 
this form, the dynamics represented by AG(s) can be 
interpreted as an "additive perturbation" to the 
Gp(s) system and used as an approximate represen- 
tation of unmodeled dynamics for use in stability 
robustness tests. Specifically, it is established 
in [8] that the unmodeled dynamics AG(s) will not 
destroy the closed- loop stability so long as 

o[ G c (s)[ I + C p (s)G c (s) l' 1 ) a [ AG ( s ) ] < 1 (11) 

for all s - Jw, o) real. Condition (11) can be 
enforced in the compensator design stage to ensure 
that closed-loop stability will be preserved for at 
least that class of unmodeled dynamics whose 

spectral norm lies below a[AG(jw)]. 
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The spillover effect on AG(s) due to the closed - 
loop compensation may be tested directly by 
applying a state-variable realization of G^Cs) to 

the full 26th- order model given by (1) and (2) and 
examining the eigenvalues of the composite system. 

A block diagram for the closed- loop configuration 
is shown in Figure 3. An approach for constructing 

G (s) to stabilize G <s) while satisfying (9) and 
c p 

(11) is presented in the next section. 


which, when used in an LQG fashion with the Kalman 
filter from Step 1 , asymptotically recovers the 
frequency response of the target loop gain over the 
low-frequency band. The resulting compensator is 
given by 


where 


G (s) - F( sl-k 
c 

T 

F - BP 

A T P + PA - PBB T P + qC T C - 0 


(17) 

(18) 
(19) 


Compensator Design Approach 

The compensator design approach employed to meet 
the foregoing design objectives can be viewed as a 
variation of the well-known Linear-Quadratic - 
Gaussian/Loop -Transfer -Recovery (LQG/LTR) algorithm 
([9), [10]). In the standard LQG/LTR approach, 
with the loop in Figure 3 broken at the output, a 
Kalman filter (G RF ) is designed to meet the 

complete set of stability -performance - robustness 
objectives. Thereafter, an optimal linear regu- 
lator is constructed such that the composite LQG 
compensator (G ) loop gain behavior asymptotically 

approaches (recovers) that of G^ F in the sense that 
C p (jo)G c (jcj) — G kf O> 

pointwise in w. Direct application of this LQG/LTR 
procedure to large space structures problems 
results in extremely conservative designs which 
cannot meet reasonable performance specifications 
[6]. However, the LQG/LTR structure still provides 
a viable approach for model-based controller 
synthesis when the standard procedure is modified 
in the following manner. 


Step 1 

Denote a state -variable realization of G p (s) by 

x - Ax + Bu (12) 

y - Cx (13) 

Select the design parameters L and p in the Kalman 
filter algorithm 

AQ + QA T + LL T - QC T CQ - 0 (14) 

H - — QC T (15) 

such that 

C KF (S) " C( Sl ‘ A } ^ (16) 

achieves a desired (target) loop gain for 
G (s)G (s) over some low-frequency band containing 
P c 

the design bandwidth. 

Step 2 

By successively increasing q > 0 in equation (19) 
(to follow), design an optimal linear regulator 


& - A-BF-HC (20) 

Step 3 

Attempt to adjust q in Step 2 until the desired 
bandwidth condition (condition (9)) is met. Also 
check stability robustness by (11). If an exces- 
sively high q (indicated by violation of (11)) is 
required to achieve the required bandwidth, turn 
down the Kalman filter gain (by increasing p in 
(14)) to "loosen" the target loop. In effect, this 
procedure reduces the target bandwidth until satis- 
faction of (11) is possible. The final design is 
accomplished by Iteratively adjusting the linear 
regulator and Kalman filter design parameters until 
an appropriate compromise is made between bandwidth 
and stability robustness. 

In large space structures applications, the 
inability (at Step 1) to meet loop gain magnitude 
over the desired bandwidth or (in Step 3) the 
production of too small a compromise bandwidth can 
often be overcome by the inclusion of additional 
flexible modes into the design model [6]. 

If, as in the SCOLE application to follow, an 
order - reduction study is performed on the resulting 
compensator, the complete set of stability- 
performance- robustness conditions needs to be 
re-evaluated with the reduced- order compensator. 


SCOLE Application 

Figure 4 shows the frequency response of the 12th- 
order (LQG) compensator, G c , resulting from an 

application of the foregoing procedure to the 12th- 
order SCOLE design model, G p . The figure indicates 

a well-behaved lead- lag structure with a 20 
db/decade roll-off. Eigenvalues of the corres- 
ponding A-HC, A-BF, and A matrices are given in 
Table II. The frequency response of G^ in equa- 
tion (8) with the 26th-order evaluation model used 
in place of G p is shown In Figure 5. Figure 5 

demonstrates the satisfaction of the 0.1 rad/sec 
bandwidth requirement and a 60 db/decade roll-off. 
Eigenvalues of the composite system resulting from 
the LQG compensator applied to the evaluation model 
are given in the first column of Table III. The 
data Indicate that the compensator was designed to 
concentrate on stabilizing the rigid body modes 
(first three entries of the A-HC block) without 
disturbing the stability of the three elastic modes 
(last three entries of the A-HC block) of the 


32 4 


I 



4 of 6 


design model. Table III also shows that there Is 
Insignificant spillover into AG. The stability 
robustness test (11) shown in Figure 6 shows more 
than 10 db robustness margin. 

The possibility of a reduced-order compensator 
(ROC) satisfying the design conditions was also 
investigated. The methods of balanced realization 
[11], Hankel-norm reduction [12] and a method based 
on stable factorization [13] were employed. The 
ratio of largest to smallest Hankel singular values 
was 0.016 so little reduction based on nearly 
uncontrollable or unobservable compensator modes 
was expected. All of the methods gave similar 
results. In each order reduction method, only a 
lOth-order ROC would stabilize the design model. 

The stable factorization results were judged to be 
(slightly) better and will be discussed herein. 

j^able II shows the eigenvalues oj^ the corresponding 
A matrix in the ROC, denoted by A^^,. Figure 7 

shows a frequency response of the ROC. A compar- 
ison of Figures 4 and 7 shows that the only 
difference between the LQG and ROC frequency 
response plots is the removal of the dip in a at 

the frequency of the third elastic mode. The 
importance of this characteristic can be seen from 
Figure 8 where the robustness condition (11) is 
evaluated using the ROC in place of G^ . Figure 8 

indicates that an effect of the-order reduction is 
a reduction of stability margin at the frequency of 
the third elastic mode. A more positive effect 
from the reduced order compensation is seen in the 
second column of Table III where the eigenvalues of 
the ROC applied to the evaluation model are presen- 
ted. The stability of the rigid-body modes from 
the LQG compensator is preserved with the auxiliary 
effect of adding stability to the first elastic 
mode. The net effect of the ROC is to enhance the 
stabilizing effect of the LQG compensator at the 
expense of a reduction of stability robustness 
margin. 


Concluding Remarks 


stability robustness margins. Since the SCOLE 
configuration Is representative of many proposed 
spaceflight experiments, the results and design 
techniques employed should potentially be 
applicable to a wide range of large space structure 
control problems . 
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TABLE I 

EIGENVALUES OF A £ 


Mode 

★ 

Eigenvalue 

1 

(-0.00524, 1.747) 

2 

(-0.00591, 1.970) 

3 

(-0.0513, 5.108) 

4 

(-0.0224, 7.449) 

5 

(-0.0387, 12.903) 

6 

(-0.0898, 29.925) 

7 

(-0.104, 34.657) 

8 

(-0.232, 77.165) 

9 

(-0.243, 80.993) 

10 

(-0.446, 148.780) 


# TABLE II 

EIGENVALUE ANALYSIS OF COMPENSATORS 


A-HC 

A-BF 

A 

Koc 

(-0.126, 0.126) 
(-0.126, 0.126) 

(-0.0563, 

(-0.0834, 

0.101) 

0.0837) 

(-0.420, 

(-1.038, 

5.125) 

2.027) 

| (-0.776. 

1.876) 

(-0.126, 0.126) 

(-0.0959, 

0.104) 

(-0.203, 

1.974) 

(-0.204, 

1.973) 

(-0.00524, 1.747) 

(-0.178, 

1.976) 

(-0.0489, 

0.184) 

(-0.0500, 

0.185) 

(-0.00591, 1.970) 

(-0.905, 

1.956) 

(-0.197, 

0.237) 

(-0.197, 

0.237) 

(-0.0153, 5.108) 

(-0.418, 

5.130) 

(-0.209, 

0.210) 

(-0.209, 

0.210) 


TABLE III 

EIGENVALUES OF FULL-ORDER SYSTEM FORCED 
BY LQG AND REDUCED -ORDER COMPENSATORS 


A-BF 


A-HC 


AG 


LQG 


(-0.0569, 0.101) 
(-0.0834, 0.0837) 
(-0.0971, 0.103) 
(-0.179, 1.976) 
(-0.905, 1.975) 
(-0.418, 5.130) 

(-0.122, 0.128) 
(-0.126, 0.126) 
(-0.126, 0.125) 
(-0.00524, 1.747) 
(-0.00591, 1.970) 
(-0.0153, 5.108) 

(-0.0243, 7.449) 
(-0.0383, 12.903) 
(-0.0898, 29.926) 
(-0.104, 34.657) 

(-0.232, 77.165) 

(-0.243, 80.993) 

(-0.446, 148.780) 


ROC 


(-0.0550, 0.0998) 
(-0.0834, 0.0837) 
(-0.0960, 0.103) 
(-0.181, 1.975) 

(-0.634, 1.823) 


(-0.122, 0.129) 
(-0.126, 0.126) 
(-0.128, 0.129) 
(-0.152, 1.752) 
(-0.00603, 1.970) 
(-0.0137, 5.109) 

(-0.0242, 7.449) 
(-0.0383, 12.903) 
(-0.0898, 29.259) 
(-0.104, 34.657) 

(-0.231, 77.165) 

(-0.243, 80.992) 

(-0.446, 148.780) 


Eigenvalues presented in 
(Real, + Imaginary) format. 
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The abjective of this work is to develop a practical sensor analytic 
redundancy management scheme for flexible spacecraft and to demonstrate 
it using the SCOLE experimental apparatus. The particular schere to be 
used is taken from previous work on the Grid apparatus by Williams and 


Reference: 

Williams, J. P. and R. C. Montgomery: Failure Detection and Accormoda- 

ion in ructural Dynamics Systems using Analytic Redundancy. 24th I EEE 
CDC, December 11-13, 1985. 


OBJECTIVE OF WORK 


DEVELOP k TEST A PRACTICAL SENSOR ARM 
SCHEME USING SCOLE 


APPROACH 

USE SCHEME PREVIOUSLY DEVELOPED FOR THE 
GRID BY WILLIAMS AND MONTGOMERY 
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OUTLINE 


The presentation is organized as follows: First, the scheme used by 

Williams and Montgomery is summarized. The scheme is based on a LOB 
design which is next described. Experimental results taken from the 
SCOLE apparatus on the performance of the Kalman filter of the LOG are 
presented and finally plans for completion of the work are given. 


OUTLINE 

SUMMARIZE THE GRID SCHEME OF WILLIAMS AND 
MONTGOMERY 

DESCRIPTION OF THE LQG DESIGN FOR THE 
SCHEME 

RESULTS FROM THE SCOLE LAB EXPERIMENT 
PLANS FOR COMPLETION OF THE WORK 




GRID ARM SCHEME - SUTTWY 


The approach of Williams and Montgomery was to use a single active 
steady state Kalman filter which is designed for the estimated failure 
state in effect. Uhder the no— failure case the sensor residuals of this 
filter should be white with zero mean. The zero-mean character of the 
estimated residuals is monitored using Wald's sequential probability 
ratio test (SPRT). SPRT is a binary test to determine if a statistical 
variable is zero-mean or has a mean, m. As data samples are gathered 
a decision variable is monitored. It is initialized at zero and is 
sequentially modified by the data samples. If it crosses either of two 
decision thresholds a decision is made. One threshold corresponds to 
the zero-mean decision while the other is for ths m— msan decision. A 
SPRT is run an each residual . If a decision of ze rcr-mean is made the 
SPRT is reinitialized and run again. If a m - mean decision is made a 
failure is declared. In event of a declared failure the failure 
signature of the sensors in the residuals are exanined to determine the 
failure state. A new LOG design for that failure state then replaces 
the current active design. 


GRID ARM SCHEME - SUMMARY 


USE SINGLE, ON-LINE, KALMAN FILTER 

USE SPRT TO CHECK THE ZERO MEAN CHARACTER 
OF THE ESTIMATED MEASUREMENT ERROR 

IF FAILURE IS DETECTED, ISOLATE USING 
FAILURE SIGNATURE IN THE ESTIMATED 
MEASUREMENT ERRORS 
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SCOLE AFtt LOG DESIGN 


The basis of the A FJ1 to be used is the LOG. Therefore the first 
order of business is to develop a sui table LOG design wherein the 
modelling errors do not defeat the zero-mean character of the residuals. 
Most of the rest of the presentation concerns this design and its 
performance. For the design model we have used a 5-mode, modal model of 
STTI E with the S CG L E platform fixed. Thus, there are no rigid body 
modes. Also the 5 modes selected are the five lowest frequency modes. 
Reaction jets are included in the filter but not in the regulator. The 
torque wheels on the other hand are used in both the filter and the 
regulator. 


SCOLE ARM LQG DESIGN 

DESIGN MODEL — MODAL MODEL 
FIXED SCOLE PLATFORM MODEL 
NO RIDID BODY MODES 
5 LOWEST FREQUENCY VIBRATION MODES 

JETS INCLUDED IN FILTER, NOT IN REGULATOR 

TORQUE WHEELS USED FOR THE REGULATOR 
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STOLE CONFIGURATION 


For the experiments repoorted herein and for the LOG design, the 
SCOLE platform rested on the ground and was considered fixed. We used 
the mid-mast and reflector accelerometers and the rate gyros on the mast 
tip. The actuators used were the reaction jets on the reflector and the 
torque wheel es at the mast tip. 


SCOLE CONFIGURATION 


SCOLE PLATFORM FIXED 
SENSORS — 

MID-MAST AND REFLECTOR ACCELEROMETERS 
RATE GYROS ON MAST TIP 

ACTUATORS — 

JETS ON REFLECTOR 
TORQUE WHEELS AT MAST TIP 
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SCOLE MODEL AND TEST RESULTS 


The next 10 slides are working charts organized in 5 pairs. 

They concern the 5 modes of the design model. The first chart of each pair 
contains the mode shape and frequency. This chart is followed by an 
experimental data record taken by manually exciting the structure 
at the natural frequency of the mode and taking free-decay data. The 
estimated mode amplitudes are indicated on the traces. 


0RBT I0NRL MODE. FREO (HZ) 


• H426 X 10 * 0 0 10= 1/2/1 
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future fuvjs 


Tasks that remain to be accomplished are the complete validation of 
the Kalman filter and regulator for both free-decay and forced response 
The SH3T must be tested on this nominal filter design and thresh o lds 
need to be set to avoid false alarms in light of the modelling errors 
inherent in the design. Possible sources of the modelling errors are 
excitation of modes not modelled and higher order and non linearities in 
the description of the sensors and actuators. The next step is to 

several failure cases for the ARM and generate appropriate i qq 
designs for each of these. Then the ARM performance can be evaluated. 
Current plans call for this to be completed by mid June 1988. This 
schectile is ambitious and may slip because of NA6A revectoring of 
resource*. 


FUTURE PLANS 

VALIDATE NOMINAL KALMAN FILTER 
TEST SPRT ON NOMINAL 
DESIGN FOR NULL FAILURES OF SENSORS 
VALIDATE FAILURE CASE DESIGNS 
TEST OVERALL ARM FDI PERFORMANCE 

TO BE COMPLETED BY BY MID JUNE '88 
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SUMMARY 


The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 
additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 

which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of the Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom- 
mittee on Large Space Structures. The physical apparatus will consist of a 

sof tly supported dynamic model of an antenna attached to the Shuttle by a 

flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 


34-8 


m 3, ms Homy*., 


i 



target, under conditions of noisy data, limited control authority and 
random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information intended to 
facilitate the analysis and control synthesis tasks. A workshop is planned 
for early December at the NASA Langley Research Center to discuss and 
compare results. 


INTRODUCTION 

Many future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex- 
ible modes will decrease and overlap the pointing system bandwidth. It 
will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 
stability, reliability and the pointing accuracy required for large, 
flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 
modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula- 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 
Concerns remain, however, because of the chronic difficulties in control- 
ling these lightly damped large-scale systems. Because of these concerns 
and because of the desire to offer a means of comparing technical 
approaches directly, an NASA/ IEEE Design Challenge is being offered. An 
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e*p er i«„ t8 l t..t article la being assembled under the cognisance of the 
Spacecraft Control Branch at the NASA Langley Research Center rith the 
advice and counsel of the IEEE (COLSS) Subcommittee on Large Space 
Structures. This Spacecraft Control Laboratory Experiment (SCOLE) rill 
serve as the focus of a design challenge for the purpose of comparing 
directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space Shuttle orblter by a flexible beam. This configuration uas 
chosen because of Its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper rill discuss the "Design 
Challenge" In term, of both a mathematical problem and a physical experi- 
mental apparatus. The SCOLE program Is not part of any flight program. 


SYMBOLS 


acceleration vector ft/sec^ 
beam cross section area 
observation matrix 

noise contaminating direction cosine matrix measurements 

line-of-sight error 

modulus of elasticity 

concentrated force expressions 

force vector 

concentrated moment expressions 
torsional rigidity 

moment of Inertia matrix for entire Shuttle/antenna configuration 
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moment of inertia matrix. Shuttle body 
moment of inertia matrix, reflector body 
1$ beam cross section moment of inertia, roll bending 

I© beam cross section moment of inertia, pitch bending 

Iy beam polar moment of inertia, yaw torsion 

L length of the reflector mast, beam 

Mj control moment applied to the Shuttle body 

M 4 control moment applied to the reflector body 

Mj) disturbance moment applied to the Shuttle body 

m mass of entire Shuttle/antenna configuration 

mi mass of Shuttle body 

m^ mass of reflector body 

P mass density of beam 

s beam position variable 

T. direction cosine matrix. Shuttle body () , = T ( ) 

1 3 earth l w Shut tie body 

T, direction cosine matrix, reflector body () , = T. () 

4 3 w earth 4 w reflector 

body 

vj inertial velocity, Shuttle body 

V 4 inertial velocity, reflector body 

lateral deflection of beam bending in y-z plane 
uq lateral deflection of beam bending in x-z plane 

uy angular deflection of beam twisting about z axis 

X,Y,Z position variables 

A displacement of proof-mass actuator 

6 line-of-sight pointing requirement 

e noise contaminating angular velocity measurements 
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pitch, roll, heading 
? damping ratio 

T noi8e contaminating acceleration measurements 

“1 angular velocity of Shuttle body 

“4 angular velocity of reflector body 

DISCUSSION 

objective of the NASA- IEEE Design Challenge concerning the cootrol 
of flexible epececr.ft Is to promote direct comperleon of dlffereot 
pproachee to control, state estimation and systems Identification. The 
design challenge has principal parts, the first using a mathematical model, 
and the second using laboratory experimental apparatus. The specific parts 

of the Spacecraft Control laboratory Experiment (SCOLE) program will be 
discussed in detail. 


wuiroi UDjectives 

The primary control task 1. to rapidly slew or change the ll„e-„f- 
slght of an antenna attached to the space shuttle orblter, and to settle or 

da*** SCrUCtUral “ the degree required for precise pointing 

of the antenna. The objective will be to minimise the time required to 
-lew and settle, until the antenna ll„e-of-slght remains within the 
8 a «. A secondary control task Is to change direction during the 
on-target" phase to prepare for the next slew maneuver. The objective Is 
change attitude and stabilise as quickly as possible, while keeping the 
line-of-sight error less than 6. 
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Math Model Dynamics 


The initial phase of the design challenge will use a mathematical 
model of the Shuttle orbiter/antenna configuration. It is necessary to 
obtain a balance, of course, between complex formulations which might be 
more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatics couples the otherwise uncoupled beam equations. 
Additional terras represent the action of two, 2-axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

"l - - ^Vl"! * H ! + t * Vl’ 
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Similarly, for the reflector body. 


“4 ' - I 4 1< “4 1 4“4 + M 4 * M b,4 ) 


; . F « 1 F ».« 

4 m. 


The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 


• T /v T 

T = - to T 
i ii 


• T T 

T 4 - - V 4 


The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapatlve terms which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 
end conditions and forces applied at the locations of the proof-mass 
actuators. 


ROLL BEAM BENDING: 
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PITCH BEAM BENDING: 
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YAW BEAM TORSION: 


2 3 

3 \ 3 \ 

pi y 71“ + 2 % i y ^ ~17 


- GI 


at 


a 8 at 


-2 

3 

* ~2 
3s 


n 3 *! 


g *,n {<8 


8 „> 


where: 


. « 

-* 1 1 at 2 


s=»0 


a 2 
d u . 




4 ) 2 2 2 

z 3t 


a 2 A 


+ m 


♦ ,2 


2 2 

* at 


s=s. 


(SHUTTLE BODY FORCE} 


3SS 


(PROOF-MASS ACTUATOR FORCE} 


{PROOF-MASS ACTUATOR} 


a 2 

3 u, 


3 1 


+ m, 


3 2 

J 3t 


s*s„ 


a 2 

3 u , 


3 1 


a 2 

3 u,. 


3t 


- I 


Jl 

3 % 

, — r 1 /32.5 + F 
ZZ ,4 a( .2 y 


s=130 


S=8. 


3 2 u r 


3 1 


3 2 A, 


+ m 


®±L 


{REFLECTOR BODY FORCE} 
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The angular velocity of the reflector body is related to the Shuttle body 
by: 



357 


The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


-*r + r f + 2 

1 1 

• 

1 

W 

1 ~ 2 

R I r r - V 

• r a| 


where: 


Rp is the feed location (3.75, 0 , 0 ) 

Rr is location of the center of the reflector (18.75, -32.5, 

-130) 

R a is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R A can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj , and the reflector body, T 4 , by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

T 1 T 4 which is a function of the structural deformations of the beam. 
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The line-of-sight error, e, is the angular difference between the 
target direction, given by the unit vector, and the line-of-sight 

direction in Earth axes, T^Rlqs* 

• - AXCSIN |D t X TjR^I or ARCSIN ^1^1 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body-beam-body 
configuration for "pitch” bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 
the program. It will be the analyst's decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge is to validate in the 

laboratory, the system performance of the more promising control system 
designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten- 

sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 

effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 


earlier. 
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Design Challenge, Part One 

For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 
configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, 6. During the time that the LOS error is 
within 5, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan- 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 

Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

The initial line— of— sight error is 20 degrees. 

e(o) = 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 

moment of inertia in roll. The Shuttle/antenna is at rest 
initially. 

4. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 

<5 * 0.02 degree 

3G0 


possible. 



Control moments can be applied at 100 Hz sampling rate to both the 
Shuttle and reflector bodies of 10,000 ft-lb for each axis. The 
commanded moment for each axis is limited to 10,000 ft-lb. The 
actual control moment's response to the commanded value is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the Shuttle body: 


-10 < 




. < 10 
command — 


- 0.1 


- 0 . 1 , 


*St,l <n + *> - • V,U) +CI V, co ». aIld (n) 


Equations for other axes and for the reflector body are similar. 
Control forces can be applied at the center of the reflector in 
the X and Y directions only. The commanded force in a 
particular direction is limited to 800 lbs. The actual control 
force's response to the commanded value is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


-800 < F v . < 800 

— Y , command — 


F y (n + 


1 ) 


e’ 0 * 1 F y (n) 


+ (1 - e" 0 * 1 ) 


F v ,(n) 

Y , command 


Equations for X-axis are similar. 

Control forces using two proof-mass actuators (each having both 
X and Y axes) can be applied at two points on the beam. The 
strokes are limited to ± 1 ft, and the masses weight 10 lbs each. 
The actual stroke follows a first-order response to limited 
commanded values. 
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For the X-axis of the proof-mass actuator at 


s 2 : 


-1 < A <1 

— X , 2 , command — 


A v 0 (n + 1) = e -0 * 1 A (n) + (1 - e' 0 * 1 ) A y .( n ) 

A > * X , 2 , command 


X,2 


Equations for other axes and locations are similar. 

8. The inertial attitude direciton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 
rate of 100 samples per second. Each element of the direction 
cosine measurement raatarix is contaminated by additive, 
uncorrelated Gaussian noise having an rms value of 0.001. The 
noise has zero mean* 



d H (») 

d , 2 ( n) 

d 13 (n> 

T .(n + 1) = T (n) + 

s, measured s,true v ' 

d 2 ,(») 

d 22 (») 

d 23 (n) 

where : 

V n> 

d 32 (n) 

d 33 (n> 

e{ d„ (n) } = 0 

E{d ij (n)d kL (n)} = 0 

for i 

* k or 

j * L 

E{dij(n)d (n + k)} = 0 

for k 

* 0 


= [ .001 ] 2 

for k 

= 0 



3G2 



9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 

measurement is contaminated by additive, Gaussian, uncorrelated 
noise having an rms value of 0.02 degree per second. The noise 
has zero mean. 

For example: 

W 1 ,X, measured^" + ^ “l ,X, f iltered^ + E 1,X^ 

^ E 1 X (n) E 1 X (n + = ® for k * 0 

■ (.02)^ for k = 0 

where 

oj * — 20 0) + 20 ai 

1,X, filtered l,X*filtered l,X,true 

10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 

axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 

through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example: 


a i,X, measured* ” + ° = 3 1 ,X, f iltered (n) + T l,X (n) 

e{t (n) x (n + k)} = 0 for k * 0 

1 1 

= (.05) 2 for k = 0 

where : 

3 1,X, filtered " 20 3 i ,X, filtered + 20 “l.X.true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 

have zero mean. 

For example: 

e{m d x (n) M D yn + k)} = 0 for k * 0 

= (100) 2 for k * 0 

In summary, the designer's task for part one is to: (1) derive a 

control law for slewing and stabilization, coded in FORTRAN; (2) select an 
Initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 
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Design Challenge, Part TVo 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 
generation capability will apply to part two, as for part one. Again, the 
analyst may select the initial alignment about the assigned initial RF 
line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer’s task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 
comparing directly different approach to controlling a flexible 
Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment (SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/ IEEE Design Challenge has the advice and counsel of the 
IEEE-COLSS Subcommittee on Large Space Structures. Workshops will be held 
to enable investigators to compare results of their research. 


3G3 


MASS CHARACTERISTICS 



205,800 1 , 132,508 7 , 007,447 7 , 113,962 - 7,555 115,202 












































The moment of inertia becomes: 
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m - 6391 .30 slugs 
nij= 6366.46 slugs 

m 2 = 0.3108 slugs 

m^= 0.3108 slugs 

m = 12.42 slugs 
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PA 
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EI n 

= A.O x 10 7 lb-ft 2 

El* * 

4.0 x 10 lb-ft 

0 

= .003 

C . 

.003 

S 0 



PI . * 

0.9089 slug-ft 

- 

4.0 x 10 7 lb-f t 2 


.003 


003 
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CONTROL MOMENTS, DISTURBANCES 






Roll bending mode # 1. Frequency= .32 Hz 



Roll bending mode # 2. Frequency= 1.29 Hz 



Roll bending mode # 3. Frequency= 4.80 Hz 



Roll bending mode # 4. Frequency= 12.29 Hz 



Roll bending mode | 5. Frequency^ 23.68 Hz 



Roll bending mode # 6. Frequency^ 38.89 Hz 



Roll bending mode # 7. Frequency— 57.90 Hz 



Roll bending mode # 8. Frequency= 80.72 Hz 


Figure 4;i.- Plots of normalized roll bending mode shapes 
for SCOLE configuration. 
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AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE 


Pitch bending mode # 1. Frequency= .29 Hz 



Pitch bending mode # 2. Frequency= 1.65 Hz 



Pitch bending mode # 3. Frequency= 4.97 Hz 



Hz 


Hz 


Pitch bending mode # 6. Frequency= 38.91 Hz 



Pitch bending mode # 7. Frequency= 57.92 Hz 



Pitch bending mode # 8. Frequency= 80.73 Hz 


Figure 4b.- Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Toraional mode # 1. Frequency= .53 Hr 



Torsional mode # 2. Frequency = 45.12 Hz 



Torsional mode # 3. Frequency^ 90.23 Hz 



Torsional mode # 4 Frequency^ 135.35 Hz 


« .snp- 

5 > 


t-L 


Torsional mode # 5. Frequency^ 180.46 Hz 



Toraional mode # 6. *requency= 225.57 Hz 


I 

8 

l 



S-.*L 

Torsional mode # 7 


£ 

* 


Frequency= 270.69 Hz 


§ 27 E 



Toraional mode # 8. 



£ 

* 


Frequency= 315.80 Hz 


Figure 4c.- Plots of normalized torsional mode shapes for 
SCOLF. configuration. 
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CONTROL DESIGN CHALLENGES of LARGE SPACE SYSTEMS and 
SPACECRAFT CONTROL LABORATORY EXPERIMENT (SCOLE) 

Gene Lin of Control Research Corporation 

Examines the Resulting Excitation due to Bang-Bang Slewing 
Effect of Excitation on Line-of-Sight Error 
Concept of Modal Dashpots 

Examines the Resulting Excitation Employing Modal Dashpot Design 
Examines Linear Velocity Feedback Force Control (Also Angular) 
Computer Simulation 


Needs to be Applied to SCOLE Experimental Apparatus 
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INITIAL TEST RESULTS on STATE ESTIMATION on the SCOLE MAST 


Dean Sparks of NASA Langley Research Center 

Modal State Estimation Tested on SCOLE (Fixed Shuttle Body) 

Kalman Filter used for State Estimation 

Six Linear Accelerometers and 3-Axis Rate Gyro 

Large Discrepencies Between Linear Finite Element Model and 
Experiment 

Sinusoidal Input Forcing Function 


Gravity Effects Very Pronounced 

A Large Number of Modes Required for Accurate Static Deflection 
Nonlinear Kinematics Can Be Significant 
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The DYNAMICS and CONTROL of the IN-ORBIT SCORE CONFIGURATION 


MINIMUM TIME ATTITUDE SLEWING MANEUVERS of a RIGID 
SPACECRAFT 

STABILITY ANALYSIS of LARGE SPACE STRUCTURE CONTROL SYSTEMS 
with DELAYED INPUT 

Peter Bainum, A.S.S.R. Reddy, Cheick M. Diarra and Feiyue Li of Howard 
University 

Examines the Changes in Modal Characteristic due to Orbital Motion 
Linear Elasticity Assumed 
Retain Nonrotating Mode Shapes 

Examines Stability for Rigid Case with Increasing Complexity 
Derives Control Law using Linear, Flexible (4 Modes) SCOLE 
Computer Simulations 


Need to Introduce Disturbances, Noise 
Need to Apply to Experimental Apparatus 
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CONTROL DESIGN APPROCHES for LaRC EXPERIMENTS 
Steve Yurkovich and Umit Ozguner of Ohio State University 

LQG, MEOP 

Computer Simulations 

Model Reference Adaptive Control 

Employ Hyperstability and Positivity Concepts 

Combine with Parameter Identification 

Computer Simulations 


Noise and Parameter Uncertainity Models Not Necessarily Realistic 
Need to Apply to Experimental Apparatus 
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SOME NONLINEAR DAMPING MODELS in FLEXIBLE STRUCTURES 


A. V. Balakrishnan of UCLA 

Nonlinear Damping Term ie lxl a ldx/dtl b dx/dt 
Uses Krylov and Bogoliubov Approximation 
"Solves" Equations for SCOLE Bending and Torsion 
Examines Multivariable Case 

Draws Analogy with Nonlinear, Boundary Feedback Control 


Useful Results for Modeling Nonlinear Damping 
Useful for Systems Identification 
Experimental Evidence Supports Findings 
Need Approximations for Modal Model 


380 



SOME NONLINEAR DAMPING MODELS in FLEXIBLE STRUCTURES 


A. V. Balakrishnan of UCLA 

Nonlinear Damping Term ie lxl a ldx/dtl b dx/dt 
Uses Krylov and Bogoliubov Approximation 
"Solves" Equations for SCOLE Bending and Torsion 
Examines Multivariable Case 

Draws Analogy with Nonlinear, Boundary Feedback Control 


Useful Results for Modeling Nonlinear Damping 
Useful for Systems Identification 
Experimental Evidence Supports Findings 
Need Approximations for Modal Model 
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INFINITE-DIMENSIONAL APPROACH to SYSTEMS IDENTIFICATION 


of SPACE CONTROL LABORATORY EXPERIMENT (SCOLE) 
S. A. Hossain and K. Y. Lee of Penn State University 

Retains Physical Parameters ie Coefficients in P.D.E.s 
Truncation problem avoided 
Weighted Least Squares 

Computer Simulations for a Few Different Examples 
Convergent after Several Iterations 

Additive Noise Used is Not Particularly Realistic 
Good to tie-in to Physical Parameters 
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OPTIMIZATION-BASED DESIGN of CONTROL SYSTEMS for FLEXIBLE 
STRUCTURES 

E. Polak, T. Baker, T-L. Wuu and Y-P. Ham of University of California at 
Berkeley 

Approach is to Use "Nonsmooth Optimization Algorithms" 

Design Closed-Loop Finite Dimensional Compensators for Flexible 
Structures 

Does Not Require Modal Decomposition, Avoiding Spillover Problem 

Applicable to Mixed Ordinary and P.D.E. Systems 

Classical Design Objectives 

Applied to Flexible Beam with Tip Mass 

Computer Simulation 


Need to Test with Disturbances, Noise 
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PLACING DYNAMIC SENSORS and ACTUATORS on FLEXIBLE SPACE 
STRUCTURES 

Gregory Norris and Robert Skelton of Purdue University 


Input/Output Cost Analysis Decomposes the Quadratic Cost into 
Contributions from Each Stochastic Input and Weighted Output. 

Closed Loop is Considered 

Past Treatment Required Perfect Sensors and Actuators 
Now Consider Dynamics of Sensors and Actuators 
Applied to SCOLE 


Important to Consider Closed Loop and Dynamics of Elements 
Order Reduction is an Important Problem 
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COMBINED PROBLEM of SLEW MANEUVER CONTROL and VIBRATION 
SUPPRESSION 

Y. P. Kakad of University of North Carolina at Charlotte 


Derives the Nonlinear Equations of Motion Suitable for Slewing and 
Vibration Control 

Analyzes the Kinematic Nonlinearities in terms of the Matrix Spectral 
Norm 

Analyzes Arbitrary Slew Maneuver 

It May be Necessary to Synthesize Control Systems for the Combined 
Problem 


Interesting Approach 

Determining "Optimal" Control Might be Difficult 
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EFFECT of ACTUATOR DYNAMICS on CONTROL of BEAM FLEXURE DURING 
NONLINEAR SLEW of SCOLE MODEL 

Shalom (Mike) Fisher of Naval Research Laboratory 

Simulation Includes Limited Proof-Mass Deflection and Time Delay 
Examines Line-of-Sight Errors due to Slewing and Settling 
NASTRAN Finite Element Model (12 Vibrational Modes) 

Nonlinear DISCOS Simulation of 20' Slew 
LQR Design 
Force Limit 
Computer Simulation 

Ready for Experimental Apparatus 
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COMPUTATIONAL EXPERIMENTS in the OPTIMAL SLEWING of FLEXIBLE 
STRUCTURES 


T. E. Baker and E. Polak 

Considers the "Swinging Flexible Arm Problem" 

45 Degree Rotation 
Torque Limit 

Minimizes Total Energy Required for Time Limit 
Same with Upper Bound on Potential (Elastic) Energy 
Computer Simulation 


Need to Consider Disturbances, Noise 


38 7 


PARAMETER IDENTIFICATION USING MODAL DATA 


L. Meirovitch and M. A. Norris of V.P.I. 

Uses Indirect Method ie Modal "Data" not Sensor Data 

Rayleigh-Ritz method 

Iterative, Sensitivity Analysis 

SCOLE Laboratory Experiment Application 

4 Degree-of-Freedom Model 

Estimated El, I(Antenna) 

Work at early stage. 

Important coupling not yet included. 
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NONLINEARITIES in SPACECRAFT STRUCTURAL DYNAMICS 

Lawrence W. Taylor and Kelly Latimer of NASA Langley Research 
Center 

Many Modes are Needed for Accurate Static and Dynamic Characteristics 

Large Amplitude Deflection Effects are Examined 

Lumped Mass and Assymptotic Approximations are Evaluated 

A 3-Dimensional Beam Equation is Derived for Large Deflections 

Linear and Nonlinear Damping Models are Examined 

Distributed Parameter Models are Seen To Reduce the Number of 
Parameters 


The Dynamic Model of the SCOLE Experimental Apparatus Needs to be 
Improved 

Effort is Needed in Solving the 3 -Dimensional Beam Equation 

The Nonlinear Damping Seen in the SCOLE is Probably Due to 
Aerodynamics 
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ANALYTICAL REDUNDANCY MANAGEMENT for SCOLE 


Ray Montgomery of NASA Langley Research Center 


Examines Problem of Control Performance when Elements are Failed 
Uses Expected Quadratic Performance Index 
Combinatorial Problem forces Cut and Try 


Important Consideration for Some Applications with Long Life/MTBF 
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